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INTRODUCTION 


I . 

The  study  of  resonance  region  electromagnetic 
interaction  with  arbitrarily  shaped  and  often  inhomogeneous 
penetrable  objects  has  increased  dramatically  in  the  past 
few  years.  This  new  concentration  of  effort  is  driven  by 
the  likelihood  of  expanded  implementation  of  high  frequency 
over-the-horizon  radars,  very  high  frequency  satellite-borne 
radars  and  medium  frequency  surface  wave  radars.  Such 
radars  have  the  potential  to  observe  "stealth"  type 
platforms  whose  low  radar  cross  section  (RCS)  design  was 
limited  to  more  conventional  radar  frequencies,  above  500 
MHz.  In  addition,  these  new  studies  of  electromagnetic 
interaction  have  proven  to  be  of  benefit  in  other  areas 
such  as  propagation  through  rain,  power  absorption  by 
biological  tissues,  focusing  by  microwave  lenses,  radiation 
from  reentry  vehicles  and  electromagnetic  pulse  (EMP) 
coupling  into  vehicles,  installations  and  weapon  systems. 
[Refs.  1  and  2] 

As  frequencies  become  sufficiently  low,  structures 
having  a  low  radar  cross  section  may  exhibit  strong 
resonance  characteristics.  Design  of  such  structures  is 
complicated  by  these  low  frequency  resonance 
characteristics,  due  to  the  lack  of  accurate  analytical 
methods.  Even  attempts  at  structure  design  utilizing 
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"shaping"  and  "edge  treatments"  are  ineffective  at  providing 
a  highly  reduced  radar  cross  section  at  resonant 
frequencies.  At  higher  frequencies,  "asymptotic" 

formulations  are  available  to  rectify  these  problems  [Ref. 
3].  A  more  comprehensive  approach  is  needed,  one  which 
allows  the  designer  to  rigorously  evaluate  the  interaction 
of  the  scatterer  with  the  illuminating  field.  At  the  same 
time,  such  a  tool  should  be  highly  interactive  with  the 
user. 

The  subject  of  this  thesis  is  the  development  of  a 
generalized,  interactive,  finite  element  scattering 
algorithm  to  be  used  in  the  analysis  and  computer  aided 
design  of  scattering  structures.  In  this,  the  first  stage 
of  effort,  the  programs  developed  for  use  on  a  personal 
computer  calculate  electromagnetic  scattering  by 
axisymmetric  inhomogeneous  structures.  These  structures  are 
composed  of  user  specified  lossy  or  lossless  penetrable 
(i.e.,  dielectric  or  magnetic)  materials.  The  resulting 
programs  provide  significant  improvement  in  computational 
speed  and  versatility  over  the  current  integral  equation 
based  computer  codes  being  used.  The  physical  structure  can 
be  quite  general,  other  than  the  requirement  for  rotational 
symmetry.  This  general  construction  avoids  the 
computational  restrictions  found  in  surface  integral 
equation  formulations,  where,  as  frequency  and  geometric 
complexity  increase,  the  time  required  to  obtain  a  solution 
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increases.  The  finite  element  algorithms,  based  upon  a 
volume  variational  formulation,  are  relatively  immune  to 
structural  complexity  driving  the  computation  time.  Another 
advantage  of  these  algorithms  is  their  method  of  matrix 
manipulation.  The  finite  element  algorithms  produce  very 
"sparse"  matrix  structures,  whereas  the  integral  equation 
formulations  produce  "full"  matrices  [Ref.  4].  Finally, 
when  considering  the  construction  of  realistic  material 
configurations,  the  computational  time  required  utilizing 
the  finite  element  algorithms  can  be  substantially  less  than 
the  time  required  by  the  use  of  integral  equation 
formulations  [Ref.  5]. 

The  development  of  the  axisymmetric  finite  element 
program  is  based  upon  Morgan's  previous  effort  [Ref.  1] , 
which  restricted  structures  to  non-metallic  configurations. 
Major  programming  modifications  to  improve  computational 
efficiency  and  to  incorporate  computer  aided  design  features 
such  as  interactive  graphics  are  now  integral  parts  of  the 
finite  element  computer  program. 

A  brief  overview  of  the  finite  element  scattering 
algorithm,  which  constitutes  the  basis  of  this  thesis,  is 
presented  in  Chapter  II.  In  addition  to  the  appropriate 
background  information,  the  concept  of  coupled  azimuthal 
potentials,  the  finite  element  algorithm  and  the  unimoment 
method  are  reviewed. 
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Having  established,  ir.  Chapter  IX,  the  necessary  theory 
utilized  in  the  construction  of  the  original  version  of  the 
computer  program  EMSCAT,  Chapter  III  presents  methods  to 
construct  program  input  data.  Included  in  this  chapter  are 
discussions  of  computer  programs  written  specifically  for  a 
given  geometric  shape,  data  file  input  and  keyboard  input. 
Basic  flow  charts  provide  an  overview  of  the  input-output 
interactions  of  the  all  computer  programs. 

The  first  half  of  Chapter  IV  discusses  the  extension  of 
the  methods  used  to  construct  program  input  data  to  those 
methods  that  utilize  graphics.  Of  special  interest  is  the 
incorporation  of  computer  aided  design  features.  The  second 
half  of  this  chapter  addresses  the  user  interactive  output 
graphics  features  provided.  These  features  allow  the  user 
to  obtain  one  of  four  different  graphs,  either  on  the 
computer  screen  or  printer,  of  the  EMCAD  output  data  of 
choice. 

In  Chapter  V,  the  validation  of  modifications  and 
improvements  to  EMCAD  is  presented.  Graphical  comparisons 
are  made  between  the  output  of  the  modified  version  of 
EMSCAT  versus  both  theoretical  calculations  and  actual 
experimental  data.  In  addition,  validation  of  new  programs 
written  specifically  to  be  used  in  conjunction  with  EMCAD  is 
presented. 

The  brief  conclusions  in  Chapter  VI  will  summarize  the 
work  previously  presented,  as  well  as  provide 
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recommendations  for  future  modifications  to  the  interactive 
package  of  computer  programs  being  used. 

The  Appendices  provide  copies  of  the  source  code  for 
all  the  computer  programs  utilized  in  this  interactive 
package.  A  sample  listing  of  the  output  of  EMCAD  is 
attached.  As  an  added  feature,  user  manuals  with  hints, 
warnings  and  directions  for  all  of  the  major  computer 
programs  are  included. 
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II.  FINITE  ELEMENT  SCATTERING  ALGORITHM 


A.  BACKGROUND 

In  recent  years,  great  attention  and  effort  has  been 
directed  toward  the  attempt  to  obtain  the  numerical 
electromagnetic  scattering  solutions  for  inhomogeneous 
material  structures.  Early  efforts  dealt  primarily  with 
improving  the  efficiency  of  existing  computational 
algorithms.  These  algorithms  utilized  surface  integral 
equation  formulations,  which  required  the  inversion  of  large 
"full"  matrices.  Accordingly,  methods  using  these 
formulations  were  computationally  slow  while  being  highly 
restricted  in  the  level  of  inhomogeneity  that  could  be 
handled.  [Refs.  1  and  2] 

A  major  advance  in  matrix  manipulation  occurred  with 
the  development  by  Mei  [Ref.  6]  of  the  unimoment  method. 
This  new  method  decoupled  the  interior  and  exterior  regions 
associated  with  radiation  or  scattering  problems,  resulting 
in  uniformly  banded  system  matrices  for  the  interior  region 
solution.  The  primary  advantage  of  uniformly  banded 
matrices  is  the  capability  to  invert  them  directly  using 
fast  banded  matrix  inversion  techniques,  such  as  the  Riccati 
transform.  Mei's  breakthrough  opened  the  door  for  the 
development  of  applications  such  as  the  finite-difference 
solution  of  radiation  by  a  biconical  antenna,  and  the 
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finite-element  solution  of  scattering  by  two-dimensional 
inhomogeneous  dielectric  cylinders.  [Refs.  1  and  2] 

Attempting  to  improve  upon  these  applications,  Morgan, 
Chang  and  Mei  [Ref.  7]  conducted  research  with  the  goal  of 
extending  the  unimoment  technique  to  problems  involving 
three-dimensional,  axisymmetric  inhomogeneous  structures. 
During  the  initial  period  of  this  research,  another  major 
breakthrough  occurred.  This  breakthrough  was  the  discovery 
of  the  coupled  azimuthal  potential  (CAP)  formulation 
[Ref.  7].  The  CAP  formulation  became  the  basis  for 
construction  of  several  numerical  algorithms.  These 
algorithms  incorporated  the  capability  to  perform  scattering 
computations  for  inhomogenous  axisymmetric  metallic  and 
penetrable  structures  [Ref.  2). 

Utilizing  a  triregional  unimoment  method,  these 
numerical  algorithms  were  incorporated  into  a  versatile 
finite-element  computer  program.  The  triregional  unimoment 
method  was  chosen  because  of  its  ability  to  minimize  the 
bandwidth  of  the  global  system  matrix,  and  its  ability  to 
minimize  the  number  of  nodes  in  the  finite-element  mesh. 
The  finite-element  computer  program  generates  scattering 
solutions  for  an  arbitrarily  shaped  inhomogeneous  lossy 
penetrable  body  of  revolution  being  impinged  upon  by 
multiple  incident  fields.  [Ref.  2] 
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B.  COUPLED  AZIMUTHAL  POTENTIALS 

The  solution  of  electromagnetic  field  problems 
involving  generalized  axisymmetric  penetrable  structures 
using  the  classical  TE  and  TM  potential  formulations  are 
restricted  to  a  handful  of  configurations,  and  have  at  most 
inhomogeneities  in  one-dimension.  By  utilizing  the  coupled 
azimuthal  potential  (CAP)  formulation,  solutions  are  valid 
in  lossy  inhomogeneous  isotropic  rotationally  symmetric 
media,  as  shown  in  Figure  2-1,  where  the  complex  material 
parameters  are  invariant  to  the  azimuthal  coordinate. 
[Ref.  1] 

Utilizing  normalized  circular  cylindrical  coordinates, 
as  also  shown  in  Figure  2-1,  the  first  step  in  developing 
the  necessary  eguations  is  to  decompose  the  total 
electromagnetic  field  into  azimuthal  modes  through  the  use 
of  an  exponential  Fourier  series, 
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where 


■:  R  i  Z  i  f  X  =  (  k  Q  p  ,  k  Q  z  »  i ) 


(2.3) 


with 


(2.4) 


and 


ng  =  1  2 0 tt  52 


(2.5) 


Two  uniformly  continuous  coupled  azimuthal  potential 
functions 


t  (  R  ,  Z  >  m  )  and  V  <  R i Z  j  ra ) 
1  2 
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can  be  used  to  generate  the  modal  electromagnetic  field 
vectors , 


e  (R,Z) 
m 


J f 

m  1 


R>  f 
r 


X 


f  — 

R 


(2.6) 


h  ( R ,  Z )  =  jf  (■**,  -  Re  *  x  9>|r  )  +  f  —  ,  (2.7) 


where  the  two-dimensional  gradient  operator  is 


v  = 


,5  A  a 

R  -  +  Z  - 

9R  91 


(2.8) 


while  the  multiplicative  function  is  given  by  [Refs.  1,  2, 
and  7  ] 


t  ( R  i  Z )  =  In  <  R )  Z )  e  <  R »  Z )  R  2  -  »2]'1 
in  r  r 


(2.9' 
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For  a  two-dimensional  cut  in  the  (R,Z)  plane  of  Figure 
2-1,  the  coupled  azimuthal  potential  functions  satisfy  a 
formally  self-adjoint  system  of  second-order  partial 
differential  equations.  Utilizing  either  Euler-Lagrange  or 
stationary-theorem  approaches,  it  can  be  shown  that  there  is 
an  easily  managed  variational  principle  which  can  be  applied 
in  addition  to  the  partial  differential  equation 
formulation. 

This  variational  formulation  is  based  upon  the 
existence  of  a  stationary  function  which  manifests  itself  as 
a  surface  integral  over  the  planar  cross  section,  S, 

1  p 

F  =  1<R,  .  *2*^  i  ^2>  dR  dZ  (2.10) 

J  S 

where  the  Lagrangian  L,  which  is  unique  to  within  an 
arbitrary  constant  multiplier  and  arbitrary  independent 
additive  function,  is  of  the  form  [Refs.  1  and  2] 
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L  =  ?'  •  (Re  vy  +  m#  x  vy  )  + 

mil  r  1  2 


•  (R>i  -  in#  x  > 

2  r  2  1 


-  [  £  *?  ♦  *  *1  1  .  (2.11) 

R  L  r  1  r  2  J 

C.  FINITE-ELEMENT  ALGORITHM 

Having  presented  an  overview  of  the  coupled  azimuthal 
potential  formulation,  it  is  now  appropriate  to  explain  its 
application  to  the  electromagnetic  scattering  problem. 
Figure  2-2  illustrates  the  finite-element  mesh  geometry 
employed  by  the  numerical  algorithm.  The  (R,Z)  planar 
configuration  shown  is  named  the  semi-annular  conformal 
(SAC)  mesh  [Ref.  1]  .  This  mesh  is  comprised  of  numerous 
linear  triangular  elements,  which  conform  to  the  surface  of 
the  axisymmetric  scattering  body.  This  surface  conformation 
simplifies  analysis  of  normal  derivative  discontinuities  at 
the  surface  boundary.  The  semi-annular  conformal  mesh 
provides  a  distinct  computational  advantage  to  those 
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FIG.  2-2.  SEMI-ANNULAR  CONFORMAL 
FINITE-ELEMENT  MESH  FOR  AN 
INHOMOGENEOUS  BODY  OF  REVOLUTION 
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scattering  structures  which  can  support  a  fairly  large 
homogeneous  spherical  core  [Ref.  2]. 

Figure  2-2  graphically  demonstrates  the  regions  of 
interest  when  calculating  the  solution  of  the  scattering 
problem.  The  homogeneous  spherical  core  occupies  Region  I, 
the  inhomogeneous  cross  sectional  surface  of  the 
axisymmetric  penetrable  scattering  body  is  in  Region  II,  and 
an  unbounded  exterior  region  is  in  Region  III.  The  surface 
of  interest  is  shown  with  the  capability  of  being  both 
banded  and  layered.  The  layers  are  depicted  as  "films"  that 
run  the  entire  length  of  the  surface,  while  bands  are 
depicted  as  isolated  "pockets".  The  possibility  exists,  as 
shown,  that  all  the  bands  and  all  the  layers  have  differing 
material  parameters,  thus  establishing  a  structure  that  is 
truly  inhomogeneous.  Note  that  the  surfaces  of  Region  I  and 
Region  III  bound  the  surface  of  interest  located  in  Region 
II.  These  surfaces,  r  =  a  and  r  =  b,  can  therefore  be  used 
to  specify  Dirichlet  boundary  conditions  [Ref.  2].  The 
numerical  solution  for  the  electromagnetic  fields  in  the 
finite-element  domain  (Region  II)  can  be  computed  for  each 
of  a  set  of  boundary  conditions.  By  requiring  continuity  of 
the  electromagnetic  fields  at  the  interfaces  where  r  =  r^ 
and  r  =  r2  ,  the  unknown  expansion  coefficients  for  the  core 
and  scattered  fields  can  be  obtained  [Ref.  1].  This  is  the 
essential  concept  behind  the  unimoment  method. 


In  Region  II  of  Figure  2-2,  the  semi-annular  conformal 
mesh  is  composed  of  nodes,  with  each  interior  node 
surrounded  by  six  triangular  regions,  known  as  "elements". 
The  azimuthal  field  components  of  the  mesh  are  described  by 

N 

§  ■  e  (R,  Z)  =  7  e  <n)u  <R, 2>  (2.12) 

tn  u  m  n 

n  =  1 

and 

H 

§  •  h  (R,Z)  =  7  h  ( n ) u  <R,Z)  (2.13) 

m  ^  m  n 

n=  1 

where  un(R,Z)  represents  the  area  coordinate  linear  pyramid 
basis  functions.  Each  of  these  basis  functions  has  linear 
variation  in  the  six  triangular  elements  around  each  node, 
while  maintaining  a  value  of  one  at  the  central  node  and  a 
value  of  zero  at  all  bounding  nodes.  Each  set  of  em(n)  and 
hm(n)  represent  one  of  the  N  complex  nodal  values  found  in 
the  azimuthal  field.  To  initiate  the  variational  finite- 
element  technique,  the  linear  shape-function  expansions  in 
Equations  2.12  and  2.13  are  substituted  into  the  surface 
integral  stationary  function  F  described  in  Equation  2.10. 
Imposing  the  null  derivative  condition  with  respect  to  each 
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of  the  unknown  coordinates  representing  the  interior  nodal 
fields,  while  constraining  the  coordinates  of  the  boundary 
nodal  fields,  allows  the  computation  of  the  stationary 
point  of  the  function  F.  The  resulting  linear  system  for 
the  nodal  fields  is  described  by 
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where  n  =  each  internal  node.  When  these  equations  are 
combined  with  the  semi-annular  conformal  finite  element  mesh 
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analysis ,  each  equation,  for  a  given  interior  node  n,  will 
relate  em(n)  and  hm(n)  to  14  adjacent  node  azimuthal  field 
values  em(i)  and  hm(i).  Some  of  these  adjacent  field  values 
may  be  specified  on  a  boundary.  The  result  of  conducting 
these  calculations  for  each  value  of  n  is  a  sparse  global 
system  matrix,  containing  complex  elements  and  displaying  a 
diagonal  tri-block  submatrix  structure*  Applying  the 
Dirichlet  boundary  conditions  to  the  semi-annular  conformal 
finite-element  mesh  shown  in  Figure  2-2,  the  numerical 
solution  for  the  electromagnetic  fields  in  Region  II  is 
obtained  by  utilizing  a  two-sweep  block -by-block  matrix 
inversion  algorithm  related  to  the  Riccati  transform  for 
difference  equations.  The  diagonal  tri-block  submatrix 
structure  is  given  by 


k  k-  1  .  k  k  kk+1  k 


(2.16) 


A^,  Bfc,  and  are  sparse  submatrices  which  constitute 
the  tri-block  diagonal  structure  of  the  global  system 
matrix.  *k  is  generated  from  Equations  2.14  and  2.15 
through  the  use  of  the  Dirichlet  boundary  conditions. 
[Ref.  2] 
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The  mesh  construction  and  matrix  generation  are 
conducted  on  a  local  basis  in  a  sequential  manner.  The 
analysis  starts  at  the  positive  Z-axis,  and  works  clockwise 
in  equal  ©-increments.  Using  a  simple  minimum  distortion 
scheme  at  each  step,  the  mesh  is  constructed  to  conform  to 
multiple  approximate  piecewise  smooth  surfaces.  The  Riccati 
transform  sparse  matrix  inversion  technique  conducts  a 
"first  sweep"  and  a  "backsweep".  On  the  "first  sweep",  the 
values  of  R^  and  are  calculated,  where 
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On  the  "backsweep",  the  values  of  the  solution  vector  are 
obtained,  where  the  solution  vector  is  given  by  [Refs.  1  and 
2] 
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A  principal  advantage  of  this  numerical  method  is  that 
the  Riccati  transform  technique  provides  a  systematic 
procedure.  This  procedure  divides  a  larger  area  into 
several  smaller  areas  defined  by  the  tri-block  matrix 
substructure,  existing  between  any  three  adjacent  radial 
rows  of  nodal  values.  Analysis  using  this  type  of  smaller 
area  procedure  allows  all  operations  including  mesh 
construction,  matrix  element  generation  and  matrix 
inversion  to  be  performed  at  the  same  time.  [Ref.  2] 

D.  UNIMOMENT  METHOD 

The  unimoment  method  is  a  computational  procedure  which 
allows  numerical  values  to  be  calculated  for  the  complex 
coefficients  of  the  incident,  scattered  and  spherical  core 
electromagnetic  field  expansions  when  dealing  with  the  semi- 
annular  conformal  finite-element  mesh.  The  scattering 
solution  is  obtained  by  representing  the  electromagnetic 
fields  of  Region  I  and  Region  III,  as  shown  in  Figure  2-2, 
in  terms  of  TE  and  TM  spherical  harmonic  expansions,  with 
unknown  core  and  scattered  field  coefficients.  Using 
superposition  of  boundary  conditions  and  interior  solutions 
for  the  azimuthal  fields,  the  core  expansion  is  established 
as  a  boundary  condition  along  r  =  a  while  the  scattered 
field  expansion  added  to  the  incident  field  is  established 
as  a  boundary  condition  along  r  =  b.  The  two  existing  sets 
of  boundary  conditions  are  modal  spherical  harmonic  core 
fields  '  along  r  =  a,  with  zero  fields  along  r  =  b,  and 
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spherical  harmonic  and  incident  fields  at  r  =  b,  with  null 
fields  at  r  =  a.  The  multiple  boundary  conditions  will 
generate  multiple  solutions  at  the  nodes  where  r  =  r^  or  r  = 
r2 ,  (i.e., along  the  boundaries).  The  method  used  to  obtain 
the  expansion  coefficients  sets  the  weighted  superposition 
of  interior  solutions  to  the  appropriate  original  field 
expansions.  Along  the  boundaries  at  r  *  and  r  =  r2 ,  this 
method  enforces  equality  utilizing  least  squares  analysis. 
[Refs.  1  and  2] 

This  method  of  analysis  is  only  valid  if  the  scattered 
electromagnetic  fields  everywhere  outside  of  the  interior 
mesh  region  can  be  represented  by  field  harmonics.  These 
representations  appear  as  infinite  series,  and  must  be 
mathematically  truncated  to  facilitate  numerical 
calculations.  Decomposing  the  plane-wave  incident  -fields 
into  azimuthal  modes  yields 
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where  em  =  1  if  m  =  0,  em  =  2  if  m  >  1 ;  M  is  the  truncation 
index,  and  the  polarization  convention  (TE/TM)  , depicted  in 
Figure  2-3,  is  used  in  the  bracketed  sinusoids.  [Ref.  2] 

During  computations,  the  numerical  scattering  solution 
algorithm  is  run  concurrently  with  the  sequential  Riccati 
transform  finite-element  solution  in  Region  II.  On  the 
forward  sweep,  boundary  conditions  for  both  the  incident  and 
modal  fields  are  generated  at  each  node  row  being  used  in 
that  particular  sweep  step.  On  the  backsweep,  boundary 
conditions  along  the  r  =  r^  and  r  =  r2  contours  are  used  to 
compute  each  nodal  value.  The  incident  field  and  spherical 
harmonic  field  for  each  node  is  formed,  and  used  to  generate 
the  nodal  residual  fields.  During  the  backsweep,  inner 
product  integrations  are  formed,  and  these  inner  product 
results  are  added  to  the  field  moment  matrix,  T,  and  the 
incident  field  driving  vector,  D  corresponding  to  each 
incident  field.  The  scattering  solution  is  obtained  by 
multiplying  the  inverse  of  the  field-moment  matrix  with  the 
driving  vector  to  obtain  the  field  expansion  coefficients. 
The  complex  scattering  amplitude  f(9,<S)  is  defined  in  terms 
of  the  far-zone  E-field, 
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(X'.Y'.Z')  COORDINATES  WITH  RESPECT  TO  INCIDENT  FIELD 
FIG.  2-3.  POLARIZATION  CONVENTION 


where  dmax  is  the  maximum  scattering  body  dimension  and 
where  the  components  of  f(e,$)  are  given  by 
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The  scattering  cross  section  is  found  to  be  [Refs.  1  and  2] 
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Utilizing  the  theoretical  principles  and  mathematical 
calculations  just  discussed,  Morgan  [Ref.  1]  wrote  the 
computer  program  EMSCAT  to  compute  electromagnetic  field 
scattering  from  layered  bodies  of  revolution.  Extensive 
interactive  computer  aided  design  features  have  been  added 
to  EMSCAT  by  Morgan  and  Connolly.  These  computer-aided 
design  features,  driving  the  change  of  name  from  EMSCAT  to 
EMCAD,  allow  structures  composed  of  multiple  dielectric 
materials  to  be  graphically  constructed  or  modified  by  the 
user,  and  then  evaluated  by  EMCAD  for  scattering.  The 
remainder  of  this  thesis  is  concerned  with  the  utilization 
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and  application  of  EMCAD,  as  well  as  with  modifications  and 
improvements  of  this  transformation. 


III.  EMCAD  INPUT  DATA  CONSTRUCTION 


As  has  been  previously  explained,  EMCAD  computes 
electromagnetic  field  scattering  from  penetrable  bodies  of 
revolution  using  the  coupled  azimuthal  potential  formulation 
in  conjunction  with  a  variational  finite-element  technique 
and  a  tri-regional  unimoment  method.  The  numerical  solution 
is  performed  using  a  two-sweep  Riccati  transform  with 
circumferential  marching.  In  order  to  provide  relatively 
"fast",  consistently  accurate  answers,  yet  be  flexible 
enough  to  allow  the  user  to  vary  any  of  a  number  of  input 
parameters  as  conditions  warrant,  EMCAD  requires  sufficient 
input  data  to  adequately  depict  the  penetrable  body  of 
revolution,  as  well  as  the  desired  incident  and  scattered 
field  orientations  of  the  user.  This  accuracy  and 
flexibility  is  accomplished  through  the  utilization  of 
keyboard  input  and/or  mouse  input  to  form  both  material  and 
structure  input  data  files.  Mouse  input  deals  with  the 
graphical  interface,  and  as  such  will  be  discussed  in 
Chapter  IV. 

A.  KEYBOARD  INPUT 

Information  obtained  by  keyboard  input  is  fairly 
diversified  in  nature.  This  data  ranges  from  input  and 
output  file  names,  to  desired  mesh  parameters,  from 
selection  of  phi  angles,  to  selection  of  incidence  angles. 
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To  facilitate  ease  in  understanding  and  running  the  program, 
EMCAD  allows  a  user  to  input  the  necessary  information 
according  to  the  level  of  expertise  needed.  For  the 
"Novice”  level,  the  program  gives  brief  explanations  of 
input  values,  describes  the  required  format,  and  gives 
typical  values  where  applicable.  For  the  other  option,  the 
"Expert"  level  assumes  the  user  is  familiar  with  the  input 
parameters  and  formats,  and  simply  prompts  for  the  required 
inputs.  A  more  detailed  explanation  of  each  facet  of  the 
keyboard  input  is  contained  in  Appendix  A,  EMCAD  USER'S 
MANUAL. 

B.  DATA  FILE  INPUT 

Once  the  information  obtained  by  keyboard  input  is 
established,  the  remaining  data  becomes  that  which  describes 
the  penetrable  body  of  revolution  being  considered.  This 
descriptive  data  can  be  separated  into  two  distinct  areas  of 
information.  The  first  area  of  information  is  the  data 
which  conveys  the  material  composition  of  each  layer  of  the 
structure.  Call  this  the  "material  file".  The  second  area 
of  information  is  the  data  which  defines  the  structural 
shape  of  each  layer  of  the  penetrable  body,  referred  to  as 
the  "structur  file".  As  shown  in  Figure  3-1,  there  are 
three  methods  of  data  construction  which  can  be  used  to 
generate  the  "Material"  file  and  the  "Structur"  file.  The 
quotation  marks  indicate  that  the  data  files  can  be  named  at 
the  user's  discretion.  Method  2  and  Method  3  utilize 
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SCATTERING  CALCULATION 
INPUT-OUTPUT  RELATIONSHIP 


computer  aided  design  features,  and  will  therefore  be 
discussed  in  Chapter  IV,  Interactive  Graphics  Features. 

Data  construction  by  way  of  Method  1  is  depicted  in 
Figure  3-2  and  is  composed  of  special  case  programs  which 
generate  data  files  for  specific  geometrically  shaped 
objects.  For  example,  program  SPHERE  allows  the  user  to 
define  a  layered  penetrable  sphere,  offset  from  center  if 
desired,  with  as  many  penetrable  layers  as  can  be  handled  by 
EMCAD  at  this  time.  The  user  enters  material  dielectric 
parameters,  layer  radii,  and  required  resolution  as 
prompted,  and  the  program  writes  the  necessary  "Material" 
and  "Structur"  data  into  files  in  the  required  EMCAD 
format.  CONE  and  CYLINDER  were  written  specifically  for 
EMCAD.  As  with  SPHERE,  they  allow  the  user  to  establish  the 
body  of  revolution  as  a  multilayered  structure.  Source  code 
for  SPHERE,  CONE,  and  CYLINDER  can  be  found  in  Appendices  J, 
K,  and  L  respectively.  Whenever  a  computer  program 
generates  data  and  data  files,  the  capability  exists  for  the 
data  to  be  generated  manually,  or  by  any  other  means 
available,  and  then  entered  into  the  appropriate  data  files 
manually  through  the  use  of  the  keyboard.  Keyboard  entry  is 
shown  as  a  viable  method  of  data  construction  in  Method  One, 
Two  and  Three. 
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METHOD  NUMBER  ONE 


Spheres,  cones,  and  cylinders  are  simple  geometric 
shapes,  and  are  extremely  useful  for  elementary 
applications;  however,  nature  is  complex,  and  not  always 
made  up  of  simple  geometry  and  elementary  shapes.  To  obtain 
realistic  results,  realistic  penetrable  bodies  are  needed. 
It  would  indeed  be  tedious  to  attempt  to  write  a  structure 
defining  computer  program  for  each  new  shape  encountered. 
As  such,  the  optimum  solution  indicates  the  need  for  the 
implementation  of  computer  aided  graphics. 
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IV.  INTERACTIVE  GRAPHICS  FEATURES 


Interactive  graphics  simply  refers  to  the  ability  to 
take  data  at  hand,  and  depict  that  data  in  some  visual  form, 
where  a  graph  is  the  most  commonly  utilized  depiction.  A 
special  subset  of  interactive  graphics,  known  as  computer 
aided  design,  also  allows  the  user  to  draw  a  picture  on  the 
screen,  and  convert  that  picture  to  data. 

EMCAD  was  written  in  the  programming  language  of  Ryan- 
McFarland  FORTRAN  (Appendix  M) .  Microcompatibles  GRAFMATIC 
(Appendix  M)  ,  and  later  Jewell  Technologies  GRAFPLUS 
(Appendix  M)  were  chosen  as  the  commercially  available 
screen  graphics  and  hard  copy  support  for  EMCAD  output 
graphics.  Because  the  output  format  that  was  used  readily 
allowed  conversion  to  an  acceptable  format,  Kest  Coast 
Consultants  CURVE  DIGITIZER  (Appendix  M)  was  chosen  as  the 
computer  aided  design  package. 

A.  COMPUTER  AIDED  DESIGN 

There  are  three  means  of  constructing  the  data  for 
penetrable  bodies  of  revolution,  as  previously  shown  in 
Figure  3-1.  Method  1  has  already  been  discussed,  and  is  the 
predecessor  to  Method  2.  In  Method  2,  the  "Material"  data 
can  be  generated  by  one  of  the  three  programs  written 
specifically  for  EMCAD,  or  by  direct  keyboard  input  into  a 
data  file  if  desired.  The  change  from  Method  1  occurs  in 
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the  construction  of  the  structure  defining  data.  As  shown 
in  Figure  4-1,  the  "Curveout"  data  is  obtained  by  using 
CURVE  DIGITIZER.  In  simple  terms,  this  means  that  the  user 
draws  the  picture  on  tne  screen  using  a  mouse  with  the 
computer  aided  design  package.  Once  the  user  is  satisfied 
with  the  single,  or  multilayered  drawing  on  the  screen,  the 
data  is  saved  in  the  format  established  by  the  computer 
aided  design  package.  An  example  of  the  two-dimensional 
cross  section  of  a  three-dimensional  body  of  revolution 
which  was  designed  using  CURVE  DIGITIZER  is  shown  in 
Figure  4-2. 

The  data  file  that  is  saved  is  a  valid  description  of 
the  structure,  but  is  not  in  a  form  which  can  be  used  by 
EMCAD.  This  format  conversion  is  accomplished  by  a  program 
written  specifically  for  this  purpose,  EMCADIN.  The 
computer  aided  design  output  data  is  written  in  X,Y  format, 
with  special  codes  interspersed  throughout.  EMCAD  requires 
radius  format  and  a  normalizing  factor  of  2n  must  be  taken 
into  account.  The  output  data  is  not  necessarily  evenly 
spaced,  whereas  EMCAD  requires  evenly  spaced  input  data.  It 
is  conceivable  that  the  data  is  of  insufficient  resolution 
for  use  by  EMCAD.  Primitive  data  should  be  of  high 
resolution  (i.e.,  361  points).  The  program  EMCADIN 
rectifies  all  of  these  problems,  as  well  as  computes  the 
minimum  and  maximum  radius  values  required.  The  output  of 
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FIG.  4-1.  SCATTERING  CALCULATION 
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SECTION  OF  BODY  OF  REVOLUTION  DESIGNED 
USING  CURVE  DIGITIZER 


EMCADIN  is  the  data  file  "Structur".  A  more  detailed 
description  can  be  found  in  Appendix  D,  EMCADIN  USER'S 
MANUAL. 

Method  3,  as  shown  in  Figure  4-3,  is  actually  a  union 
of  Method  1  and  Method  2.  "Material”  data  is  generated  as 
described  in  both  previous  methods.  Structural  data  is 
generated  as  in  Method  1,  but  is  then  imported  to  the 
computer  aided  design  package  for  modifications.  The  data 
to  be  modified  must  be  changed  by  CRVDIGIN,  written 
especially  for  this  purpose,  into  a  format  which  can  be 
accepted  as  input  by  CURVE  DIGITIZER.  CRVDIGIN  is  explained 
in  greater  detail  in  Appendix  F,  CRVDIGIN  USER'S  MANUAL. 
Once  the  data  has  been  read  into  the  computer  aided  design 
package,  the  remainder  of  the  development  is  essentially 
Method  2 . 

B.  INTERACTIVE  OUTPUT  GRAPHICS 

Having  constructed  the  "Material"  and  the  "Structur" 
input  data  files  using  one  of  the  three  methods  discussed, 
it  is  now  possible  to  conduct  the  scattering  calculations. 
Figure  4-4  shows  that  the  execution  of  EMCAD  yields  output 
data  and  output  graphics. 

As  also  shown  in  Figure  4-4,  EMCAD  generates  four  data 
files  which  are  used  as  input  to  the  graphics  program 
MCADGRAF.  It  is  significant  to  note  that  the  extensions  of 
these  filenames  have  special  meaning.  TEP  refers  to  TE 
incidence  (Figure  2-3),  F-Phi;  TET  refers  to  TE  incidence, 
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DATA 


FIG.  4-3.  SCATTERING  CALCULATION 
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FIG.  4-4.  SCATTERING  CALCULATION 

OUTPUT  AND  GRAPHICS 


F-Theta;  TMP  refers  to  TM  incidence  (Figure  2-3),  F-Phi;  TMT 
refers  to  TM  incidence,  F-Theta. 

MCADGRAF  is  a  menu  driven  output  graphics  program. 
Appendix  H,  MCADGRAF  USER'S  MANUAL,  gives  examples  of  all 
the  different  input  menus.  The  menu  structure  gives  the 
user  some  general  information  concerning  the  use  of  the 
program,  as  well  as  user  choice  of  aspect  ratio,  data  file 
type  selection,  data  file  name  input,  incidence  angle 
selection,  and  phi  angle  selection.  Once  these  selections 
have  been  made,  MCADGRAF  produces  four  output  graph  types, 
giving  the  user  the  capability  to  make  a  printout  of  any  or 
all  of  these.  The  first  graph  type  is  magnitude  versus 
theta  in  a  rectangular  coordinate  system,  as  shown  in  Figure 
4-5.  The  second  graph  type  is  phase  versus  theta  in  a 
rectangular  coordinate  system,  as  shown  in  Figure  4-6.  The 
third  graph  type  is  magnitude  versus  theta  in  a  polar 
coordinate  system,  shown  in  Figure  4-7.  The  final  type  is 
simply  an  axis  rotated  version  of  the  third  type,  for  a 
different  perspective,  as  shown  in  Figure  4-8. 
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FIG.  4-5.  MCADGRAF  GRAPHICAL  OUTPUT 
MAGNITUDE  VERSUS  THETA 


Hardcopy  — )  P  or  p 


FIG.  4-6.  MCADGRAF  GRAPHICAL  OUTPUT 


PHASE  VERSUS  THETA 


Offset  Sphere  Validation  Run 
Scattered  Field  Coup  lex  Amplitudes 
IE  Incidence  =  90.  Degrees  Phi  =  0.  Degrees 

Polar  Plot i  F-Phi  Magnitude  Versus  Theta 


FIG.  4-7.  MCAOGRAF  GRAPHICAL  OUTPUT 
MAGNITUDE  VERSUS  THETA.  POLAR  PLOT 
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Composite  Body 

Scattered  Field  Conplex  Amplitudes 
TM  Incidence  =  135.  Degrees  Phi  =  180.  Degrees 
Polar  Plot,  F-Theta  Magnitude  Versus  Theta 
Axis  Rotated  Minus  Ninety  Degrees 
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FIG.  4-8.  MCADGRAF  GRAPHICAL  OUTPUT 
MAGNITUDE  VERSUS  THETA,  POLAR  PLOT 
AXIS  ROTATED  MINUS  NINETY  DEGREES 
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V.  VALIDATION 


In  order  to  transform  the  existing  program  EMSCAT  to 
the  modified  version  entitled  EMCAD,  it  was  imperative  to 
ensure  that  both  the  accuracy  and  the  reliability  of  EMCAD 
were  at  least  as  good  as  that  of  EMSCAT.  It  was  also 
essential  to  establish  flawless  compatibility  between  EMCAD 
and  any  program  written  to  support  EMCAD. 

The  initial  step  used  to  validate  the  operation  of 
EMCAD  consisted  of  a  comparison  between  the  results  obtained 
from  EMCAD  and  the  exact  theoretical  values  of  a  Mie  series 
far-field  complex  amplitude  E-plane  pattern  for  an  offset 
lossless  (er  *  2.6)  penetrable  dielectric  sphere.  Run 
parameters  chosen  included  broadside  incidence  (incidence 
angle  =  90  degrees),  sphere  radius  equal  to  0.8  (as  a  factor 
of  wavelength),  and  z-axis  offset  =  0.16  (as  a  factor  of 
wavelength).  The  series  truncations  at  M=7,  NMI=8  and 
Nmii»15  were  chosen  for  the  variable  mesh  scheme.  These 
values  were  the  same  as  those  used  by  Morgan  [Ref.l]  in  his 
initial  validation  of  EMSCAT.  The  results  of  this 
comparison  are  shown  in  Figure  5-1  for  a  phi  angle  of  0 
degrees,  and  in  Figure  5-2  for  a  phi  angle  of  90  degrees. 

Having  demonstrated  the  close  resemblance  between  the 
results  obtained  from  EMCAD,  and  those  predicted  by  theory, 
the  next  step  used  to  validate  EMCAD  required  a  comparison 
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OFFSET  SPHERE  SCATTERING  (#2) 


between  the  results  obtained  from  EMCAD  and  actual  measured 


data.  Morgan  [Ref.  1]  validated  the  original  version  of 
EMSCAT  with  extensive  experimental  verification.  Having  the 
results  of  this  experimental  verification  readily  available 
presented  the  unique  opportunity  to  validate  the  operation 
of  EMCAD  with  the  same  experimental  data  used  to  validate 
EMSCAT.  For  each  axisymmetric  body  constructed,  parameters 
chosen  included  TE  and  TM  incidence  with  incidence  angles  of 
0,  45,  90,  135  and  180  degrees.  The  series  truncations  at 
M=9,  NMI=15  and  N^II=2  0  were  chosen  for  the  variable  mesh 
scheme.  All  experimental  bodies  were  constructed  of 
plexiglass  (polymethylmethacrylate)  with  material  parameters 
defined  by  er  =  2.59  -  j 0.017.  [Refs.  1  and  2] 

The  decision  was  made  to  demonstrate  the  versatility 
and  reliability  of  the  program  by  utilizing  a  number  of 
significant  shapes  which  could  be  readily  fabricated  for 
experimental  measurements.  Since  a  sphere  had  already  been 
used,  and  shown  to  work,  a  bicone  was  chosen  as  an 
appropriate  figure.  A  representation  of  the  solid  bicone 
utilized  is  depicted  in  Figure  5-3.  The  comparison  for  the 
solid  bicone  validation  is  shown  in  Figure  5-4.  The  hollow 
bicone  of  Figure  5-5  demonstrates  increasing  complexity. 
The  sphere  shown  inside  the  bicone  is  a  free  space  cavity  in 
the  object,  although  it  could  be  made  of  any  material.  The 
comparison  for  the  hollow  bicone  is  presented  in  Figure  5-6. 
The  next  figure  chosen  was  that  of  a  hollow  cylinder.  The 
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HOLLOW  BICONE  SCATTERING 


internal  spherical  cavity  had  the  same  dimension  as  that 
found  in  the  hollow  bicone.  The  height  of  the  cylinder  was 
1.9  (in  terms  of  wavelength)  and  the  radius  of  the  base  of 
the  cylinder  was  .945  (in  terms  of  wavelength).  The 
validation  comparison  for  a  phi  angle  of  0  degrees  is 
depicted  in  Figure  5-7  and  the  validation  comparison  for  a 
phi  angle  of  90  degrees  is  presented  in  Figure  5-8. 

To  this  point,  all  validation  runs  have  correlated 
extremely  well  with  the  experimental  data,  and  all 
experimental  figures  have  consisted  of  combinations  of 
simple  geometric  shapes  such  as  spheres,  cones,  and 
cylinders.  These  are  extremely  useful  for  elementary 
applications,  however,  nature  is  complex,  and  not  always 
made  up  of  nice  geometry  or  elementary  shapes.  To  prove 
that  EMCAD  could  handle  more  realistic  bodies,  the  composite 
structure  of  Figure  5-9  was  constructed.  Since  a  program 
was  not  written  to  generate  the  "structur"  data  required  for 
EMCAD  input,  the  decision  was  made  to  utilize  CURVE 
DIGITIZER  (CAD  package)  to  produce  the  required  information. 
The  representation  of  the  composite  body  of  Figure  5-9, 
required  as  input  for  EMCAD,  is  depicted  in  Figure  5-10. 
Note  that  this  representation  is  composed  of  data  defining  a 
two  dimensional  axisymmetric  body.  The  actual  comparison 
for  this  structure  is  shown  in  Figure  5-11.  As  was  the  case 
before,  the  data  generated  by  EMCAD  correlates  extremely 
well  with  the  experimental  data. 
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SCATTERING  CALCULATION  VALIDATION 
(COORDINATES  IN  TERMS  OF  WAVELENGTH) 


Composite  Body  Validation  Run 

Incidence  Angle  —135  Degrees 
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FIG.  5-11.  COMPOSITE  BODY  SCATTERING 


In  addition  to  demonstrating  the  validity  of  the 
present  version  of  EMCAD,  the  composite  body  of  Figure  5-9 
helped  demonstrate  program  compatibility.  Since  the  output 
data  of  CURVE  DIGITIZER  was  not  in  the  format  required  for 
use  by  EMCAD,  the  computer  program  EMCADIN  was  developed  and 
written  to  transform  the  output  data  of  CURVE  DIGITIZER  into 
a  format  which  could  be  used  by  EMCAD.  Figures  5-12,  5-13 
and  5-14  present  figures  which  were  generated  by  the 
programs  SPHERE,  CONE,  and  CYLINDER,  transformed  by  CRVDIGIN 
into  proper  format,  and  imported  into  the  computer  aided 
design  package  CURVE  DIGITIZER.  The  pictures  shown  are  the 
graphical  representations  of  the  imported  data  as  seen  by 
CURVE  DIGITIZER.  This  generation  and  transformation  not 
only  shows  the  validation  of  the  programs  written  to 
generate  specific  geometric  shapes,  but  also  demonstrates 
the  compatibility  of  these  programs  with  CRVDIGIN,  and  the 
compatibility  of  CRVDIGIN  with  CURVE  DIGITIZER. 
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MULTI-LAYERED.  MULT  I -HE  I  GUT 
MULTI-RADIUS  CYLINDER 
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DATA  OF  PROGRAM  CYLINDER 


VI .  CONCLUSIONS 


This  thesis  has  described  a  very  useful  numerical 
algorithm  for  computing  resonance  region  electromagnetic 
interaction  with  arbitrarily  shaped  and  even  inhomogeneous 
axisymmetric  penetrable  objects.  The  background  material 
associated  with  the  development  of  this  interactive,  finite 
element  scattering  algorithm  was  presented  in  Chapter  11. 
Included  in  this  background  information  was  the  concept  of 
coupled  azimuthal  potentials,  the  finite  element  algorithm 
and  the  unimoment  method.  Chapter  III  then  presented 
procedures  for  constructing  input  data  for  the  scattering 
program.  Computer  programs  written  specifically  to  generate 
a  given  geometry,  and  to  allow  data  file  input  and  keyboard 
input  were  examined  in  detail.  The  first  half  of  Chapter  IV 
discussed  the  extension  of  the  methods  used  to  construct 
program  input  data  to  those  methods  utilizing  graphics.  Of 
special  interest  was  the  incorporation  of  computer  aided 
design  features.  The  second  half  of  the  chapter  presented 
user  interactive  output  graphics  features.  Validations  of 
the  modifications  and  improvements  to  EMCAD  were  presented 
in  Chapter  7.  Extensive  graphical  comparisons  were  made 
between  the  output  of  the  modified  version  of  EMSCAT,  and 
both  theoretical  calculations  and  actual  measured 
experimental  data.  The  impressive  accuracy  of  the 
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comparisons  validated  the  changes  and  modifications  to 
EMSCAT  used  in  developing  EMCAD.  Inter-program 
compatibility  was  also  extensively  demonstrated.  The 
Appendices  provide  detailed  users'  manuals  for  four  of  the 
major  computer  programs  developed. 

This  was  only  the  first  stage  of  a  multistage  project. 
Future  developments  of  the  axisymmetric  code  will  proceed 
with  major  programming  modifications  to  allow  metallic 
surfaces  to  be  embedded  in  the  inhomogeneous  penetrable 
material  of  the  scattering  object.  Also,  the  near  fields  of 
the  object  will  be  made  available  for  display.  This  is 
useful  for  identifying  "hot  spots"  and  field  caustics. 

A  second  effort  has  already  begun  to  develop  a  fully 
three-dimensional  scattering  object  capability.  The  basis 
of  this  development  was  the  discovery  of  a  generalized 
coupled  azimuthal  potential  formulation  which  does  not 
require  an  axisymmetric  material  configuration  [Ref.  8]. 
Such  a  computer  algorithm  will  ultimately  employ  three- 
dimensional  graphics  in  a  powerful  computer  aided  design 
package  for  user  design  of  scattering  objects  as  well  as  the 
output  of  scattering  patterns  and  near  field  quantities. 


APPENDIX  A 


EMCAD  USER'S  MANUAL 

Welcome  to  EMCAD.  This  program  computes 
electromagnetic  field  scattering  from  layered  bodies  of 
revolution  using  the  coupled  azimuthal  potential  formulation 
in  conjunction  with  a  variational  finite  element  technique 
and  a  tri-regional  unimoment  method.  The  numerical  solution 
is  performed  using  a  two-sweep  Riccati  transform  with 
circumferential  marching. 

A .  BACKGROUND 

EMCAD  is  based  upon  an  earlier  code  named  EMSCAT  which 
was  written  during  the  period  1974-1976  by  Dr.  M.  A.  Morgan 
while  at  the  University  of  California  at  Berkeley. 
Extensive  interactive  computer  aided  design  (CAD)  type 
features  have  been  added  in  this  effort.  These  CAD  features 
allow  structures  composed  of  multiple  dielectric  materials 
to  be  graphically  constructed  or  modified  by  the  user,  and 
then  evaluated  by  EMCAD  for  scattering.  EMCAD  has  been  run 
in  both  Ryan-McFarland  FORTRAN  (Appendix  M)  and  Microsoft 
FORTRAN  (Appendix  M) ,  with  the  current  version  using 
Microsoft.  Microcompatibles  GRAFMATIC  (Appendix  M)  and 
later  Jewell  Technologies  GRAFPLUS  (Appendix  M)  were  chosen 
as  the  commercially  available  screen  graphics  and  hard  copy 
support  for  EMCAD  output  graphics.  Because  its  output 
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format  allowed  easy  conversion  to  an  acceptable  format,  West 
Coast  Consultants  CURVE  DIGITIZER  (Appendix  M)  was  chosen  as 
the  computer  aided  design  package.  Any  questions  concerning 
colors,  formats,  changes  to  the  program,  etc...  should  be 
referred  to  the  respective  user's  manuals,  Dr.  M.  A.  Morgan 
or  LT.  E.  M.  Connolly  (Appendix  M) . 

B.  COMPUTER  SETUP 

EMCAD  was  written  for  the  following  computer 
configuration: 

*  IBM  AT  or  compatible,  including  80386  based 
compatibles 

*  hard  disk 

*  EGA  graphics  card 

*  EGA  monitor 

*  math  coprocessor  -  not  required,  but 
highly  recommended 

*  memory  -  64 OK  RAM  minimum 

EMCAD  can  be  placed  in  any  directory,  at  the  user's 
discretion,  however,  if  a  RAM-drive  is  to  be  used,  is  highly 
recommended  that  during  each  run  the  user  allocate  a  minimum 
of  4M  of  RAM  for  EMCAD  and  EMCAD  input/output  files.  EMCAD 
input  is  taken  from  a  variety  of  sources.  The  data  files 
from  these  sources  are  required  to  be  in  the  same 
directory/ sub-directory  as  EMCAD. 
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C.  PROGRAM  INPUT 


EMCAD  requires  both  keyboard  input  and  data  file  input. 
As  shown  in  Figure  A-l,  two  input  data  files,  "MATERIAL" . DAT 
and  "STRUCTURE" . DAT,  are  required  as  input  to  EMCAD.  The 
quotation  marks  indicate  that  the  data  files  can  be  named  at 
the  user's  discretion.  These  two  data  files  contain 
information  which  describes  the  penetrable  body  of 
revolution  being  considered.  The  "material  file"  contains 
the  data  which  conveys  the  material  composition  of  each 
layer  of  the  structure.  The  "structure  file"  contains  the 
data  which  defines  the  structural  shape  of  each  layer  of  the 
penetrable  body.  The  three  methods  which  are  used  to 
generate  these  data  files  are  discussed  in  much  greater 
detail  in  both  Chapter  II  and  Chapter  III. 

D.  PROGRAM  EXECUTION 

EMCAD  is  a  menu  driven  program.  Once  the  correct 
directory  has  been  selected,  EMCAD  can  be  invoked  by  typing 

EMCAD 


at  the  DOS  prompt.  This  is  shown  in  the  top  line  of  Figure 
A- 2 .  The  program  will  load  itself,  and  the  user  will  see 
the  next  two  statements  shown  in  Figure  A-2  requesting  the 
names  of  the  "material  file"  and  the  "structure  file". 
Sample  keyboard  inputs  are  shown  for  these  two  requests.  It 
must  be  pointed  out  that  the  extensions  are  included  (.dat) 
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SCATTERING  CALCULATION 
INPUT-OUTPUT  RELATIONSHIP 


C:\IHESISNNSFOBDEHCAD 


ENTER  HATER  I AL  PARAMETER  FILE  (D: FILENAME. EXTENSION) 
ndatl.dat 


ENTER  STRUCTURE  DATA  FILE  (D : FILENAME . EXTENSION ) : 
sdatl.dat 

Offset  Sphere  R=8.8  D=8.16 
Er=(2.6»0. )  Ur=(l.,8.) 


PLEASE  PRESS  ANY  HEY  TO  CONTINUE. 


FIG.  A-2.  INPUT  DATA  FILES 
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as  required  by  EMCAD.  The  two  lines  beginning  with  "Offset 
Sphere..."  refer  to  two  label  lines  which  were  previously 
inserted  into  the  "material  file"  to  allow  the  user  to 
differentiate  this  set  of  data  from  every  other  set  of  data. 
By  choosing  labels  carefully,  the  user  can  discern  from  the 
very  first  screen  whether  the  proper  set  of  data  is  being 
loaded.  This  becomes  very  important  because  as  figures 
increase  in  complexity,  run  time  for  the  program  increases. 
It  is  not  uncommon  to  have  runs  which  take  several  hours  on 
an  AT  class  computer.  It  would  be  very  disturbing  to  run 
the  program  for  seven  or  eight  hours  only  to  find  out  that 
the  incorrect  set  of  data  had  been  used.  Upon  completion  of 
this  screen,  the  user  depresses  ENTER  to  see  the  following 
screens: 

*  Welcome  To  EMCAD  (Figure  A-3) 

This  screen  presents  some  brief  background  information. 

*  EMCAD  Level  of  Expertise  (Figure  A-4) 

To  facilitate  ease  in  understanding  and  running  the  program, 
EMCAD  allows  a  user  to  input  the  necessary  information 
according  to  the  level  of  expertise  needed.  For  the 
"Novice"  level,  the  program  gives  brief  explanations  of 
input  values,  describes  the  required  format,  and  gives 
typical  values  where  applicable.  The  "Expert"  level 
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******  WELCOME  TO  ENCAD  ****** 

THIS  PROGRAM  COMPUTES  EH  FIELD  SCATTERING  FROM  PENETRABLE 
BODIES  OF  REVOLUTION  USING  THE  COUPLED  AZIMUTHAL  POTENTIAL 
(CAP)  FORHULATIOM  IN  CONJUNCTION  WITH  A  VARIATIONAL  FINITE- 
ELEMENT  TECHNIQUE  AND  A  TRI-REG IONAL  UNIMOMENT  METHOD .  THE 
NUMERICAL  SOLUTION  IS  PERFORMED  USING  A  TWO-SWEEP  RICCATI 
TRANSFORM  WITH  CIRCUMFERENTIAL  MARCHING. 

UPDATED  VERSION  OF  ORIGINAL  U.C.  BERKELEY  CODE. 
MODIFICATIONS  BY  H.A.  HORGAN  MAR  1987  -  FEB  1988 

E.H.  CONNOLLY  JUL  1987  -  APR  1988  . 


PLEASE  PRESS  ANY  KEY  TO  CONTINUE. 


FIG.  A-3.  WELCOME  TO  EMCAD 


*******  ENCAD  INPUT  »»****» 


ENCAD  ALLOWS  A  USEE  TO  INPUT  THE  NECESSARY  INFORNATION 
ACCORDING  TO  THE  LEVEL  OF  HIS  EXPERTISE. 

1  NOVICE  LEVE  -  GIVES  BRIEF  EXPLANATIONS/DESCRIPTIONS 

OF  INPUT  VALUES.  DESCRIBES  THE  REQUIRED 
FORMAT  FOR  THE  INPUT  VALUE.  GIVES 
TYPICAL  VALUES  WHERE  APPLICABLE. 

2  EXPERT  LEVEL  -  ASSUMES  THE  USER  IS  FAMILIAR  WITH  THE 

INPUT  PARAMETERS  AND  FORMATS.  SINPLY 
PROMPTS  FOR  REQUIRED  INPUTS. 

3  EXII  ENCAD 

PLEASE  SELECT  THE  LEVEL  OF  EXPERTISE  ENTERING  1,  2,  OR  3. 


FIG.  A- 4.  EMCAO  LEVEL  OF 
USER  EXPERTISE 
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assumes  the  user  is  familiar  with  the  input  parameters  and 
formats,  and  simply  prompts  for  the  required  inputs.  The 
remaining  menus  will  use  the  "Novice”  level  to  cover  all 
ranges  of  user. 

*  Input  Name  for  Output  Data  File 

(Figure  A-5) 

This  menu  prompts  the  user  for  the  filename,  WITHOUT 
EXTENSION,  that  will  be  used  to  hold  the  output  data  of 
EMCAD.  The  program  will  append  the  extension  .OUT  to  the 
user's  filename. 

*  Input  Name  for  Output  Graphics  Data  Files 

(Figure  A- 6) 

EMCAD  produces  four  output  data  files  for  input  to  the 
graphing  routines  of  MCADGRAF.  All  four  data  files  have  the 
same  name,  but  each  has  a  distinct  extension  related  to  the 
type  of  data  the  file  contains.  This  menu  prompts  the  user 
for  the  name  of  the  output  graphics  data  files,  WITHOUT 
EXTENSION.  This  filename  may  be  the  same  as  that  used  for 
the  output  data. 

*  Graphics  Caption  ,  (Figure  A-7) 

This  menu  affords  the  user  the  opportunity  to  input  a 
graphics  caption  to  help  identify  this  set  of  data  from  all 
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THE  OUTPUT  DATA  FILE  IS  THE  FILE  CONTAINING  THE  OUTPUT 
RESULTS  OF  ALL  NUMERICAL  CALCULATIONS  CONDUCTED  BY  ENCAD. 
THE  FORMAT  FOR  THIS  INPUT  IS  FILENAME  ONLY.  NO  EXTENSION  IS 
REQUIRED  OR  DESIRED.  ENCAD  AUIOHAIICALLY  APPENDS  AN 
EXTENSION  OF  .OUT  10  YOUR  FILENAME. 

PLEASE  ENTER  THE  FILENAME  OF  THE  OUTPUT  DATA  FILE. 


FIG.  A-5.  INPUT  NAME  FOR 


OUTPUT  DATA  FILE 
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THE  OUTPUT  GRAPHICS -DAIA  FILE  IS  THE  FILE  CONTAINING  THE 
OUTPUT  DAIA  FRON  ENCAD  TO  BE  USED  AS  INPUT  TO  GRAPHING 
ROUTINES.  THE  FORMAT  FOR  THIS  INPUT  IS  FILENAHE  ONLY.  NO 
EXTENSION  IS  REQUIRED  OR  DESIRED.  ENCAD  AUTOMATICALLY 
APPENDS  AN  EXTENSION  OF  .IHT,  .TNP,  .TET,  AND  .TEP  TO  YOUR 
FILENAME  AS  II  PRODUCES  FOUR  OUTPUT  FILES  FOR  GRAPHICS. 

.TNT  - >  TN  INCIDENCE,  F-THETA 

.IMP  - )  TN  INCIDENCE,  F-PHI 

.TET - >  TE  INCIDENCE,  F-THETA 

.TET - >  TE  INCIDENCE,  F-PHI 

PLEASE  ENTER  THE  FILE  NAME  OF  THE  OUTPUT  DATA  FILE. 


FIG.  A-6.  INPUT  NAME  FOR 
OUTPUT  GRAPHICS  DATA  FILES 
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THE  GRAPHICS  CAPTIOH  IS  A  PERSONALIZED  CAPTION  ALLOWING  THE 
USER  TO  IDENTIFY  THIS  SET  OF  GRAPHS  FRON  ALL  OTHER  SETS. 

THE  NAXIHUfl  LENGTH  OF  THIS  CAPTION  IS  64  CHARACTERS, 

NOTE:  WHEN  USED  WITH  THE  GRAPHICS  PACKAGE,  THE  PROGRAN  IS 
ABLE  TO  DIFFERENTIATE  BETWEEN  UPPER  CASE  AND  LOWER  CASE 
CHARACTERS. 

PLEASE  ENTER  ANV  GRAPHICS  CAPTION  YOU  DESIRE 


FIG.  A-7.  GRAPHICS  CAPTION 
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other  sets.  The  graphics  caption  appears  on  the  first  line 
of  the  header  of  each  graph  produced  by  MCADGRAF. 

*  Mesh  Density  (Figure  A-8) 

As  discussed  in  Chapter  II,  EMCAD  utilizes  a  semi-annular 
conformal  finite-element  mesh  in  order  to  conduct  scattering 
calculations.  This  menu  prompts  the  user  for  the  minimum 
and  maximum  mesh  density  parameters.  At  azimuthal  modes  M  = 
0  and  1,  is  used  as  the  density  parameter.  The  density 

decreases  linearly  to  at  the  final  mode  step.  Notice 

that  these  values  are  REAL. 

*  Incidence  Angles  (Figure  A-9) 

EMCAD  allows  up  to  five  incident  fields  to  impinge  on  the 
object  of  interest.  This  menu  asks  for  the  number  of 
incidence  angles,  and  then  the  value  for  each  in  degrees. 

*  Number  of  Scattering  Field  Theta  Points 

(Figure  A-10) 

EMCAD  conducts  scattering  calculations  by  starting  at  6  =  0 , 
and  progressing  clockwise  to  9  =  180  in  equal  increments. 
This  menu  allows  the  user  to  determine  the  size  of  the 
increments,  and  indirectly  the  accuracy. 
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DMIM  AND  DflAX  ARE  PARANETERS  OF  HESH  DENSIIY  IN  TERNS  OF 
ELEMENTS/ INTERIOR  LAMBDA.  INPUT  VALUES  ARE  EXPECTED  TO  BE 
REAL,  I.E.  THE  DECIMAL  POINT  MUST  BE  INCLUDED. 

TYPICAL  VALUES  ARE  DMIN  =  10.  AND  DMAX  =  15. 

PLEASE  INPUT  DMIN 

10. 


PLEASE  INPUT  DMAX 
15. 


FIG.  A-8.  MESH  DENSITY 
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THE  HUMBER  OF  INCIDENT  FIELD  ANGLES  IS  THE  TOTAL  NUMBER  OF 
INCIDENT  FIELDS  THAT  IMPINGE  ON  THE  OBJECT  OF  INTEREST. 
THIS  PROGRAM  ALLOWS  A  MAXIMUM  OF  FIUE  INCIDENT  FIELD 
ANGLES.  WHEN  ENTERING  TOUR  ANSWER  PLEASE  DO  NOT  INCLUDE  A 
DECIMAL  BECAUSE  THE  INPUT  MUST  BE  IN  INTEGER  FORMAT,  I.E. 


NA  =  3  ,  NA  =  1 


PLEASE  INPUT  THE  NUMBER  OF  INCIDENT  FIELD  ANGLES. 
2 


ENTER  INC  FLD  ANGLE  (DEG)  FOR  «  1 

8. 


ENTER  INC  FLD  ANGLE  (BEG)  FOR  »  2 
90. 


FIG.  A- 9.  INCIDENCE  ANGLES 
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THE  NUMBER  OF  SCATTERING  FIELD  THETA  POINTS  DETERMINES  THE 
SPACING  BETWEEN  THETA  POINTS  DURING  ENSCAT  ITERATIONS  AND 
CALCULATIONS. 

DELTA  THETA  =  188  /  (NUMBER  THETA  POINTS  -  1)  SO. .. 

NUMBER  THETA  POINTS  =  3?  - >  DELTA  THETA  =  5  DEGREES 

NUMBER  THETA  POINTS  =  19  - >  DELTA  THETA  =  10  DEGREES 

WHEN  ENTERING  ¥OUR  ANSWER,  PLEASE  DO  NOT  INCLUDE  A  DECIMAL 
BECAUSE  THE  INPUT  MUST  BE  IN  INTEGER  FORMAT.  I.E.  NT  =  19 

PLEASE  INPUT  THE  NUMBER  OF  SCATTERING  FIELD  THETA  POINTS 


FIG.  A-10.  NUMBER  DF 
SCATTERING  FIELD  THETA  POINTS 
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*  Phi  Angles 


(Figure  A-ll) 


EMCAD  allows  up  to  eight  phi  angles  which  define  the 
location  from  which  the  user  views  the  scattering  for  all  of 
the  incident  angles.  This  menu  asks  for  the  number  of  phi 
angles,  and  then  the  value  for  each  in  degrees. 

*  Expansion  Modes  (Figure  A-12) 

EMCAD  partitions  the  incident  and  scattered  fields  into 
azimuthal  modes,  indexed  by  "M" .  The  solutions  for  these 
modes  are  accumulated  during  the  "M-loop".  If  0  degrees  and 
180  degrees  are  the  only  incidence  angles  of  interest,  then 
the  M  =  1  mode  is  the  only  mode  that  is  needed.  The  values 
of  N1  and  N2  represent  the  number  of  spherical  harmonic 
modes  in  the  spherical  core  and  the  external  region. 

E .  PROGRAM  OUTPUT 

As  previously  mentioned,  EMCAD  produces  text  output  data 
and  output  data  for  graphics.  The  graphics  data  is 
discussed  in  Appendix  H,  MCADGRAF  User's  Manual,  and  the 
text  output  data  is  demonstrated  in  Appendix  B,  Portion  of 
an  EMCAD  Sample  Output. 
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THE  NUMBER  OF  PHI  ANGLES  IS  THE  TOTAL  NUMBER  OF  PHI  ANGLES. 
THIS  PROGRAM  ALLOWS  A  MAXIMUM  OF  EIGHT  PHI  ANGLES.  WHEN 
ENTERING  MUR  ANSWER  PLEASE  DO  NOT  INCLUDE  A  DECIMAL 
BECAUSE  THE  INPUT  MUSI  BE  IN  INTEGER  FORMAT.  I.E. 


NP  =  3  ,  NP  =  1 


PLEASE  INPUT  IHE  NUMBER  OF  PHI  ANGLES. 
2 


ENTER  PHI  FLD  ANGLE  (DEG)  FOR  I  1 


90. 


ENTER  PHI  FLD  ANGLE  (DEG)  FOR  i  2 
180. 


F I G.  A- 1 1 .  PHI  ANGLES 
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ENTER  flSTOP  (.LE.  13) 


ESTIHAIED  "HININUH"  UALUE  IS:  1 

1 

ENTERING  NO.  OF  INIERNAL  AND  EXTERNAL  EXPANSION  NODES 
ESTIMATED  (KI*HHIN)  "NINIHUH"  INTERNAL  N1  IS:  12 
ENTER  N1  (.GE.  MSTOP) 

12 

ESTIMATED  (MIN)  ’’MINIMUM”  INTERNAL  N2  IS:  21 
ENTER  N2  (.GE.  MSTOP) 

21 

FIG.  A- 1 2.  EXPANSION  MODES 


83 


APPENDIX  B 


PORTION  OP  AN  ENCAD  SAMPLE  OUTPUT 

This  Appendix  demonstrates  a  portion  of  an  EMCAD 
sample  output  for  the  specific  case  chosen  where 
M  indexes  from  0  to  7  inclusive.  Only  the  M  =  1 
case  is  fully  shown,  as  well  as  the  accumulated 
solution  for  all  eight  modes. 

xxxxxxxxxxxxxxxxxxxx  EMCAD  OUTPUT  DATA  xxxxxxxxxxxxxxxxxxxxxxxxxxx 
COMPLEX  Er(n)  AND  Ur(n) 

1  (  2 . 610E+00  -1.000E-06)  (  1.000E+00  -1.000E-06) 

INCIDENT  FIELD  ANGLES 
N  ALPHA(N) 

1  0 .  DEG 

2  90.  DEG 

3  180.  DEG 

SERIES  SOLUTION  PARAMETERS 
MST  ART  =  0 

M3TQP  =  7 

Nl=  3 

N2=  15 
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**************  M 


1  ************** 


MESH  PARAMETERS 
RADII  (WAVELENGTHS) 

RADMIN=  . 64000E+00 

RADMAX=  .96000 E+00 

RADA=  . 55241 E+00 

RADB=  . 1 0476  E+0 1 

INTEGER  CONSTANTS 
NDS-  13 

IMX=  59 

NR1  =  1 

NR2-  1 

MESH  DENSITY  AND  ELEMENT  SIDE  RATIOS 
DENSITY/ LAMB DA -OUT  =  24.233 

DENSITY/LAMBDA-IN  =  15.000 

DSA/DR  =  .725 

DSB/DR  =  1.375 


CCATI 

MARCHING 

1 

.  603E  +  01 

1 

.603E+01 

2 

.  602E+01 

3 

.602E+01 

4 

.  6  00E  +  01 

5 

.  599E+01 

6 

.  597E  +  01 

7 

.  595E+01 

S 

.  592E  +  01 

9 

.  589E  +  01 

10 

.  586E+01 

11 

.  583E+01 

12 

.579E+01 

13 

. 575E+01 

14 

.571E+01 

15 

.566E+01 

16 

.  562E+01 

17 

.  557E+01 

18 

.  552E+01 

19 

.  547E+01 

20 

.  542E+01 

21 

.  537E+01 

22 

. 531 E+  0 1 

23 

.526E+01 

24 

. 520E+01 

25 

. 515E+01 

26 

.  509E  +  01 

27 

.  503E+01 

28 

.  498E+01 

29 

.  492E+01 

30 

. 437  E+D 1 

31 

. 482E+01 

32 

477E+01 

33 

. 471E+01 

34 

.466E+01 

35 

.461E+01 

36 

.457E+01 

37 

.452E+01 

38 

.  448E+01 

39 

.443E+01 

40 

.  439E+01 

41 

.  435E+01 

42 

.  432E+01 

43 

.  428E+01 

44 

.  425E+01 

45 

.  422E+01 

46 

.  419E+01 

47 

.416E+01 

43 

.  414E+01 

49 

.  412E+01 

50 

.  410E+01 

51 

.  408E+01 

52 

.4Q6E+01 

53 

.  405E+01 

54 

.  404E+01 

55 

.  403E  +  01 

56 

.  403E  +  01 

57 

.  402E  +  0 1 

58 

.402E+01 

59 

.492E+01 

INTERIOR  AND  SCATTERED  FIELD  EXPANSION  COEFFICIENTS 


AZIMUTHAL  MODE  NUMBER  M-  1 

INCIDENCE  ANGLE  ALPHA  =  0.  DEGREES 


*******  INTERIOR  COEFFICIENTS  ******* 


STAGGERED 

COEFFICIENTS 

:  TM 

INCIDENCE 

N 

ACN) 

B(  N ) 

I 

.  500E+00 

. 266E+00 

-  .  318E+00 

.  298E+00 

2 

.  209E+00 

- 

. 773E-01 

.  233E+0Q 

.  192Ef00 

3 

I03E+00 

- 

. 177E+00 

.  172E+00 

-  .  60OE-OI 

4 

- . 207E+00 

.582E-01 

.  163E-01 

-.  147E+00 

5 

. 996E-02 

. 169E-01 

-.135E+00 

.  9S1E-02 

6 

. 886  E  —  0 1 

. 390E-01 

109E+00 

. 599 E-0 1 

7 

.614E-01 

. 138E-01 

-  373E-01 

.  437  E-01 

8 

. 197E-01 

— 

. 325E-02 

-  .  625E-02 

. 1 36  E-0 1 

STAGGERED 

COEFFICIENTS 

*  TE 

INCIDENCE 

N 

ACN) 

B(N) 

1 

. 266E+00 

- 

. 500E+00 

.  298E+00 

.  318E+00 

2 

-.778E-01 

- 

. 209E+00 

.  192E+00 

-  .  233E  +  00 

3 

- . 177E+00 

. 103E+00 

-  .  660E-01 

-  .  1  ?  2  E+  0  0 

4 

.582E-01 

. 2O7E+0Q 

-  .  147E+00 

-  163E-01 

5 

.  169E-01 

- 

996E-02 

98 1 E-0  2 

.  1  35  E->-  0  0 

6 

.390E-01 

- 

386E-01 

. 599E-01 

.  1  0  9  £■>-  0  0 

7 

.  158E-01 

- 

.614E-01 

. 4  37  E-0 1 

.  37  3E-01 

8 

~ . 325E-02 

. 1 97  E-0 1 

. 136E-01 

■625E-02 

*******  SCATTERED 

COEFFICIENTS 

******* 

STAGGERED 

COEFFICIENTS 

:  TM 

INCIDENCE 

N 

CCN) 

DC  N ) 

1 

902E-01 

- 

. 144E-01 

. 104E+00 

.  224E  +  00 

2 

- . 115E-01 

- 

. 195E+00 

.  955E-02 

-  .  123E+00 

3 

- . 103E+Q0 

. 157E+00 

- . 994E-01 

-  .  350E-01 

4 

.  138E+00 

. 131E+00 

- . 928  E-0 1 

.  696E-01 

5 

.472E-01 

- 

. 954E-01 

. 6  7  0  E- 0 1 

.  158E+00 

6 

.  168E-01 

- 

.616E-01 

. 5 1 0  E-0 1 

.  659E-01 

7 

.361E-02 

- 

. 1 7 1 E-0  1 

. 159E-01 

.  135E-01 

8 

.  466E-03 

- 

. 306E-02 

. 298E-02 

.  160E-02 

9 

.372E-04 

- 

.416E-03 

. 398  E-0  3 

.  103E-03 

10 

.  192E-05 

- 

. 497E-04 

. 451E-04 

-  .239E-C5 

II 

- . 151E-06 

- 

. 692E-05 

.636E-05 

-  .  951E-06 

12 

-.606E-07 

- 

. 1 19E-05 

. 1 11E-05 

-  .  422E-06 

13 

- . 349E-07 

- 

. 207  E- 06 

. 255E-06 

-.  1 1 7  E- Ob 

14 

150E-07 

- 

, 433E-07 

.  528E-07 

-  .  342E-07 

15 

-.336E-08 

- 

. 776E-08 

. 112E-07 

-  .  626E-08 

STAGGERED 

COEFFICIENTS 

:  TE 

INCIDENCE 

N 

CCN) 

DC  N ) 

1 

144E-01 

.902E-01 

.  224E+00 

-  .  104E+00 

2 

- . 195E+00 

.115E-01 

123E+00 

-  .  955E-02 

3 

.  157E+00 

. 103E+00 

- . 850E-01 

.  994E-01 

4 

.  131E+00 

- 

. 138E+00 

. 696E-01 

.  928E-01 

5 

-  .  954E-01 

- 

. 472E-01 

.  158E+00 

-  .670E-01 

6 

-  .  6 16  E-0 1 

- 

. 168E-0I 

.659E-01 

- . 5 1 0  E-0 1 

37 


7 

8 
9 

10 
1  1 
12 
1  3 

14 

15 


-.171E-01  -.361E-02 
- . 3  06  E- 0  2  -.466E-03 
-.416E-03  - . 37  2E-04 
- . 497  E-C4  - . 192E-05 
- . 6  9 2£- 0 5  .151E-06 
-.119E-05  .oOoC-07 
- . 207E-06  .3495-07 
-.433E-07  .150E-07 
-.776E-03  .336E-08 


■135E-01  -.159E-01 

. loOE-02  -.29SE-02 

.103E-03  -.39SE-03 

-.2395-05  -.451E-04 

-  ■  95  i  E-06  -  .  s  5o  E-0  5 

-.422E-06  -.1115-05 

-.1175-06  -.255E-06 

-.3425-07  - . 528  E- 0  7 

-  .  6  26  E- 08  -. 1 1 2  E  —  0  7 
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INTERIOR  AND  SCATTERED  FIELD  EXPANSION  COEFFICIENTS 


AZIMUTHAL  MODE  NUMBER  M=  1 
INCIDENCE  ANGLE  ALPHA  =  -30.  DEGREES 


*******  INTERIOR  COEFFICIENTS 

******* 

STAGGERED 

COEFFICIENTS 

:  TM 

INCIDENCE 

N 

ACN) 

B  ( N  3 

1 

.457E+00 

. 212E+00 

.  555E-02 

.  83QE-01 

2 

.913E-02 

- 

. 474E-01 

-  .  246E+00 

-  .  152E+00 

3 

.  354E-01 

. 375E-01 

- . 262E-01 

.  429E-01 

4 

.260E-03 

- 

. 239E-02 

- . 538E-01 

- . 703E-01 

5 

.650E-02 

. 110E-01 

-.653E-01 

- . 473E-01 

6 

.  782E-02 

■751E-02 

. 204E-02 

-  .  190E-01 

7 

.  198E-02 

. 290E-02 

.651E-02 

-  .  626E-02 

8 

.  7 08E-Q3 

. 104E-02 

, 371E-02 

- . 27  7  E- 02 

STAGGERED 

COEFFICIENTS 

:  TE 

INCIDENCE 

N 

ACN) 

BIN) 

1 

.662E-01 

- 

. 344E-01 

. 624E+00 

.  267E-01 

2 

.  294E+00 

. 353E-01 

. 397E-01 

-  .  579E-01 

3 

.  833E- 0 1 

- 

. 504E-01 

. 342E-01 

.  253E-01 

4 

.  445E-01 

. 346E-01 

. 219E-01 

-. 115E-01 

5 

.653E-01 

.237E-01 

-. 173E-02 

.  1 19E-01 

6 

- .  126E-01 

. 127E-01 

- . 561E-03 

.  916E-02 

7 

151E-01 

. 504E-02 

- . 896E-03 

.  338E-02 

3 

- . 574E-02 

. 1996-02 

- . 384E-03 

.  100E-02 

*******  SCATTERED 

COEFFICIENTS 

******* 

STAGGERED 

COEFFICIENTS 

:  TM 

INCIDENCE 

N 

CCN) 

D(N) 

1 

.901E-01 

- 

.2296+00 

. 536E-01 

-  .  5 1 7  E  —  0 1 

2 

- . 376E-01 

.2006-01 

- . 105E+00 

.  169E+00 

3 

.411E-01 

- 

.3186-01 

.400E-01 

-.698E-02 

4 

-.387E-03 

- 

. b  06  6-0  3 

-.847E-01 

.  700E-01 

5 

. 127E-01 

- 

.368E-02 

- . 342E-01 

.  348E-01 

6 

. 497E-02 

- 

. 40SE-02 

- . S95E-D2 

-  .  529  E- 02 

7 

.  103E-02 

- 

.306E-03 

-. 195E-02 

-  .  393E-02 

8 

. 164E-03 

- 

. 112E-03 

- . 297E-03 

-  .  105E-02 

9 

.  124E-04 

- 

. 197E-05 

-.402E-04 

-  .  181E-03 

10 

.  339E-05 

. 149E-06 

- . 346E-05 

-  .  26  5 E- 0  4 

1  1 

- . 352E-Q6 

. 528  E- 06 

- . 403E-06 

-  .  292E-05 

12 

.  1 36  E- 06 

. 523E-07 

- . 977E-07 

- . 256  E- 06 

13 

- . 314E-07 

. 132E-07 

.357E-03 

- . 139E-07 

14 

.795E-Q8 

. 301E-08 

- . b02E-03 

-  .  652E-09 

15 

- . 137  E-OS 

- 

. 2S2E-09 

. 730E-09 

-  .  169E-09 

STAGGERED 

COEFFICIENTS 

:  TE 

INCIDENCE 

N 

CCN) 

D(N) 

1 

-.758E-02 

. 291E-01 

. 722E-02 

.  198E-03 

2 

. 260E-01 

- 

.267F-01 

- . 152E-01 

.  294E-01 

3 

532E-01 

. 473E-01 

. 212E-01 

-  .  237E-01 

4 

.726E-01 

- 

.648E-01 

- . 108E-01 

-  .  879E-02 

5 

. 301E-01 

- 

. 344E-01 

. 157E-01 

.  513E-Q2 

6 

. 763E-02 

. 348E-02 

.667E-02 

.  194E-02 
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7 

.  162E-02 

. 324E-02 

.  143E-02 

. 572E-03 

3 

.  236  E  —  0  3 

. 343  E-0  3 

.  222E-03 

. 1 14E-Q3 

9 

.  330E-04 

. 135E-03 

.  137E-04 

.261E-04 

10 

. 241 E- 05 

. ISoE-04 

.  394E-05 

. 294E-05 

1 1 

. 301E-Q6 

. 171E-05 

-  .  36  3E-06 

. 3  32 E-0 6 

12 

. 53  3 E- 07 

.  955E-07 

.  1 7 6  E- 06 

. 7  03  E-07 

1  3 

-.925E-03 

-  .  222E-03 

374E-07 

.115  E- Q  7 

14 

. 457  E-03 

-  .  20oE-08 

.  7  57E-0S 

.  1 7  b  E  -  0  3 

15 

-.360E-09 

-  .  157E-09 

-  .  154E-0S 

-  .  52SE-09 

90 


ACCUMULATED  FORWARD  AND  BACK  SCATTERED  AMPLITUDES 


M  =  1 

STAGGERED 

AMPS :  INCIDENT  TM 

ALPHA 

(DEG) 

0  . 

90  . 

ISO  . 

FWD 

- . 805E+01  -.177E+02 
-.294E+01  - . 359E+01 
-  .  776E+01  - . 17  2E+02 

139E+01 
. 200E+01 
- . 371E+01 

BK 

550E+01 
.  141E  +  01 
.  474E+01 

STAGGERED 

AMPS-  INCIDENT  TE 

ALPHA 

(DEG) 

0  . 

90  . 

ISO  . 

FWD 

-  .  805E+01  -.177E+02 
-.213E+01  - . 295E+01 
-.776E+01  - .  172E+02 

-  .  139E+01 
.  196E-01 
-.371E+01 

BK 

- . 550  E+  0 1 
-.247E+00 
.474E+01 

SCATTERING  CROSS  SECTION  CONVERGENCE 


M=  1 
NINC  3  1 


N 

NEW  TERM 

PCT  CHANGE 

1 

6E-01 

100 . 0000 

n 

<- 

.  2E+00 

80.6270 

3 

. 7E+00 

69  .  3592 

4 

.  1E+01 

58.7450 

5 

.2E+01 

47  .  0760 

6 

- 1 £+01 

17 . 3336 

7 

. 1E+0Q 

1.7682 

8 

.  4E-02 

.0730 

9 

.  9E-04 

.0017 

10 

.  2E-05 

.  0000 

11 

.  4E-07 

.  0000 

12 

.  2E-08 

.0000 

13 

•IE-09 

.  0000 

14 

. 6  E-l  1 

.0000 

15 

.  3E-12 

.  0000 

NINC  =  2 

N 

NEW  TERM 

PCT  CHANGE 

1 

.  6E-01 

100  .  0000 

2 

.  2E+00 

77.2227 

3 

.6E-01 

18.8104 

4 

.  3E+00 

52 . 9934 

5 

. 1 E+00 

17  .4571 

6 

. IE-01 

1.6261 

7 

.  3E-02 

.  3495 

3 

.  2E-03 

.  030  1 

9 

.  9E-05 

.0012 

10 

.  3E-06 

.0000 

11 

.  4E-08 

.  0000 

12 

. 7  E-l 0 

.  0000 

13 

.  IE-11 

.  0000 

14 

.IE-12 

.  0000 

15 

.  3E-14 

.  0000 
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* 


NINC  = 


N 

NEW  TERM 

PCT 

1 

.  5E-01 

100 

-> 

<- 

.  3E+00 

85 

T 

o 

.  1E+01 

7  6 

4 

.  3E+01 

6  3 

5 

. 1E+01 

22 

6 

. 2E+00 

4 

7 

. IE-01 

8 

. 3E-03 

9 

. IE-04 

10 

.  4E-06 

11 

. 9E-08 

12 

.  4E-09 

13 

. 2E-10 

14 

.IE-11 

15 

. 5E-13 

NIMC  =  4 


N 

HEW  TERM 

PCT 

1 

.  6  E — 0 1 

100 

2 

.  2E+00 

80 

3 

.  7  E+0  0 

69 

4 

. 1E+01 

58 

5 

.  2E  +  01 

47 

6 

. 1E+01 

17 

7 

.  1E  +  00 

1 

8 

.  4E-02 

9 

. 9E-04 

10 

. 2E-05 

11 

.4 E-07 

12 

. 2E-08 

13 

•IE-09 

14 

.  6  E—  1 1 

15 

- 3E-12 

NINC  =  5 


N 

NEW  TERM 

PCT 

1 

.8E-03 

100 

2 

.IE-01 

93 

3 

. 9E-01 

87 

4 

.  3E+00 

73 

5 

. 1E+00 

24 

6 

.IE-01 

2 

7 

.  2E-02 

8 

•2E-03 

9 

•6E-05 

10 

. 1 £-06 

11 

.2E-08 

12 

. 3E-10 

13 

.IE-11 

14 

■8E-13 

15 

. 4E-14 

CHANGE 
0000 
.3396 
4532 
6074 
9305 
1769 
1771 
.  0063 
0002 
0000 
.  0000 
0000 
.0000 
.  0000 
.0000 


CHANGE 
.  0000 
6270 
.  3592 
,  7450 
.  076  0 
.  3856 
.7682 
.07  30 
.  0017 
.  0000 
.  0000 
.  0000 
.  0000 
.  0000 
.  0000 


CHANGE 
.  0000 
.3376 
.6869 
.  3983 
.8194 
.  0363 
.  4062 
.0320 
.0011 
.  0000 
.  0000 
.  0000 
.  0000 
.  0000 
.  0000 


modal  cross  section  evolution 
NORMALIZED  W.R.T.  WAVELENGTH  SQUARED 

INCIDENCE  ANGLE  =  Q.  DEGREES 


M 

0 

INCIDENT  TM  POLARIZATION 

INCIDENT 

TE  POLARIZATION 

SCT 

. OOOE+OO 

ABS 

.  OOOE+OO 

TOT 

. OOOE+OO 

SCT 

. OOOE+OO 

ABS 

. OOOE+OO 

TOT 

. OOOE+OO 

1 

. 558E+01 

.  705E-01 

•565E+01 

. 558E+01 

. 705E-01 

.  565E+01 

2 

. 558E+01 

. 7  05E-01 

.  565E+01 

. 558E+01 

•705E-01 

.  5&5E  +  01 

3 

.  558E+01 

. 7  05E-0I 

.  565E+OI 

. 558E+01 

■705E-01 

565E+01 

4 

. 558E+01 

. 7  05E-01 

.  565£i 01 

. 558E+01 

•705E-01 

. S65E+01 

5 

.  558E+01 

. 705E-01 

■  565E+01 

. 558E+01 

•705E-01 

. 565E+01 

6 

•  558E+01 

. 7  05E-01 

■  565E+01 

. 558E+01 

•705E-01 

•565E+01 

7 

. 558E+01 

.  7  05E-01 

. 56  5E+01 

.  558  E-*-Cl 

.705E-0I 

•  565E+01 

INCI 

DENCE  ANGLE 

5  90.  DEGREES 

M 

0 

INCIDENT 

TM  POLARIZATION 

INCIDENT 

TE  POLARIZATION 

SCT 

■247E+00 

ABS 

■ 100E+00 

TOT 

.  347E+00 

SCT 

. 326  EAQ0 

ABS 

.  105E+00 

TOT 

•  4  31 E+  0  0 

1 

.  104E+01 

•  998 E- 01 

.  114E+01 

• 834E+00 

.  105E+00 

.  940E+00 

2 

. I 97E+G1 

. I10E+OO 

•  208E+01 

. 202E+01 

.  110E+0Q 

•  21 3  E+ 0  1 

3 

.  317E+0I 

. 106E+00 

. 328E+0I 

. 285E+01 

.  106E+00 

. 296E+01 

4 

•  440E+01 

1 13E+00 

. 45IE+01 

. 396E+0I 

. I 18 E+0 0 

. 407  E+  0  J 

5 

■  521 E+0 1 

.  144E+00 

. 536E+01 

• 5I9E+01 

• 199E+00 

■  539E  +  01 

6 

. 525E+01 

.  166E+00 

•542E+0I 

. 523E+0I 

■216E+00 

•545E+01 

7 

. 525E+01 

. 166E+Q0 

.  54ZE+01 

• 523E+0 1 

•  21 6 E+ 0  0 

.  54  5E+  0 1 
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INCIDENCE  ANGLE  =  180.  DEGREES 

INCIDENT  TM  POLARIZATION 


INCIDENT  TE  POLARIZATION 


M 

SCT 

ABS 

TOT 

SCT 

ABS 

TOT 

0 

. 220E-12 

- . 220E-12 

. OOOE+OO 

. 2Q9E-12 

-  .  209E-12 

.OOOE+OO 

1 

.546E+01 

. 255E-0 1 

. 548E+01 

. 546E+01 

.  255E-01 

. 546E+01 

2 

. 546E+01 

•255E-01 

. 548E+01 

.546E+01 

. 255E-01 

. 548E+01 

3 

. 546E+01 

. 255E-01 

. 543E+01 

•546E+01 

. 255E-01 

. 548E+  01 

4 

.  546  E-<  0 1 

.  255C-U1 

. 548E+01 

. 546E+01 

. 255E-0I 

.  548 E -i  01 

5 

•546E+01 

.255E-01 

. 548E+01 

. 546E+01 

. 255E-0 1 

. 548E+  01 

6 

•546E+01 

. 255E-01 

. 548E+01 

. 546E+01 

•255E-0I 

. 548E+01 

7 

.546E+01 

.  255E-01 

. 548E+01 

. 546E+01 

.255E-01 

. 548E+01 

I 

I 
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*****  scattered  field  complex  amplitudes  ***** 


INCIDENCE  ANGLE  = 
PHI =  90.  DEGREES 

STAGGERED  RAD  FL  DS ; 


0.  DEGREES 


TM  INCIDENCE 


THETA  F-THETA  F-PHI 


MAG 

DEG 

MAG 

DEG 

0  . 

. 851E-06 

.  00 

.195E+02 

65.60 

5  . 

. 821E-Q6 

.00 

.  187E+02 

65.19 

IQ  . 

. 735E-06 

.  00 

. 166  E* 02 

63. 93 

15. 

.603E-06 

.  00 

.  133E+02 

61.71 

20  . 

.461E-Q6 

.  00 

. 925E+01 

58 . 17 

25  . 

. 322E-06 

.00 

. 499E+01 

51.69 

30  . 

. 229E-06 

.00 

.  114E+01 

19.66 

35. 

•207E-06 

.00 

. 253E+01 

-112.81 

40. 

. 225E-06 

.  00 

.  484E+0 1 

-124.99 

45. 

. 239E-06 

.  00 

.6Q3E+01 

-131.99 

50. 

. 235E-06 

.00 

.60 8E+01 

-138.30 

55. 

. 217E-06 

.00 

,  515E+01 

-144.73 

60  . 

. 197E-06 

.  00 

.  353E+01 

-151.82 

65. 

. 183E-06 

.00 

.  159E+01 

-161 .86 

70  . 

. 173E-06 

.  00 

.  382E+00 

49 . 46 

75. 

. 160E-06 

.  00 

.  187E+01 

18 . 57 

80  . 

. 138E-06 

.  00 

.281E+01 

10.27 

85. 

-111 E— 06 

.  00 

.  307E+01 

3.63 

90. 

.933E-07 

.00 

. 27  QE+01 

-2.13 

95. 

.993E-07 

.00 

.  186E+01 

-6.50 

loo . 

•120E-06 

.00 

.812E+00 

-5.21 

105. 

•137E-06 

.00 

. 29GE+00 

113.65 

no. 

■138E-06 

.  00 

.  996E+00 

140.86 

115. 

. 121E-06 

.  00 

.  133E+01 

138.49 

120  . 

.889E-07 

.  00 

.  121E+01 

133.24 

125. 

.557E-07 

.  00 

.716E+00 

123.28 

130. 

. 533E-07 

,  00 

,  189E+00 

39 . 52 

135. 

. 939E-07 

.00 

. 7  54E+Q0 

-30 .93 

140. 

. 129E-06 

.  00 

.  125E+01 

-35 . 32 

145. 

153E-06 

.00 

.  142E+01 

-29 .65 

150  . 

. 165E-06 

,00 

.  134E+01 

-9.95 

155. 

. 173E-06 

.00 

.  149E+01 

26 . 89 

160. 

. 184E-06 

.  00 

.231E+Q1 

55.01 

165. 

. 202E-06 

.  00 

.348E+01 

67  .  53 

170  . 

. 224E-06 

.  00 

.460E+01 

72.91 

175. 

•241E-06 

.  00 

.  539E+G 1 

75.20 

130. 

.248E-06 

.00 

.568E+01 

75.85 
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*****  SCATTERED  FIELD  COMPLEX  AMPLITUDES  ***** 

INCIDENCE  ANGLE  =  0.  DEGREES 

PHI =  90.  DEGREES 

STAGGERED  RAD  FLDS  =  TE  INCIDENCE 


THETA  F-THETA 


F-PHI 


MAG  DEG 


0  . 

.  195E+02 

-114. 

40 

5. 

.  188E+02 

-115. 

57 

10  . 

.  168E+02 

-119. 

26 

15. 

.  1 39E+02 

-126  . 

13 

20  . 

.  105E+02 

-j.37  . 

69 

25. 

.  7  38E+0 1 

-157  . 

34 

30  . 

.  524E+01 

169  . 

80 

35. 

.  473E+0 1 

128  . 

60 

40  . 

.  515E+01 

95. 

19 

•45. 

.  547  £-*-01 

69  . 

76 

50  . 

. 537E+U1 

46  . 

52 

55. 

.497E+01 

21  . 

84 

60  . 

. 451 £+0  1 

-5. 

65 

65. 

. 419E+01 

-34. 

83 

70  . 

.396E+01 

-63. 

23 

75. 

.366E+01 

-90  . 

.06 

80. 

. 316E+01 

-117 

.25 

85. 

. 255E+01 

-149 

.31 

90. 

.214E+01 

168 

.63 

95. 

.227E+01 

123 

.71 

100. 

. 275E+01 

89 

,  95 

105. 

. 314E+01 

65 

.78 

110. 

. 317E+01 

45 

.  98 

115. 

. 277E+01 

26 

.74 

120  . 

.203E+01 

3 

.  03 

125. 

.127E+01 

-39 

.  54 

130  . 

. 133E+01 

-107 

.69 

135. 

. 215E+01 

-147 

.  00 

140  . 

. 295E+0 1 

-169 

.  25 

145  . 

.  349E+01 

172 

.65 

150  . 

. 378E+01 

154 

.61 

155. 

.396E+01 

135 

.25 

160. 

.421E+01 

115 

.51 

165  . 

.462E+01 

98 

.  19 

170  . 

.512E+01 

85 

.58 

175. 

. 552E+0 1 

78 

.  23 

180  . 

. 568E+01 

75 

.85 

MAG 

.  851E-06 
.819E-06 
.724E-06 
. 58  OE-06 
. 404E-06 
.218E-06 
. 497E-07 
. 1 I 0E-06 
.2I2E-06 
. 26  4E-06 
.266E-06 
. 225E-06 
. 154E-06 
. 693E-07 
. 167E-07 
.818E-07 
. 123E-06 
. 134E-06 
.  1 18E-06 
.8I4E-07 
.  355E-07 
.127E-07 
. 435E-07 
.  579E-07 
.527E-07 
.  313E-07 
. 828E-08 
.329E-07 
.  545E-07 
.619E-07 
.  585E-07 
. 6  53E-07 
. 101E-06 
.152E-06 
.  201E-06 
.236E-06 
.  248E-06 


DEG 
.  00 
.00 
.00 
.00 
.  00 
.00 
.00 
.00 
.  00 
.  00 
.  00 
.  00 
.00 
.  00 
.00 
.00 
.  00 
.00 
.00 
.00 
.  00 
.  00 
.  00 
.00 
.00 
.  00 
.  00 
.  00 
.  00 
.  00 
.  00 
.  00 
.  00 
.00 
.  00 
.00 
.  00 
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nnnnonnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnoonnn 


APPENDIX  C 
EMCAD  SOURCE  CODE 


PROGRAM  EMCAD 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

THIS  PROGRAM  COMPUTES  EM  FIELD  SCATTERING  FROM  LAYERED 
BODIES  OF  REVOLUTION  USING  THE  COUPLED  AZIMUTHAL  POTENTIAL 
(CAP)  FORMULATION  IN  CONJUNCTION  WITH  A  VARIATIONAL  FINITE- 
ELEMENT  TECHNIQUE  AND  A  TRI-REGIQNAL  UNIMOMENT  METHOD.  THE 
NUMERICAL  SOLUTION  IS  PERFORMED  USING  A  TWO-SWEEP  RICCATI 
TRANSFORM  WITH  CIRCUMFERENTIAL  MARCHING.  UNIT  1  IS  USED  FOR 
SEQUENTIAL  UNFORMATTED  STORAGE  OF  GENERATED  RICCATI  ARRAYS  TO 
DISK.  UPDATED  VERSION  OF  CDC-7600  U.C.  BERKELEY  CODE  "EMSCAT" 

WITH  LAYERED  MEDIA  ADDITIONS. 

MODS  BY  M . A .  MORGAN  MAR  87-FEB  88 
E.M.  CONNOLLY  MAR  87-FEB  88. 

xxxxxxxxxxxxxxxxxxxx  EMCAD  PARAMETERS  xxxxxxxxxxxxxxxxxxxxxxxxx 

xx  NA  -  NO.  OF  INCIDENT  ANGLES,  DALPHA(J)  (DEGREES) 

XX  DMIN, DMAX  -  MIN  AND  MAX  MESH  DENSITY  ( EL  EM/LAMBDA ) 

xx  NR1 , NR2  -  RADIAL  NODE  SPACING  BETWEEN  RA-R1  8  R2-RB 

xx  NLAY  -  NO.  OF  LAYERED  SURFACES,  NLAY+1  IS  FREE  SPACE 

XX  ER(NLAY),UR(NLAY)  -  COMPLEX  ER  AND  UR  OF  PENETRABLE  LAYERS 

xx  MSTART ,MSTOP  -  AZIMUTHAL  MODE  RANGE  ON  "M” 

xx  N1.N2  -  UPPER  LIMITS  ON  INTERIOR  8  EXTERIOR  MODES  "N" 

xx  NT , NP  -  NO.  OF  FIELD  PRINTOUT  THETA  8  PHI  POINTS 

XX  NB , NBMX  -  NO.  OF  INTERIOR  NODAL  POTENTIALS/THETA  STEP 

XX  NCMAX  -  NO.  OF  INTERIOR  +  EXTERIOR  EXPANSION  COEF'S 

XX  NMMAX  -  NCMAX  +  2*NA  (NO.  UNIMOMENT  BC’S) 

xx  NINC, NINMX  -  2XNA  (2  POLARIZATIONS/INC  ANGLE) 

XX  BC  -  0  OR  1  :  Z-AXIS  HOMO  DIRICHLET  OR  CAUCHY  BC'S 
XX  I MX ,  I B IG  -  NO.  OF  THETA  STEPS  IN  MARCHING  8  DATA  FILE 

XXXXXXXXXXXXXXXXXXXX  ARRAY  DIMENSIONS  xxxxxxxxxxxxxxxxxxxxxxxxxx 

XX  R , RI NV  -  (NBMX, NBMX)  RICCATI  MATICES 

XX  S , SINV  -  (NBMX, NMMAX)  S-VECTOR  ARRAYS 

XX  FMOM, DUM  -  (NCMAX, NCMAX)  COEFFICIENT  MOMENT  ARRAYS 

XX  DMOM.COEF  -  (NCMAX, NINMX)  DRIVING  8  COEF  ARRAYS  (SAME) 

XX  FT  1  TO  FP2  -  ( NT , NP , NA )  SCATTERED  E-FIELDS 

xx  RAD  -  (2,5)  RADII(THETA)  OF  LAYERS  AT  I  8  1+1  STEPS 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

CHARACTERX1  BEL L , FF , L EVEL 
CHARACTERX8  GR FDA T , PRT DAT 

CHARACTERX12  GRFDT1 , GRFDT2 , GRFDT 3 , GRFDT4 , PRTDT 1 
CHARACTERX64  DUMMY , GRA FL AB , STRDA T , HDR , MATPAR 
COMPLEX  A, B.G,P,E1 ,H1 , E2, H2, EB1 ,HB1 , EB2.HB2 
COMPLEX  F30,F21,F12,F20,F1I,F10,F1,F2,FT1,FP1,FT2,FP2 
COMPLEX  R(58,58),S(58,80),RINV(58,58),SINV(58,80),ER(6),UR(6) 
COMPLEX  FMOM ( 7  0,70),COEF(70,10) 

INTEGER  BC, ELND, IDUMMY 
REAL  RAD( 2 , 5 ) 

COMMON/ ONE/X( 62 ),Y( 62) 

COMMON/ TWO/ EL  ND( 6  0,3),NDEL(62,3),NDCT(62) 

COMMON/ THREE/ EBI(3,80),HB1(3,80),EB2(3,80),HB2(3,80) 
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COMMON/ FOUR/A  C 58,4),BC58,6),GC58,4),PC4,80).KA(58),KBC58),KGC58) 
COMMON/FI  VE/E1  C60,3,3),HlC60,3,3),E2C60,3,3),H2Cci0,3,3) 

COMMON/ SI  X/ FT 1 ( 57 ,3, 5 ) . FP 1 ( 57 , S , 5 ) , FT2 ( 57 , 3  >  5 ) . FP2 ( 57 , 8 , 5  ) 
COMMON/S  EVEN/  SCATC 13, 10),ABSBC13.10),TOTLC15,10) 

COMMON/ EIGHT/ ALPHA! 5 ) . PH  I C 3 ) , DAL PHA C  5 ) , DPH I C 3 ) 

COMMON/TEN/ F30 , F21 , F12 , F20 , FI  1 , FI  0 , P30 , P21 , P 12 , P20 , PI 1 , P 1 G 
EQUIVALENCE  CRINV.SINV) 

BELL -CHARI  7 ) 

FF  =  CHAR ( 12) 

SETTING  MAXIMUM  MATRIX  DIMENSIONS 

NCMAX=7  0 

NMMAX  =  8  0 

NBMX=53 

NINMX3 1 0 

PI  =  3. 1415927 

DTR^PI/180  . 

READING  INC.  FIELD,  MESH,  AND  SOLUTION  PARAMETERS 
HRITEC  *, 137 ) 

READC  *,101)  MATPAR 
0PENC2,FILE=MATPAR) 

HRITEC*, 112) 

READC  * , 1 0 1 )  STRDAT 
OPENC 10 , FILE'STRDAT) 

READC  2,101)  HDR 
HRITEC*, 101) 

HRITEC*, 101)  HDR 
READC  2 , 1 0 1 )  HDR 
WRITEC  *, 101 )  HDR 
HRITEC*, 25) 

READC*, 101)  DUMMY 
HRITEC*, 101)  FF 
HRITEC*, 10) 

READC*, 101)  DUMMY 
HRITEC*, 101)  FF 
HRITEC*, 12) 

READC*, 101)  LEVEL 
HRITEC*, 101)  FF 
IFC  LEVEL . EQ.  '3'  )  GO  TO  99 
IFCLEVEL  EQ. '2'  )  THEN 
HRITEC*, 1131 
READC*. 101)  PRTDAT 
ELSE 

HRITEC*, 13) 

READC*, 101)  PRTDAT 
ENDIF 

PRTDT1  =  PRTDAT//' .OUT' 

0PENC3,FILE=PRTDT1) 

IFCLEVEL  EQ. • 2 ' )  THEN 
HRITEC*, 114) 

READC*, 101)  GRFDAT 
ELSE 

HRITEC*, 14) 

READC*, 101)  GRFDAT 
ENDIF 

GRFDT1  =  GRFDAT//' . TMT ' 

GRFDT2  =  GRFDAT//' . TMP ' 

GRFDT3  =  GRFDAT//' . TET' 

GRFDT4  =  GRFDAT//' .TEP' 

OPENC  4 , FIL  E=GRFDT1 ) 

0PEN(7,FILE=GRFDT2) 

0PENC8 , FI LE=GRFDT  3 ) 

OPENC  9 , FILE=GRFDT4 ) 

IFCLEVEL. EQ. '2')  THEN 
HRITEC*, 115) 

READC*, 101)  GRAFLAB 
ELSE 


MRITEC*. 15) 

READC*, 101)  GRAFLAB 
EHDIF 

HR I TEC  4 , 1 01 )  GRAFLA8 
HR  I  TEC  7,101)  GRAFLAB 
WRITEC3, 101 )  GRAFLAB 
MR  I  TEC  9 , 101  )  GRAFLAB 
I FC  LEVEL .EQ. >2' )  THEN 
HR  I  TEC  *,116  ) 

READC*,*)  DMIN 
MRITEC*, 117) 

READ  C  * ,  * )  DMAX 
ELSE 

MRITEC*, 16) 

READC  * ,  * )  DMIN 
HR  I  TEC  * ,  17  ) 

READC*,*)  DMAX 
ENDIF 

HR I  TEC  3,101)  FF 
MRITEC  3, 100 ) 

HRITEC  3, * ) 

HRITEC3, 111) 

READC2, 102)  NLAY 
DO  22  L  =  1 , NLAY 
READC  2,103)  ERC  L ) , URC  L  ) 
WRIT  EC  3 , 1 0  7  )  L , ERC  L  ) , URC  L  ) 
22  CONTINUE 

READC 10,138)  IBIG 
READC 10,103)  RMIN , RMAX 
IFCLEVEL . EQ. '2' )  THEN 
MRITEC  *, 118  ) 

READC*,*)  NA 
ELSE 

MRITEC*, 18) 

READC*,*)  NA 
ENDIF 
NINC=2*NA 
SAMX=0 . 0 
MRITEC  3, 109) 

DO  33  1=1, NA 
MRITEC*, 23)  I 
READC*.*)  DALPHACI) 
WRITEC3.110)  I, DALPHACI) 
ALP HA CI)=DTR* DALPHACI) 

SA  =  SINCALPHACI)  ) 
IFCSA.GE.SAMX)  SAMX=SA 
33  CONTINUE 

IFCLEVEL . EQ. '2')  THEN 
MRITEC*, 119) 

READC*,*)  NT 
ELSE 

MRITEC*, 20) 

READC*,*)  NT 
ENDIF 

IFCLEVEL . EQ . '2' )  THEN 
MRITEC*, 120) 

READC*,*)  NP 
ELSE 

MRITEC*, 21) 

READC*,*)  NP 
ENDIF 

DO  44  1=1, NP 
MRITEC*, 24)  I 
READC*,*)  DPHICI) 

44  PHICI)=DTR* DPHICI) 

MSTART  =  0 
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MST0P=INTCRMAX*SAMX)+1 
WRITEC*, 121 )  MSTOP 
REA  DC  * , * )  MSTOP 
IFCSAMX.LT. 1 . E-6 )  MSTART-1 

C  SETTING  RADIAL  STEPS  BETWEEN  RA  to  RI  AND  R2  to  RB 

NR1  =  1 
NR2  =  1 

C  DEFINING  OUTSIDE  PERMITTIVITY  WITH  SLIGHT  LOSS 

URCNLAY+1  )  =  C 1 . 0, -1 . 0  E-6 ) 

ERCNLAY+1 )=C 1 . 0, -1 . OE-6 ) 

C  ESTIMATING  NUMBER  OF  CORE  AND  EXTERIOR  MODES 
N1=INT(RMIN*CABSCCSQRT(2 . 0*ER( 1 )*URC 1 ) ) ) )  +  l 
N2=INTCRMAX*CABS(CSQRTC2. 0*ERC NLAY  )  *URC NLAY ) ) ) )  +  l 
WRITEC*, 122)  N1 
REA  D ( * , * )  Ni 
WRITEC  * , 123 )  N2 
READ!*,*)  N2 
WRITEC  * , * )  BELL 

PAUSE  'Check  that  N1+N2  is  .LE.  35  ....  Otherwise  Abort* 

IFCCN1+N2)  . GT  .  35  )  STOP 

WRITEC  3.105)  MSTART, MSTOP, N1,N2 

WRITEC  *,124) 

Ml =MSTART  + 1 
M2=MSTCP+1 
MM =MSTOP~MS TART-1 
IFCMM.LE.O)  MM=1 
DO  77  MD=M1,M2 

OPENC1 . FILE='RS. DAT' , STATUS= ' UNKNOWN ' , ACC ESS = ' SEQUENTIAL 
1  FORM= ' UNFORMATTED ' ) 

DD=DMAX-DMIN 

M=MD-1 

WRITEC  3> 101 )  FF 
WRITEC3,204)  M 
WRITEC*,*)  BELL 
WRITEC*, 125)  M 
BC  =  0 

IF  CM.EQ.i)  BC  =  1 
MF  -  0 

IFCM.EQ.O)  MF -1 

C  COMPUTING  MODAL  EXPANSION  TRUNCATION  LIMITS 

NM1 =N1+1 -M-MF 
IFCNM1.LT. 1)  NM1 = 1 
NM2-N2+1-M-MF 
I F( NM2 . L T . 1 )  NM2- 1 
NCOEF-2*  (  NM1-*  NM2  ) 

NMODES=NCQEF+NINC 

C  COMPUTING  MESH  PARAMETERS  AND  RADCI)  FOR  GIVEN  MODE  -  M 
DENS=DMAX-DD*CM+MF-1 . 0)/MM 
WRITEC*, 126) 

CALL  MESH  C DENS , ER , UR , NR1 , NR2 , IMX , NDS , LMX, NB , DR, DT, 

1  RA , Rl , R2 , RB , RMIN , RMAX , NLAY ) 

C  CHECKING  MAXIMUM  MESH  PARAMETERS 
IFCNB.LE.NBMX)  GO  TO  55 
WRITEC3, 104) 

STOP 

55  CONTINUE 

C  LOADING  AND  MARCHING  CCR))  AND  CCS)) 

WRITEC  3, 101 )  FF 
WRITEC3,201) 

WRITEC*, 127)  IMX 
DO  66  1=1, IMX 
IR=IMX-I 

WRITEC*, 128)  I , IR 
WRITEC*, 129) 

CALL  LODER  C I , M, BC, ER , UR, NMODES , NM1 , NM2, NINC , NB , IMX, I  BIG 
1  NDS,RA,RB,R1,R2, DR, DT , NR1 , NR2, RAD, NLAY ) 
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IF(I.GT.l)  GO  TO  57 

WRITE!  3, 30  0  )  I , ( RAD( 1 , L ) , L  =  1 , NlAY ) 

57  CONTINUE 

WRITE! 3, 300)  I , ( RADC 2,  L ) , L =1 , NLAY ) 

WRITE! *.130) 

CALL  MARCH  ( I , IMX , BC , MB , NBMX, NMODES , NMMAX , R , S , RI NV ) 

66  CONTINUE 

WRITE!*. 131) 

CALL  SWEEP  (M, IMX, BC. NM1 , NM2 , NINC , NMODES , NMMAX, NCOEF, NCMAX, 

1  NR1 , NR2, NB , NBMX, R1 , R2, DT , ER, UR , R, S , SINV , FMOM, COEF, NINMX  ) 
WRITE!*, 132) 

CALL  MODECO  ( NM1 , NM2, NINC, NCOEF, NCMAX, FMOM, COEF, M, NINMX ) 

WRITE!*, 133) 

CALL  RADFLD  (M, NM1 , NM2, NT , NP, NA, NCMAX, COEF, NINMX ) 

WRITE!*, 134) 

CALL  CRSSEC  ! M, NM1 , NM2 , NA, NCMAX, COEF , NINMX) 

REWIND  10 

READ! 10, 138)  IBIG. IDUMMY 
READ! 10,103)  RMI N , RMAX 
77  CONTINUE 

CLOSE! 1 , STATUS3 ' DELETE ' ) 

CL  OSE! 2 ) 

CLOSE! 1 0  ) 

WRITE!*, 135) 

CALL  DATOUT  ! MSTART , MSTOP , NT , NP , NA ) 

WRITE!*,*)  BELL 
WRITE!*, 136) 

10  FORMAT (////Z////////////7X, '******  WELCOME  TO  EMCAD  ******', //7 X , 

1 'THIS  PROGRAM  COMPUTES  EM  FIELD  SCATTERING  FROM  PENETRABLE' ,/7X, 
2'BODIES  OF  REVOLUTION  USING  THE  COUPLED  AZIMUTHAL  POTENTIAL ' ,/7X, 

3 ' ! CAP )  FORMULATION  IN  CONJUNCTION  WITH  A  VARIATIONAL  F I N I T E- ' , / 7 X , 
4 ' ELEMENT  TECHNIQUE  AND  A  TRI-REGIONAL  UNIMOMENT  METHOD.  THE',/7X, 
5 ' NUMERICAL  SOLUTION  IS  PERFORMED  USING  A  TWO-SWEEP  RICCATI ' , /7 X , 

6 'TRANSFORM  WITH  CIRCUMFERENTIAL  MARCHING . ' , /7X , 

7 ' UPDATED  VERSION  OF  ORIGINAL  U.C.  BERKELEY  CODE.',/7X, 

8' MODIFICATIONS  BY  M . A .  MORGAN  MAR  1987  -  FEB  1988', /7X, 

9'  E.M.  CONNOLLY  JUL  1987  -  APR  1988  .  S//////7X, 

1 'PLEASE  PRESS  ANY  KEY  TO  CONTINUE.',/) 

12  FORMAT (//////////// ////7X, '*******  EMCAD  INPUT  ******»’, //7X , 

1 'EMCAD  ALLOWS  A  USER  TO  INPUT  THE  NECESSARY  INFORMATION' ,/7X, 

2 'ACCORDING  TO  THE  LEVEL  OF  HIS  EXPERTISE ,//7X, 

3 ' 1  NOVICE  LEVEL  -  GIVES  BRIEF  EXPLANATIONS/DESCRIPTIONS' ,/7X, 

A '  OF  INPUT  VALUES.  DESCRIBES  THE  REQUIRED' ,/7X, 

5'  FORMAT  FOR  THE  INPUT  VALUE.  GIVES  ',/7X, 

6'  TYPICAL  VALUES  WHERE  APPLICABLE. ',//7X, 

7 ' 2  EXPERT  LEVEL  -  ASSUMES  THE  USER  IS  FAMILIAR  WITH  THE',/7X, 

3'  INPUT  PARAMETERS  AND  FORMATS.  SIMPLY’, /7X, 

9'  PROMPTS  FOR  REQUIRED  INPUTS .', //7X , 

1 ' 3  EXIT  EMCAD', //7X, 

2 ' PLEASE  SELECT  THE  LEVEL  OF  EXPERTISE  BY  ENTERING  1,  2,  OR  3.', 

3//7X) 

13  FORMAT! ////////////////////// 7 X , 

1 ' THE  OUTPUT  DATA  FILE  IS  THE  FILE  CONTAINING  THE  OUTPUT', /7X, 

2 ' RESULTS  OF  ALL  NUMERICAL  CALCULATIONS  CONDUCTED  BY  EMCAD. ’,/7X, 

3 ' THE  FORMAT  FOR  THIS  INPUT  IS  FILENAME  ONLY.  NO  EXTENSION  IS',/7X, 
4 ’ REQUIRED  OR  DESIRED.  EMCAD  AUTOMATICALLY  APPENDS  AN  ',/7X, 

5 ' EXTENSION  OF  .OUT  TO  YOUR  FI L ENAME .  ' , / '7 X , 

6 '  PL  EASE  ENTER  THE  FILENAME  OF  T  r.Z  OUTPUT  DATA  FILE.',//7X) 

14  FORMAT !//////////////////////7X, 

1 ' THE  OUTPUT  GRAPHICS  DATA  FILE  IS  THE  FILE  CONTAINING  THE  ',/7X, 

2 ' OUTPUT  DATA  FROM  EMCAD  TO  BE  USED  AS  INPUT  TO  GRAPHING  ' ,/7X, 

3 ' ROUTINES .  THE  FORMAT  FOR  THIS  INPUT  IS  FILENAME  ONLY.  NO',/7X, 

4 ' EXTENSION  IS  REQUIRED  OR  DESIRED.  EMCAD  AUTOMATICALLY  ',/7X, 

5 'APPENDS  AN  EXTENSION  OF  . TMT ,  . TMP ,  .JET,  ‘HD  .  TEP  TO  YOUR'./7X, 

6 ' FI L  ENAME  AS  IT  PRODUCES  FOUR  OUTPUT  FILES  FOR  GRAPHICS .' ,//7X, 

7'. TMT  - >  TM  INCIDENCE,  F-THETA ' , /7 X , 


S'.TMP  - >  TM  INCIDENCE,  F-PHI',/7X, 

9'.TET  - >  TE  INCIDENCE,  F-THETA' ,/7X, 

l'.TET  - >  TE  INCIDENCE,  F-°HI',//7X, 

2 1  PL  EASE  ENTER  THE  FILE  NAME  OF  THE  OUTPUT  DATA  FILE.’,//7X) 

15  FORMAT (///////////////////////IX. 

1  1  THE  GRAPHICS  CAPTION  IS  A  PERSONALIZED  CAPTION  ALLOWING  THE’,/7 
2 1  US ER  TO  IDENTIFY  THIS  SET  OF  GRAPHS  FROM  ALL  OTHER  SETS.',/7X, 

3 ' THE  MAXIMUM  LENGTH  OF  THIS  CAPTION  IS  64  CHARACT ERS . ' , / 7 X . 

4 ' NOTE :  WHEN  USED  WITH  THE  GRAPHICS  PACKAGE,  THE  PROGRAM  IS',/7X, 
5'ABLE  TO  DIFFERENTIATE  BETWEEN  UPPER  CASE  AND  LOWER  CASE’,/7X, 

6  ’CHARACTERS .  ’ ,//7X, 

7'PLEASE  ENTER  ANY  GRAPHICS  CAPTION  YOU  DESI RE » , // 7 X ) 

16  FORMAT (//////////////////////7X, 

l’DMIN  AND  DMAX  ARE  PARAMETERS  OF  MESH  DENSITY  IN  TERMS  0F’,/7X, 
2’ ELEMENTS/ INTERIOR  LAMBDA.  INPUT  VALUES  ARE  EXPECTED  TO  BE’,/7X, 
3  ’  REAL ,  I.E.  THE  DECIMAL  POINT  MUST  BE  INCLUDED ,/7X, 

4’TYPICAL  VALUES  ARE  DMIN  =  10.  AND  DMAX  =  15.’,//7X, 

4 ’ PL  EASE  INPUT  DMIN',//7X) 

17  F0RMATC//7X, ’PLEASE  INPUT  DMAX’,//7X) 

IS  FORMAT (////////////////////// 7X, 

1  ’  THE  NUMBER  OF  INCIDENT  FIELD  ANGLES  IS  THE  TOTAL  NUMBER  OF’,/7) 

2  ’  I  NCI  DENT  FIELDS  THAT  xMPINGE  ON  THE  OBJECT  OF  INTEREST .’ ,/7X, 

3  ’  THIS  PROGRAM  ALLOWS  A  MAXIMUM  OF  FIVE  INCIDENT  FIELD  ’,/7X, 

4 ’ANGLES .  WHEN  ENTERING  YOUR  ANSWER  PLEASE  DO  NOT  INCLUDE  A’,/7) 
5 ’ DECIMAL  BECAUSE  THE  INPUT  MUST  BE  IN  INTEGER  FORMAT,  I.E.’,//7) 
5’  NA  =  3  ,  NA  =  1  ’ ,//7X, 

o ’ PL  EASE  INPUT  THE  NUMBER  OF  INCIDENT  FIELD  ANGL ES . ’ , //7X ) 

19  FORMAT (/////////// ///////////7X, 

1MHE  INCIDENT  ANGLES  DEFINE  THE  DIRECTION  FROM  WHICH  THE  ’,/7X, 

2 ’ INCIDENT  FIELDS  ORIGINATE .’, //7X , 

3 ’ 0  DEGREES  - >  +  Z  DIRECTION’ ,/7X, 

4 ’ 9 0  DEGREES  - >  +  X  DIRECTION ’ ,/7X, 

5*180  DEGREES  - >  -  Z  DIRECTION* ,//7X) 


20  FORMAT (/////// 


'///IX, 


1  *  THE  NUMBER  OF  SCATTERING  FIELD  THETA  POINTS  DETERMINES  THE’,/7> 
2 ' SPACING  BETWEEN  THETA  POINTS  DURING  EMCAD  ITERATIONS  AND’,/7X, 
3’CALCULATIONS. * ,//7X, 

4' DELTA  THETA  =  180  /  (NUMBER  THETA  POINTS  -  1)  SO...',//7X, 

5  *  NUMBER  THETA  POINTS  =  37  - >  DELTA  THETA  =  5  DEGREES ',//? X , 

6 '  NUMBER  THETA  POINTS  =  19  - >  DELTA  THETA  =  10  DEGREES’,//?)' 

7 ’WHEN  ENTERING  YOUR  ANSWER,  PLEASE  DO  NOT  INCLUDE  A  DECIMAL',/?? 
8 ' BECAUSE  THE  INPUT  MUST  BE  IN  INTEGER  FORMAT.  I.E.  NT  =  19’,///? 
9 'PLEASE  INPUT  THE  NUMBER  OF  SCATTERING  FIELD  THETA  POINTS ’,// 7 X J 
21  FORMAT (/////////////////// ///7X, 

l’THE  NUMBER  OF  PHI  ANGLES  IS  THE  TOTAL  NUMBER  OF  PHI  ANGLES.',/? 
2 ' Tills  PROGRAM  ALLOWS  A  MAXIMUM  OF  EIGHT  PHI  ANGLES.  WHEN  ’,/7X, 
3 ’ ENTERING  YOUR  ANSWER  PLEASE  DO  NOT  INCLUDE  A  DECIMAL  ',/7X, 

4 ' BECAUSE  THE  INPUT  MUST  BE  IN  INTEGER  FORMAT.  I.E.',//7X, 

5’  NP  =  3  ,  NP  =  1  '  ,//7X, 

6 ' PL  EASE  INPUT  THE  NUMBER  OF  PHI  ANGL ES . ' , //7X ) 

23  FORMATC/7X, 'ENTER  INC  FLD  ANGLE  (DEG)  FOR  *  *,I3,//7X) 

24  FORMATC/7X, 'ENTER  PHI  FLD  ANGLE  (DEG)  FOR  #  *,I3,//7X) 

25  FORMAT(//////7X, 'PLEASE  PRESS  ANY  KEY  TO  CONTINUE.',////) 

100  FORMAT!/, ' XXXXXXXXXXXXXXXXXXXX  EMCAD  OUTPUT  DATA  xxxxxxxxxxxxxx* 
1 XXXXXXXXXXX ' ) 

101  FORMAT ( A ) 

102  FORMAT (15) 

103  FORMAT ( 4 ( El 4 . 6 ) ) 

104  FORMAT ( '  xx*  PROGRAM  ABORTED  BECAUSE  NB  >  NBMX  ***') 

105  FORMAT(/7X, 'SERIES  SOLUTION  PARAMET ERS '/7X, ’ MSTART= ' , 1 5 , /7X , 
l'MSTOP=  *,I5,/7X, 'Nl=  ' , I8,/7X, 'N2=  ',18) 

107  FORMAT (//7X, I  5 , 2 ( 2 X , ' ( ' , 1  PEI  1 . 3 , 2X , 1PE1 1 . 3 , ' ) ' ) ) 

109  FORMAT(//7X, ' INCIDENT  FIELD  ANGL ES ’ , //9X , ' N ' , 1  OX, ' ALPHA( N ) ' ) 

110  FORMAT (//7X,I3,7X,F5.0, '  DEG') 

111  FORMAT(/7X, 'COMPLEX  Er(n)  AND  Ur(n)') 

112  FORMAT(//7X, 'ENTER  STRUCTURE  DATA  FILE  ( D ; FILENAME . EXTENSION ) : 
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113  F0RMATC//7X, 'ENTER  OUTPUT  DATA  FILENAME  WITHOUT  EXTENSION' ,//7X, 

1 ' FOR  PRINTER  OUTPUT  ENTER  "LPT1";  ') 

114  FORMAT(//7X, 'ENTER  GRAPHICS  DATA  FILENAME  WITHOUT  EXTENSION  ') 

115  FORMATC//7X, ' ENTER  CAPTION  OR  LABEL  ') 

llo  FORMATC//7X,  'ENTER  DMIN  ( ELEMENTS/  INTERIOR  LAMBDA)  -.  ') 

117  FORMATC//7X, ' ENTER  DMAX  ( EL EMENTS/ I  NT ERI OR  LAMBDA):  ') 

113  FORMAT(//7X, 'ENTER  NUMBER  OF  INCIDENT  FIELD  ANGLES  (.LE.  5)  ') 

119  FORMAT(//7X, ' ENTER  NUMBER  OF  THETA  POINTS  FOR  SCAT  FIELD  (.LE  37.) 
1 ' , 4X) 

120  FORMATC//7X, ' ENTER  NUMBER  OF  PHI  ANGLES  (.LE.  8):  ') 

121  F0RMAT(//7X, 'ENTER  MSTOP  (.LE.  13)’,//7X, 

1 ' ESTIMATED  "MINIMUM"  VALUE  I S : ' , 4X , 1 4 , //7X ) 

122  FQRMAT(//7X, 'ENTERING  NO.  OF  INTERNAL  AND  EXTERNAL  EXPANSION  MODES 

1 ' ,//7X. 

2 ' ESTIMATED  (XIXRMIN)  "MINIMUM"  INTERNAL  N1  =  '  ,2X, I4,//7X, 

3 ' ENTER  N1  ( . GE .  MSTOP);  ',//7X) 

123  FORMATC//7X, 'ESTIMATED  (KOxRMIN)  "MINIMUM"  EXTERNAL  N2  »',2X,I-., 

1//7X,  'ENTER  N2  (.GE.  MSTOP)-.  ',// 7X) 

124  FORMATC//7X, 'INDEXING  PROGRAM  THROUGH  VALUES  OF  M') 

125  FORMATC7X,  'M-LOOP  ....  M  =  M5) 

126  FORMAT(7X, 'CALL  MESH') 

127  FORMAT(7X,  'ENTER  I -LOOP ;  NO  STEPS  TO  COMPLETE:  M6) 

128  FORMAT (7X, 'I=',2X»I5,7X, 'NO  STEPS  TO  G0:',2X,I5) 

129  F0RMATC7X, 'CALL  LODER’) 

130  FORMAT(7X, 'CALL  MARCH') 

131  FORMAT (  7X, 'EX  I-LOOP,  SWEEP') 

132  FORMAT (  7X, 'CALL  MODECO') 

133  FORMAT (7X, 'CALL  RADFL D '  ) 

134  F0RMAT(7X, 'CALL  CRSSEC’) 

135  FORMAT ( 7X , 'EX  M-LOOP,  DATOUT') 

136  FORMAT(//7X, ' xxxxxxxxxxxxxxx  EMCAD  COMPLETED  xxxxxxxxxxxxxx ' ) 

137  F0RMATC//7X, ' ENTER  MATERIAL  PARAMETER  FILE  ( D: FILENAME. EXTENSION) 

1  ;  '  ) 

138  FORMAT (2(15)) 

201  FORMAT (//7X ,  '  RICCATI  MARCHING  STEPS') 

204  FORMAT C//7X , '  xxxxxxxxxxxxxx  M  =  ',15,'  xxxxxxxxxxxxxx '// ) 

300  F0RMAT(7X,I5,5E12.3) 

99  STOP 
END 


The  above  cons-ti tutes  only  the  main  program  of  EMCAD.  The 
remaining  subroutines  would  take  approximately  26  more  pages  to 
print.  The  entire  program  code  can  be  obtained  by  sending  an 
IBM  PC  compatible  formatted  5  1/4  "  floppy  disk  to 

Prof.  M.  A.  Morgan 
Code  62  MW 

Naval  Postgraduate  School 
Monterey,  Ca .  93943 
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APPENDIX  D 


EMCADIN  USER'S  MANUAL 

Welcome  to  EMCADIN.  This  program  is  designed  to  accept 
properly  formatted  input  data  from  the  West  Coast 
Consultants  commercial  shareware  product  CURVE  DIGITIZER 
(Appendix  M)  and  convert  it  to  a  form  which  can  be  used  by 
the  scattering  calculation  program  EMCAD.  CURVE  DIGITIZER 
is  a  computer  aided  design  (CAD)  package  which  will  allow 
user  manipulation  of  the  data. 

A.  BACKGROUND 

EMCADIN  was  written  in  the  programming  language  of  Ryan- 
McFarland  FORTRAN  (Appendix  M)  .  Any  questions  concerning 
formats,  changes  to  the  program,  etc...  should  be  referred 
to  the  respective  user's  manuals,  Dr.  M.  A.  Morgan  or  LT.  E. 
M.  Connolly  (Appendix  M) . 

B.  COMPUTER  SETUP 

EMCADIN  was  not  written  for  any  specific  computer 
configuration,  and  as  such,  should  run  on  most  IBM 
compatible  machines. 

C.  PROGRAM  INPUT  AND  OUTPUT 

* 

EMCADIN  obtains  its  input  from  a  previously  generated 
data  file.  The  input  data  file  for  EMCADIN  must  be  located 
in  the  same  directory/subdirectory  as  EMCADIN.  This  data 


file  is  produced  through  the  utilization  of  the  computer 
aided  design  package,  CURVE  DIGITIZER.  As  shown  in  Figure 
D-l ,  the  output  of  CURVE  DIGITIZER  is  the  input  to  EMCADIN. 
For  EMCADIN  and  CURVE  DIGITIZER  to  be  compatible,  it  was 
necessary  to  establish  and  adhere  to  the  predefined  data 
format  of  CURVE  DIGITIZER.  In  order  for  EMCADIN  to  be 
utilized  successfully,  it  is  imperative  that  the  user 
understand  the  need  for  the  particular  format  being  used,  as 
well  as  the  format  itself.  The  transition  from  CURVE 
DIGITIZER  to  EMCADIN  provides  more  of  an  opportunity  for 
catastrophic  error  than  at  any  other  point  in  this  entire 
package  of  interactive  computer  programs. 

Understanding  the  format  of  the  input  of  the  data  for 
EMCADIN  necessitates  a  thorough  understanding  of  some  of 
the  idiosyncracies  of  CURVE  DIGITIZER,  referred  to 
henceforth  for  brevity  as  CD.  Computer  aided  graphic  design 
deals  in  the  basic  business  of  handling  points.  Any  graphic 
figure  can  be  represented  as  some  combination  of  points. 
Lines  can  be  formed  by  connecting  two  points.  Arcs  can  be 
represented  by  placing  more  points  closer  together  and  by 
approximating  the  curvature  with  sufficiently  small  line 
segments  between  points.  The  required  input-output  format 
for  CD  is  utilized  to  insure  that  adequate  information  is 
present  to  direct  CD  how  to  connect  the  points.  A  similar 
process  may  be  that  of  the  simple  children's  game  connect- 
the-dots.  A  series  of  dots  exists  which  define  a  picture, 
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if  connected  properly.  As  more  and  more  dots  are  used, 
more  complex  pictures  with  better  resolution  can  be 
constructed.  The  fundamental  assumption  employed  by  CD  is 
that  each  point  should  be  connected  to  the  next  point  listed 
in  the  data  file  unless  some  code  is  present  to  indicate 
otherwise.  For  simplicity  and  unrestricted  accessibility, 
CD  input  and  output  data  are  represented  in  terms  of  XY 
rectangular  coordinates.  CD  flags  the  break  between  two 
data  points  by  what  shall  be  referred  to  as  a  separation 
code.  Since  two  coordinates,  X  and  Y,  are  required  to 
establish  a  data  point,  two  separation  codes,  one  for  the  X 
data  and  one  for  the  Y  data  are  required  to  separate  two 
distinct  data  points.  The  separation  code  employed  is  of 
the  form  999990.  The  last  two  digits,  90,  may  change  to 
reflect  scaling,  color,  or  any  of  a  number  of  other  things. 
It  was  not  possible  to  determine  the  specifics  from  the  user 
manual.  In  the  specific  case  of  EMCADIN  input  data,  999990 
999990  codes  are  required  to  distinguish  the  structural 
data  describing  one  layer  from  that  describing  another 
layer. 

CD  stores  the  data  in  the  exact  order  in  which  the  data 
points  were  drawn.  There  are  other  methods  of  generating 
the  simple  geometric  shapes  demonstrated  in  the  next  two 
examples,  however,  these  examples  were  specifically  chosen 
to  illustrate  potential  pitfalls  that  could  occur  if  extreme 
care  is  not  taken.  For  the  first  example,  the  simple  square 
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of  Figure  D-2  was  constructed  utilizing  CD.  The  base  of  the 
square  was  drawn  by  defining  point  number  1,  and  then  by 
drawing  a  line  to  point  number  2.  The  top  was  likewise 
established  by  defining  point  number  4,  and  then  by  drawing 
a  line  to  point  number  3.  The  two  sides  were  then  formed  by 
connecting  point  number  4  to  point  number  1  with  a  line,  and 
by  connecting  point  number  2  to  point  number  3  with  a  line. 
This  is  one  of  many  methods  of  connecting  the  four  data 
points  that  define  this  square.  Figure  D-3  lists  the  output 
data  file  which  contains  the  CD  defining  data  for  this 
square.  Notice  there  are  significantly  more  than  four  data 
points  describing  this  object  which  requires  only  four  data 
points.  Recall  from  earlier  discussion  that  EMCADIN 
utilizes  the  999990  codes  to  distinguish  structural  data 
describing  one  layer  from  that  describing  another  layer. 
Figure  D-2  depicts  only  one  layer,  yet  EMCADIN  using  the 
data  file  of  Figure  D-3  would  interpret  the  data  as  four 
distinct  layers.  The  999990  codes  of  Figure  D-3  are  present 
due  to  the  order  in  which  the  data  points  were  entered.  Had 
the  square  been  generated  by  defining  point  number  3, 
connecting  it  to  point  number  2,  then  to  1,  then  to  4, 
there  would  have  been  no  999990  codes  separating  the  data, 
and  the  data  file  of  Figure  D-3  would  consist  of  only  four 
pairs  (XY  coordinates)  of  data. 

The  second  and  third  major  errors  that  might  occur  can 
also  be  demonstrated  using  Figure  D-2.  The  sole 
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-2,-2 


2,-2 

999990,  999990 
-2,  2 
2,  2 

999990,  999990 
-2,  2 
-2,-2 

999990,  999990 
2,-2 
2,  2 
END, END 


-5.5 

5.5 

-4.5 

11 

8.5 

4 


FIG.  D-3.  EXAMPLE  OF  AN  OUTPUT 
DATA  FILE  FOR  A  CD  GENERATED 
SQUARE 
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justification  for  utilizing  CD  was  its  ability  to  generate 
data  to  be  used  ultimately  by  the  scattering  program  EMCAD. 
The  main  body  of  the  thesis  discusses  the  requirement  of 
EMCAD  for  axisymmetric  bodies  of  revolution.  Any  three 
dimensional  axisymmetric  body  of  revolution  can  be  defined 
by  knowing  the  axis  of  revolution,  and  then  by  defining  the 
body  in  planar  terms.  As  an  illustration,  a  given  sphere  is 
a  three  dimensional  body  of  revolution.  The  sphere  has 
axial  symmetry  and  can  be  described  by  a  semi-circle  in  any 
meridian  plane  intersecting  the  axis.  Rotating  this  planar 
semicircle  about  the  axis  would  produce  the  original 
sphere.  EMCAD  requires  this  meridian  representation  of  the 
three  dimensional  body  as  its  input.  For  simplicity's  sake, 
the  CD  representation  of  the  meridian  plane  will  be  referred 
to  as  the  XY  plane,  with  X  being  the  horizontal  axis  and  Y 
being  the  vertical  axis,  as  well  as  the  axis  of  revolution. 
This  coordinate  system  is  for  discussion  purposes  only,  and 
is  not  to  be  confused  with  the  three  dimensional  coordinate 
system  of  EMCAD  theory.  Viewing  the  body  of  Figure  D-2  as  a 
two  dimensional  body,  the  Y  axis  is  the  required  axis  of 
revolution.  Therefore  all  CD  planar  representations  for 
EMCAD  should  be  constructed  with  values  on  the  X  axis 
restricted  to  positive  values. 

The  third  problem  occurs  when  removal  of  data  points  is 
attempted  in  order  to  obtain  a  proper  planar  representation. 
EMCADIN  performs  a  valuable  service  to  the  user  of  improving 
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the  resolution  of  input  data.  It  suffices  to  say  that  due 
to  the  inherent  mathematics  involved,  EMCAD  requires  that 
two  data  points  be  defined,  both  on  the  Y  axis.  The  first 
data  point  defines  the  beginning  of  the  object,  at  X  =  0, 
and  the  last  data  point  defines  the  end  of  the  object,  also 
at  X  *  0.  In  the  example  of  the  square,  the  beginning  data 
point  would  be  (0,2),  and  the  ending  data  point  would  be 
(0,-2).  Without  these  beginning  and  ending  data  points  for 
each  layer  of  data,  the  EMCADIN  interpolation  will  "bomb", 
and  will  produce  erroneous  results,  if  any  results  at  all. 

For  the  second  example,  the  simple  circle  of  Figure  D-4 
was  chosen.  CD  allows  the  user  to  construct  a  circle  by 
simply  defining  the  coordinates  of  the  center  of  the  circle, 
and  by  specifying  the  desired  radius.  In  every  instance,  CD 
generates  the  data  points  by  starting  at  point  3  on  Figure 
D-4,  progressing  to  point  2,  to  point  5,  to  point  4,  and 
back  to  point  3  in  a  counter  clockwise  pattern.  A  sample  of 
the  output  data  file  for  this  circle  is  shown  in  Figure  D-5. 
This  data,  for  a  chosen  radius  of  2,  reflects  the  same 
counter  clockwise  pattern,  which  is  directly  opposite  to  the 
pattern  required  by  EMCADIN  input  data.  In  EMCADIN,  theta 
is  defined  in  a  clockwise  pattern,  starting  on  the  positive 
y  axis  and  progressing  by  equal  step  sizes  clockwise  to  the 
negative  y  axis.  By  using  this  method  of  circle 
construction,  the  user  must  find  someway  to  edit  and 
rearrange  the  data  before  it  can  be  used  by  CD.  The  other 
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2,  0 

1.997964,  .0902255 
1.99186,  .1802673 
1.981699,  .269942 
1.967504,  .3590671 
1.949302,  .447461 
1.927131,  .5349439 
1.901037,  .6213375 
1.871071,  .7064659 
1.837295,  .7901558 
1.799778,  .8722369 
1  758597,  .9525418 
1.713835,  1.030907 
1.665583,  1.107173 
1.61394,  1.181185 
1.55901,  1.252792 
1.500906,  1.321848 
1.439746,  1.388212 
1  375654,  1.451749 
1.308761,  1.512331 
1.239203,  1.569833 
1.167122,  1.624138 
1  092665,  1.675137 
1  015982,  1.722725 
.9372309,  1.766805 
.8585713,  1.807287 
.7741675,  1.844089 
.6901873,  1.877137 
.6048018,  1.906362 
.5181847,  1.931705 
.4305126,  1.953115 
.3419638,  1.970549 
.2527188,  1.983969 
.1629591,  1.99335 
7  286744E-02 ,  1.998672 
-1 . 737241E-02 ,  1.999925 
-.1075769,  1.997105 
-.1975623,  1.990218 


FIG.  D-5.  EXAMPLE  OF  A  PORTION 
OF  AN  OUTPUT  DATA  FILE  FOR  A  CD 
GENERATED  CIRCLE 
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problem  that  arises,  particularly  with  circles,  is  that  as 
the  radius  becomes  smaller  and  smaller,  fewer  and  fewer 
points  are  used  to  construct  the  circle,  causing  the  data 
points  to  become  less  and  less  accurate. 

For  the  case  of  multiple  layered  bodies,  EMCADIN 
expects  that  the  data  point  values  for  each  layer  have  been 
entered  in  the  input  data  file  starting  with  the  layer 
closest  to  the  center  of  the  object,  labeled  as  Layer  1  in 
Figure  D-6,  then  advancing  outward  from  the  center  until  all 
layers  have  been  accounted  for.  It  has  already  been  stated 
that  data  points  are  listed  in  the  order  in  which  they  are 
entered,  which  means  that  not  only  might  data  points  in  a 
given  layer  be  out  of  proper  order,  but  data  points  in 
differing  layers,  or  even  the  layers  themselves  might  be 
interspersed.  At  the  present  time,  there  appears  to  be  only 
two  methods  to  preclude  the  abundance  of  possible  errors. 
The  first  of  these  dictates  the  use  of  extreme  care  when 
entering  the  data  points.  The  values  must  be  entered  in 
order,  starting  at  0  =  0  degrees  and  progressing  clockwise 
to  e  =  180  degrees  for  the  layer  closest  to  the  center 
point,  and  then  entered  for  the  next  closest  layer,  and  then 
the  next  closest  layer  .  .  .  until  all  layers  have  been 
described.  Even  a  small  band  defined  by  a  different 
material  parameter  must  be  defined  for  all  9  from  0  degrees 
to  180  degrees.  If  this  process  is  followed  exactly,  there 
is  no  need  for  the  second  method.  The  second  method 
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EMCADIN  VIEW  OF  LAYERED  DATA 


involves  the  use  of  some  text  editor  to  rearrange  the  data 
as  necessary,  remove  unnecessary  999990  codes,  insure  data 
for  layers  as  well  as  layers  are  in  the  proper  order,  etc. . . 
The  prudent  method  would  include  extreme  caution  coupled 
with  a  data  review  utilizing  a  text  editor. 

The  EMCADIN  program  output  data  consists  of  the  number 
of  layers  of  the  body  being  considered,  the  minimum  and 
maximum  values  of  radius,  and  the  listing  of  radius  values 
in  a  clockwise  manner  for  9  =  0  to  9  =  180  degrees  for  each 
layer. 

D.  PROGRAM  EXECUTION 

Once  the  correct  disk  or  directory  has  been  selected, 
EMCADIN  can  be  invoked  by  typing 

EMCADIN 

at  the  DOS  prompt.  The  program  will  load  itself  and  the 
user  will  see  the  following  screens: 

*  Welcome  to  EMCADIN  (Figure  D-7) 

This  screen  presents  some  brief  background  information. 

*  Input  Filename  (Figure  D-8) 

This  screen  prompts  the  user  for  the  name  of  the  file 
containing  the  input  data  for  EMCADIN.  If  the  file  is 
obtained  from  CD,  the  file  must  have  an  extension  of  .dat, 
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UELCOHE  10  ENCAD IN 

THIS  PROGRAM  ACCEPTS  PROPERLV  FORMATTED  INPUT  DATA 
FROM  CUHUE  DIGITIZER  AND  CONVERTS  IT  TO  A  FORM 
WHICH  CAN  BE  USED  B¥  ENCAD. 


PLEASE  PRESS  ENTER  TO  CONTINUE. 


FIG.  D-7.  WELCOME  TO 
EMCADIN  MENU 
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INPUT  FILENAME 


INPUIFN  IS  THE  INPUT  FILENAME  OF  THE  FILE  CONTAINING 
THE  DATA  OF  INTEREST  WHICH  WILL  BE  INTERPOLATED  ON  AND 
CONUERTED  TO  A  FORM  WHICH  CAN  SUBSEQUENTLY  BE  USED  BY 
ENCAD. 


PLEASE  ENTER  THE  NAME  OF  THE  INPUT  FILE.  THE  EXTENSION 
OF  THE  FILENAME  MUST  BE  .DAT,  I.E.  FILENANE.DAT 


FIG.  D-8.  EMCADIN  PROMPT 
FOR  INPUT  FILENAME 
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and  the  extension  must  be  entered  by  the  user  as  part  of  the 
filename.  The  input  data  file  must  be  in  the  same 
directory/ subdirectory  as  EMCADIN. 

*  Output  Filename  (Figure  D-9) 

This  screen  prompts  the  user  for  the  name  of  the  file  that 
will  hold  the  output  data  from  EMCADIN.  The  extension  of 
this  file  is  left  to  user  discretion,  but  must  be  included 
as  part  of  the  output  filename. 

*  Resolution  (Figure  D-10) 

EMCADIN  uses  linear  interpolation  to  provide  the  capability 
to  generate  output  data  with  better  resolution  than  the 
input  data.  This  screen  prompts  the  user  for  the  desired 
output  data  resolution,  in  the  units  of  degrees. 
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OUTPUT  FILENAHE 


OUTPUTFN  IS  THE  OUTPUT  FILENAME  OF  THE  FILE  CONTAINING 
THE  DATA  WHICH  HAS  BEEN  INTERPOLATED  AND  CONVERTED  TO  A 
FORN  WHICH  WILL  BE  USED  BY  ENCAD. 


PLEASE  ENTER  THE  NAME  OF  THE  OUTPUT  FILE.  THE  EXTENSION 
OF  THE  FILENAME  MUST  BE  INCLUDED  I.E.  FILENANE.DAT 


FIG.  D-9.  EMCADIN  PROMPT 
FOR  OUTPUT  FILENAME 
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RESOLUTION 


DELTHE  IS  THE  USER  INPUT  VALUE  OF  THE  DESIRED  THETA 
RESOLUTION  IN  DEGREES. 

PLEASE  ENTER  THE  DESIRED  DELIA  THETA  VALUE  IN  DEGREES. 


FIG.  D- 1 □.  EMCADIN  PROMP 
FOR  OUTPUT  DATA  RESOLUTION 
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APPENDIX  E 


EMCADIN  SOURCE  CODE 


PROGRAM  EMCADIN 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

WELCOME  TO  EMCADIN  WRITTEN  BY  LT  E.  M.  CONNOLLY. 

EMCADIN  IS  A  PROGRAM  DESIGNED  TO  ACCEPT  PROPERLY  FORMATTED 
INPUT  DATA  FROM  CURVE  DIGITIZER  AND  TO  CONVERT  IT  TO  A 
FORM  WHICH  CAN  BE  USED  BY  THE  PROGRAM  EMCAD,  WRITTEN  BY 
DR.  M.  A.  MORGAN. 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
VARIABLE  DEFINITIONS 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


INPUTFN 

OUTPUTFN 

I,  J,  K, 
L,  N 
RADMIN 
RADMAX 
PI 

DELTHE 

DELTHETA 

THETA(N) 

DEGTORAD 

NUMPTSIN 

NUMSETS 

DUMMYI 

DUMMY2 

XARRAYCI) 

YARRAY ( I ) 

A  ( I ,  J  ) 

B(  I , J  ) 

NUMPTQUT 

THETAI 

THETAK 

XFOUND 

YFOUND 


THE  INPUT  FILENAME  OF  THE  FILE  CONTAINING 
THE  DATA  OF  INTEREST  WHICH  WILL  BE  INTERPOLATED 
ON  AND  CONVERTED  TO  A  FORM  WHICH  CAN 
SUBSEQUENTLY  BE  USED  BY  EMCAD 
THE  OUTPUT  FILENAME  OF  THE  FILE  CONTAINING 
THE  DATA  WHICH  HAS  BEEN  INTERPOLATED  AND 
CONVERTED  TO  A  FORM  WHICH  WILL  BE  USED  BY 
EMCAD 

INTEGER  COUNTERS 

THE  MINIMUM  RADIUS  VALUE 
THE  MAXIMUM  RADIUS  VALUE 
A  CONSTANT 

THE  USER  INPUT  VALUE  OF  THE  DESIRED  THETA 
RESOLUTION  IN  DEGREES 

THE  USER  INPUT  VALUE  OF  THE  DESIRED  THETA 
RESOLUTION  IN  RADIANS 

AN  ARRAY  HOLDING  THE  DESIRED  VALUES  OF  THETA 
AS  DETERMINED  BY  DELTHE 

A  CONVERSION  FACTOR  FROM  DEGREES  TO  RADIANS 
THE  NUMBER  OF  POINTS  HELD  IN  A  GIVEN  INPUT 
DATA  SET 

THE  HUMBER  OF  INPUT  DATA  SETS 

AN  INPUT  VARIABLE  ALLOWING  SCREEN  ADVANCING 

DURING  MENU  SELECTION 

A  CHARACTER  STRING  USED  TO  LOCATE  THE  END  OF 
EACH  INPUT  DATA  SET 

AN  ARRAY  CONTAINING  THE  VALUES  OF  X  FOR  ALL 
THE  INPUT  DATA  SETS 

AN  ARRAY  CONTAINING  THE  VALUES  OF  Y  FOR  ALL 
THE  INPUT  DATA  SETS 

AN  ARRAY  CONTAINING  THE  VALUES  OF  X  FOR  ALL 
THE  J  INPUT  DATA  SETS 

AN  ARRAY  CONTAINING  THE  VALUES  OF  Y  FOR  ALL 
THE  J  INPUT  DATA  SETS 

THE  NUMBER  OF  POINTS  OUT,  DETERMINED  BY  THE 
DESIRED  RESOLUTION 

THE  VALUE  OF  THETA  DETERMINED  BY  THE  I-TH  INPUT 
X  AND  Y  VALUES 

THE  VALUE  OF  THETA  DETERMINED  BY  THE  K-TH  INPUT 
X  AND  Y  VALUES 

THE  DESIRED  VALUE  OF  X  CORRESPONDING  TO  A  GIVEN 
VALUE  OF  THETA ( N ) 

THE  DESIRED  VALUE  OF  Y  CORRESPONDING  TO  A  GIVEN 
VALUE  OF  THETA ( N ) 
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RADIUS!  N,  J  )  AN  ARRAY  CONTAINING  THE  OUTPUT  CALCULATED 

VALUES  OF  RADIUS  FOR  ALL  THE  INPUT  DATA  SETS 
M  SLOPE  OF  THE  LINE  CONTAINING  THE  I-TH  AND  K-TH 

DATA  POINTS 

3B  Y-INTERCEPT  OF  THE  LINE  CONTAINING  THE  I-TH  AND 

K-TH  DATA  POINTS 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
VARIABLE  DECLARATIONS 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
DOUBLE  PRECISION  XINPUT! 2000  ) ,  YINPUT! 2000  ) ,  RADIUS! 2000,5) 

DOUBLE  PRECISION  DELTHE,  DELTHETA,  THETA(2000),  THETAI ,  THETAK 
DOUBLE  PRECISION  A(2000,5),  B(2000,5) 

DOUBLE  PRECISION  DEGTORAD,  INTER1,  INTER2,  RADMIN,  RADMAX 
DOUBLE  PRECISION  PI,  DUMMY,  XFOUND,  YFOUND 
DOUBLE  PRECISION  BB,  M 

INTEGER  I,  J,  K,  L,  N,  NUMPTSIN,  NUMPTOUT,  NUMSETS 
CHARACTERX3  DUMMY1 ,  DUMMY2 
CHARACTERX12  INPUTFN,  OUTPUTFN 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
INITIAL  VALUES 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
PI  =  3.14159 
J  =  0 
K  =  0 
L  =  0 

NUMPTSIN  =  0 
NUMSETS  =  1 
DEGTORAD  =  PI  /  180. 

RADMIN  =  10000 
RADMAX  =  -10000 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
INPUT  DATA  OF  INTEREST 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITE!  *,110) 

READ(X.IOO)  DUMMY 1 
WRITE!*, 120) 

READ!*, 100)  INPUTFN 
WRITE!*, 130) 

READ!*, 100)  OUTPUTFN 
WRITE!*, 140) 

READ!*,*)  DELTHE 

***XXX***X*XXX**XXX*X*XXXXXXXX*XX****X**X**X*XXX*XXX****X*****X*X* 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
DETERMINE  LOCATION  OH  END  OF  DA [ A 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
OPEN! 7 ,  FILE  =  INPUTFN) 

OPEN! 8 .FILE  =  OUTPUTFN) 

10  READ  (7,100)  DUMMY2 

NUMPTSIN  =  NUMPTSIN  +  1 

IF! DUMMY2  .NE.  'END')  GOTO  10 

CLOSE! 7 ) 

OPEN! 7 ,  TILE  =  INPUTFN) 

N  =  1 
J  =  1 

DO  20  1=  1, NUMPTSIN  -  1 

READ  ( 7 , x )  XINPUT(I),  YINPUT(I) 

IKABS(XINPUTd)  -  999990)  .GT.  20.)  THEN 
A!  ,1,  J  )  =  XINPUT ( I  ) 

B ( N , J )  =  Y INPUT ( I ) 

NUMSETS  =  J 
ELSE 
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J  =  J  +  1 
N  -  0 
ENDIF 
N  =  N  +  1 
20  CONTINUE 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
INTERMEDIATE  CALCULATIONS 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
NUMPTOUT  =  INTC1S0.  /  DELTHE)  +  I 
DELTHETA  =  DELTHE  *  DEGTORAD 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
INTERPOLATION  FOR  THETA  RESOLUTION 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
DO  30  J  =  1,  NUMSETS 
DO  40  N  =  1,  NUMPTOUT 

DETERMINE  THETA  WHERE  DATA  POINTS  X  AND  Y  ARE  TO  BE  FOUND 
THETACN)  =  C  N  -  1  )  *  DELTHETA 
DO  50  I  =  1,  NUMPTSIN  -  1 
K  =  1-1 

IFC  I  .EQ.  1)  BC  K, J  )  =  BCI,J) 

IFC ABSC BC I , J  )  )  .LT.  .001)  THEN 
THETAI  =  PI/2 

ELSE  IFC  BC I, J  )  .LT.  0)  THEN 
THETAI  =  PI  -  DATANCDABSCACI,J)/BCI,J))) 

ELSE 

THETAI  =  DATANCABSCACI, J)/B(I, J))) 

ENDIF 

DETERMINE  IF  K-TH  DATA  PAIR  FORMS  ANGLE  GREATER  THAN  90  DEGREES 
ASSUMES  CENTER  OF  COORDINATE  SYSTEM  AT  X  =  0  Y  =  0 
IFCABSCBCK, J) )  . LT .  .001)  THEN 
THETAK  =  PI/2 

ELSE  I FC  BC  K , J  )  .LT.  0)  THEN 
THETAK  =  PI  -  DATANCDABSCACK, J)/BCK, J)) ) 

ELSE 

THETAK  =  DATAN(ABS(ACK,J)/BCK,J))) 

ENDIF 

CHECK  TO  SEE  IF  K-TH  POINT  IS  THE  DESIRED  DATA  POINT 
IFC  DABSCTHETACN)  -  THETAK)  .LT.  .001)  THEN 
XFOUND  =  A  C  K , J  ) 

YFOUND  =  B C  K  ,  J  ) 

RADIUSCN,J  )  =  2*PI*DSQRTCXF0UND**2  +  YF0UND*X2 ) 

I FC  RADIUSC  N . J  )  .LT.  RADMIN)  RADMIN  =  RADIUSCN,J) 

IFCRADIUSCN, J)  .GT.  RADMAX)  RADMAX  =  RADIUSCN.J) 

GOTO  40 

CHECK  TO  SEE  IF  I-TH  POINT  IS  THE  DESIRED  DATA  POINT 
ELSE  IFC DABSCTHETACN)  -  THETAI)  .LT.  .001)  THEN 
XFO"ND  =  A C  I , J  ) 

YFOUND  =  BCI.J) 

RADIUSC  N . J  )  =  2*PI*DSQRTCXF0UND*X2  +  YF0UND**2) 

IFCRADIUSCN, J)  .LT.  RADMIN)  RADMIN  =  RADIUSCN,J) 

IFCRADIUSCN, J)  .GT.  RADMAX)  RADMAX  =  RADIUSCN,J) 

GOTO  40 

CHECK  TO  SEE  IF  THE  DESIRED  DATA  POINT  IS  BETWEEN  THE  K-TH 
DATA  POINT  AND  THE  I-TH  DATA  POINT 
ELSE  IFCTHETACN)  .GE.  THETAK  .AND. 

1  THETACN)  . L  E .  THETAI)  THEN 

C  CHECK  FOR  VERTICAL  LINE 

IFCDABSCACI, J)-ACK, J))  .LT.  .0001)  THEN 
XFOUND  =  A C I ,  J  ) 

YFOUND  =  XFOUND  /  DTANC  THETACN)  ) 

RADIUSC  N, J )  =  2*PI*DSQRTCXF0UND**2  +  YP0UND**2) 

IFCRADIUSCN, J)  .LT.  RADMIN)  RADMIN  =  RADIUSCN,J) 

IFCRADIUSCN, J)  .GT.  RADMAX)  RADMAX  =  RADIUSCN,J) 

GOTO  40 
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C  CHECK  FOR  HORIZONTAL  LINE 

ELSE  IFCDABSCBCI, J)-BCK, J))  .LT.  .0001)  THEN 
YFOUND  =  BCI,J) 

XFOUND  =  YFOUND  *  DTANC  THETA(N)  ) 

RADIUSCN.J)  =  2XPI*DSQRTCXF0UND**2  +  YF0UND*X2) 

I FC  RADIUS C  N , J  )  . LT  .  RADMIN)  RADMIN  =  RADIUSCN.J) 

I FC  RADIUSC  N , J  )  .GT.  RADMAX)  RADMAX  =  RADIUSCN.J) 

GOTO  40 

C  SLOPE  OF  LINE  CONTAINING  I-TH  AND  K-TH  DATA  POINTS 
ELSE 

M  =  (  BCI,J)-BCK,J)  )  /  (  ACI,J}-ACK, J)  ) 

C  Y-INTERCEPT  OF  LINE  CONTAINING  I-TH  AND  K-TH  DATA  POINTS 
BB  =  BCI.J)  -CM*  AC  I ,  J  )  ) 

RECALL  EQUATION  OF  ANY  LINE  IS  Y  =  MX  +  B 

THE  LINE  CONTAINING  THE  I-TH  AND  K-TH  DATA  POINTS  INTERSECTS 
THE  LINE  DRAWN  FROM  THE  ORIGIN  THROUGH  THE  DESIRED  DATA  POINT 
AT  THE  DESIRED  DATA  POINT  XFOUND  AND  YFOUND.  THE  EQUATION  FOR 
THE  LINE  CONTAINING  THE  I-TH  AND  K-TH  DATA  POINTS  MUST  THEN 
ALSO  WORK  AT  THE  DESIRED  DATA  POINT.  WRITING  XFOUND  AND 
YFOUND  IN  TERMS  OF  THEIR  ANGLE  THETACN),  AND  SUBSTITUTING 
YFOUND  =  BB  /  C  1  -  C  M  *  DTANC  THETACN)  )  )  ) 

XFOUND  =  YFOUND  *  DTANC  THETACN)  ) 

RADI  US  C  N , J )  =  2*PI*DSQRTCXFQUND**2  +  YF0UND**2) 

IFCRADIUSCN, J)  .LT.  RADMIN)  RADMIN  =  RADIUSCN,J) 

IFCRADIUSCN, J  )  . GT  .  RADMAX)  RADMAX  =  RADIUSCN.J) 

GOTO  40 
ENDIF 

NEED  TO  TRY  A  NEW  I-TH  AND  K-TH  DATA  PAIR 
ELSE 

DUMMY  =  0 
ENDIF 

50  CONTINUE 
40  CGNTINUE 
30  CONTINUE 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

OUTPUT  TO  DATA  FILE 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITEC8, 103)  NUMPTOUT,  NUMSETS  +  1 
WRITEC  8 , 102)  RADMIN,  RADMAX 
DO  400  N  =  1, NUMPTOUT 

WRITEC 8 , 101 )  N, CRADIUSCN, J),  J=l, NUMSETS) 

400  CONTINUE 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
CLOSEC7) 

CL0SEC8 ) 

100  FORMAT C  A ) 

101  FORMAT  (I5.5CE14.6)) 

102  FORMAT C  2C  E14 . 6 ) ) 

103  FORMAT C  2C 15 ) ) 

110  FORMATC////////////////////5X, 'WELCOME  TO  EMCADIN' ,//5X, 

1 ' THIS  PROGRAM  ACCEPTS  PROPERLY  FORMATTED  INPUT  DATA',/5X, 

2 ' FROM  CURVE  DIGITIZER  AND  CONVERTS  IT  TO  A  FORM. ' */5X, 

3 ' WHICH  CAN  BE  USED  BY  EMCAD.  ',/5X, 

4////5X, 'PLEASE  PRESS  ENTER  TO  CONTINUE .' ,//5X) 

120  F0RMATC////////////////////5X, ' INPUT  FILENAME' ,//////5X, 

1 ’ INPUTFN  IS  THE  INPUT  FILENAME  OF  THE  FILE  CONTAINING' ,/5X, 

2 ' THE  DATA  OF  INTEREST  WHICH  WILL  BE  INTERPOLATED  ON  AND',/5X, 

3 'CONVERTED  TO  A  FORM  WHICH  CAN  SUBSEQUENTLY  BE  USED  BY',/5X, 

4 ' EMCAD . ’ ,///5X, 

5 ' PLEASE  ENTER  THE  NAME  OF  THE  INPUT  FILE.  THE  EXTENSION' ,/5X, 

6 ' OF  THE  FILENAME  MUST  BE  .DAT,  I.E.  FI L ENAME . DAT ' , //5X ) 

130  F0RMATC////////////////////5X, ’OUTPUT  FILENAME' ,//////5X, 

1 ' OUTPUTFN  IS  THE  OUTPUT  FILENAME  OF  THE  FILE  CONTAINING’ ,/5X, 

2 'THE  DATA  WHICH  HAS  BEEN  INTERPOLATED  AND  CONVERTED  TO  A’,/5X, 
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3 ' FORM  WHICH  WILL  BE  USED  BY  EMCAD. ' ,///5X, 

A 'PLEASE  ENTER  THE  NAME  OF  THE  OUTPUT  FILE.  THE  EXTENSION' ,/5X, 
6'OF  THE  FILENAME  MUST  BE  INCLUDED  I.E.  FI L ENAME . DAT ' , //5X ) 

140  FORMAT (//////////////////// 5X , 'RESOLUTION', ////// 5X , 

1 ' DEL  THE  IS  THE  USER  INPUT  VALUE  OF  THE  DESIRED  THETA', /5X, 
'RESOLUTION  IN  DEGREES ///5X , 

'PLEASE  ENTER  THE  DESIRED  DELTA  THETA  VALUE  IN  DEGREES .', //SX ) 
END 


APPENDIX  F 


CRVDIGIN  USER'S  MANUAL 

Welcome  to  CRVDIGIN.  This  program  is  designed  to 
accept  properly  formatted  input  data,  and  to  convert  it  to  a 
form  which  can  be  used  by  the  West  Coast  Consultants 
commercial  shareware  product  CURVE  DIGITIZER  (Appendix  M)  . 
CURVE  DIGITIZER  is  a  computer  aided  design  (CAD)  package 
which  will  allow  user  manipulation  of  the  data.  Upon 
completion  of  using  CURVE  DIGITIZER,  the  user  may  employ  the 
program  EMCADIN  to  transform  the  data  into  a  form  which  can 
be  used  by  the  scattering  calculation  program  EMCAD. 

A.  BACKGROUND 

EMSCAT  was  written  by  Dr,  M.  A.  Morgan  while  at  the 
University  of  California  at  Berkeley.  Extensive  interactive 
computer  aided  design  type  feature  have  been  added  to 
EMSCAT,  which  consider  scattering  from  inhomogeneous 
axisymmetric  objects.  These  computer  aided  design  features, 
driving  the  change  of  program  name  from  EMSCAT  to  EMCAD, 
allow  structures  composed  of  multiple  dielectric  materials 
to  be  graphically  constructed  or  modified  by  the  user,  and 
then  evaluated  by  EMCAD  for  the  scattering.  CRVDIGIN  was 
written  in  the  programming  language  of  Ryan-McFarland 
FORTRAN  (Appendix  M)  .  Any  questions  concerning  formats, 
changes  to  the  program,  etc.  should  be  referred  to  the 
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respective  user's  manuals,  Dr.  M.  A.  Morgan  or  LT.  2.  M. 
Connolly  (Appendix  M) . 

B.  COMPUTER  SETUP 

CRVDIGIN  was  not  written  for  any  specific  computer 
configuration,  and  as  such,  should  run  on  most  IBM 
compatible  machines. 

C.  PROGRAM  INPUT  AND  OUTPUT 

CRVDIGIN  obtains  its  input  from  a  previously  generated 
data  file.  The  input  data  file  for  CRVDIGIN  must  be  located 
in  the  same  directory/subdirectory  as  CRVDIGIN.  This  data 
file  may  be  produced  by  one  of  two  methods.  Tne  first 
method  consists  of  special  case  computer  programs  which 
generate  data  files  describing  the  penetrable  body  of 
revolution  for  the  specific  geometrically  shaped  object 
being  considered.  A  few  examples  of  these  special  case 
programs  are  SPHERE,  CONE  and  CYLINDER  listed  in  Appendices 
J,  K  and  L  respectively.  Each  of  these  programs  generates 


two  descriptive 

data 

files 

which 

can  be  separated  into 

distinct  areas 

of 

information. 

The  first 

area  of 

information  is 

the 

data 

which 

conveys  the 

material 

composition  of  each  layer  of  the  structure.  This  data  is  of 
no  concern  at  this  particular  time,  however,  it  will  be 
required  before  running  EMCAD.  The  second  area  of 
information,  which  is  the  data  required  as  input  for 
CRVDIGIN,  is  the  data  which  defines  the  structural  shape  of 
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each  layer  of  the  penetrable  body.  This  data  is  referred  to 
as  "STRUCTUR" .  DAT,  and  is  shown  in  the  flow  diagram  of 
Figure  F-l.  The  second  method  of  input  data  file 
construction  consists  of  user  generation  through  direct 
keyboard  input. 

As  shown  in  Figure  F-l,  the  output  of  CRVDIGIN,  shown 
as  "CURVEIN" .  DAT,  is  the  input  to  the  CAD  package  CURVE 
DIGITIZER.  For  CRVDIGIN  and  CURVE  DIGITIZER  to  be  able  to 
accept  input  data  from  a  variety  of  sources,  it  was 
necessary  to  establish  and  adhere  to  a  predefined  data 
format.  The  special  case  computer  programs  were  written  to 
incorporate  this  format.  In  order  for  the  second  method  to 
be  utilized  successfully,  as  well  as  for  other  special  case 
programs  to  be  written,  it  is  imperative  that  the  user 
understand  the  need  for  the  particular  format  being  used,  as 
well  as  the  format  itself.  The  input  format  for  CRVDIGIN  is 
relatively  uncomplicated.  An  example  of  a  portion  of  an 
input  data  file  for  CRVDIGIN  is  shown  in  Figure  F-2 .  The 
first  line  of  data  contains  two  numbers  in  integer  format. 
The  first  integer  represents  the  number  of  data  points  used 
to  describe  any  layer.  The  second  integer  lists  the  total 
number  of  layers  +  1.  The  '+  1'  was  used  to  correct  some 
minor  discrepancies  during  the  design  of  the  code.  The  data 
defining  the  physical  structure  of  the  boundary  for  each 
layer  is  written  in  radius  versus  theta  format.  The  value 
of  the  angle  theta  associated  with  each  radius  is  determined 
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SCATTERING 


CRVD1GIN  AND 
FLOW  DIAGRAM 


361  4 

0 . 250300E+01  0 . 382760E+02 

1  0 . 188496E+02  0.251327E+02  0.314159E+02 

2  0 . 176975E+02  0.241175E+02  0.304952B+02 

3  0 . 166793E+02  0.231827E+02  0 . 296292E+02 

4  0 . 157731E+02  0.223194E+02  0.288130E+02 

5  0 . 149613E+02  0 . 215196E+02  0.280428E+02 

6  0 . 142301E+02  0 , 207767E+02  0.273146E+02 

7  0 . 135680E+02  0 . 200848E+02  0.266253E+02 

8  0 . 129657E+02  0.194390E+02  0.259719E+02 

9  0. 124155E+02  0.188348E+02  0.253516E+02 

10  0 . 119109E+02  0 . 182884E+02  0.247621E+02 

11  0 . 1 14466E+02  0 . 177363E+02  0.242013E+02 

12  0 . 110180E+02  0 . 172357E+02  0 . 236670E+02 

13  0 . 106211B+02  0 . 167638E+02  0.231575E+02 

14  0 . 102525E+02  0.183182E+02  0.228712E+02 

15  0 . 990939E+01  0.158969E+02  0.222066E+02 

16  0 . 958921E+01  0.154980E+02  0.217622E+02 

17  0 . 928975E+01  0.151197E+02  0.213369B+02 

18  0 . 900910E+01  0 . 147605E+02  0.209294E+02 

19  0 . 874555E+01  0.144191E+02  0.205387E+02 

20  0 . 849761E+01  0.140941E+02  0.201838E+02 

21  0 . 826396E+01  0.137845E+02  0.198040E+02 

22  0. 80434 1E+01  0.134893E+02  0.194582E+02 

23  0.783490E+01  0.132074E+02  0.191256E+02 

24  0 . 763750E+01  0.129380E+02  0.188057E+02 

25  0 . 745035E+01  0.128803E+02  0 . 184978E-H02 

26  0 . 727270E+01  0.124336E+02  0.182009E+02 

27  0 . 710385E+01  0.121972E+02  0.179149B+02 

28  0 . 694318E+01  0.119706E+02  0.176390E+02 

29  0 . 679012E+01  0.117531E+02  0.173728E+02 

30  0 . 664416E+01  0.115442E+02  0.171158E+02 


FIG.  F-2.  PORTION  OF  SAMPLE 
INPUT  DATA  FOR  CRVDIGIN.EXE 
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by  starting  at  theta  equal  to  0  degrees  and  progressing  to 
theta  equals  180  degrees  in  step  sizes  determined  by  the 
number  of  points  available.  The  program  expects  that  the 
values  of  the  radius  for  each  layer  have  been  entered 
starting  with  the  layer  closest  to  the  center  of  the 
object,  then  advancing  outward  from  the  center  until  all 
layers  have  been  accounted  for.  The  second  line  of  data 
contains  two  real  values.  The  first  real  value  is  the 
minimum  radius  value  present,  and  the  second  value  is  the 
maximum  radius  value  present.  Following  the  description  for 
layer  entry  above,  this  implies  that  the  minimum  radius 
value  should  come  from  the  first  layer,  and  the  maximum 
radius  value  should  come  from  the  last  layer.  Figure  F-3 
demonstrates  the  indexing  of  layers.  The  next  A  lines  of 
data  consists  of  B  columns  of  data.  The  value  of  A  is  the 
number  of  points,  and  the  value  of  B  is  the  number  of  layers 
+  1.  The  first  column  value  is  an  index  value  starting  at  1 
and  incrementing  by  1.  The  second,  third,  ...  columns 
contain  values  of  the  radius  associated  with  theta  =  (index 
-  1)  *  deltatheta  for  the  first  layer,  the  second  layer  . . . 
For  each  successive  row  the  index  changes,  the  value  of 
theta  changes  and  the  radii  associated  with  that  given  value 
of  theta  are  listed. 

The  output  of  CRVDIGIN,  used  as  an  input  to  CURVE 
DIGITIZER,  is  just  as  simple  in  appearance  as  the  input  to 
CRVDIGIN,  yet  much  more  difficult  in  theory.  The 


133 


MULTI-LAYERED.  MULTI-HEIGHT 


FIG.  F-3.  GRAPHICAL  REPRESENTATION  OF 
ASSIGNMENT  OF  LAYER  BOUNDARIES 


explanation  of  this  difficulty  actually  is  embedded  in  the 
operation  of  CURVE  DIGITIZER.  Although  these  comments  may 
apply  to  other  CAD  packages,  they  are  specifically  directed 
towards  CURVE  DIGITIZER,  referred  to  henceforth  for  brevity 
as  CD.  Computer  aided  graphic  design  deals  in  the  basic 
business  of  handling  points.  Any  graphic  figure  can  be 
represented  as  some  combination  of  points.  Lines  can  be 
formed  by  connecting  two  points.  Arcs  can  be  represented  by 
placing  more  points  closer  together  and  by  approximating  the 
curvature  with  sufficiently  small  line  segments  between 
points.  The  required  input-output  format  for  CD  is 
utilized  to  insure  that  adequate  information  is  present  to 
direct  CD  how  to  connect  the  points.  A  similar  process  may 
be  that  of  the  simple  children's  game  connect-the-dots.  A 
series  of  dots  exists  which  define  a  picture,  if  connected 
properly.  As  more  and  more  dots  are  used,  more  complex 
pictures  with  better  resolution  can  be  constructed.  The 
fundamental  assumption  employed  by  CD  is  that  each  point 
should  be  connected  to  the  next  point  listed  in  the  data 
file  unless  some  code  is  present  to  indicate  otherwise.  For 
simplicity  and  unrestricted  accessibility,  CD  input  and 
output  data  are  represented  in  terms  of  XY  rectangular 
coordinates.  CD  flags  the  break  between  two  data  points  by 
what  shall  be  referred  to  as  a  separation  code.  Since  two 
coordinates,  X  and  Y,  are  required  to  establish  a  data 
point,  two  separation  codes,  one  for  the  X  data  and  one  for 


the  Y  data  are  required  to  separate  two  distinct  data 
points.  The  separation  code  employed  is  of  the  form  999990. 
The  last  two  digits  may  change  to  reflect  scaling,  color,  or 
any  of  a  number  of  other  things.  It  was  not  possible  to 
determine  the  specifics  from  the  user  manual.  In  the 
specific  case  of  CRVDIGIN  output  data,  999990  999990  codes 

are  used  to  distinguish  the  structural  data  describing  one 
layer  from  that  describing  another  layer.  Other  than  the 
insertion  of  separation  codes  at  the  appropriate  locations, 
CRVDIGIN  performs  the  simple  mathematical  task  of  converting 
the  input  data  generated  in  radius-theta  format  to  output 
data  represented  in  XY  format. 

D.  PROGRAM  EXECUTION 

Once  the  correct  disk  or  directory  has  been  selected, 
CRVDIGIN  can  be  invoked  by  typing 

CRVDIGIN 

at  the  DOS  prompt.  The  program  will  load  itself  and  the 
user  will  see  the  following  screens: 

*  Welcome  to  CRVDIGIN  (Figure  F-4) 

This  screen  presents  some  brief  background  information. 
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UELCOHE  70  CFUDIGIN 


THIS  PROGRAM  ACCEPTS  PROPERLV  FORMATTED  INPUT  DATA 
AND  CONUERTS  IT  TO  A  FORM  WHICH  CAN  BE  USED  BY  THE 
COMMERCIAL  SHAREWARE  PRODUCT  CURUE  DIGITIZER. 

CURVE  DIGITIZER  IS  A  CAD  PACKAGE  WHICH  WILL  ALLOW  USER 
MANIPULATION  OF  THE  DATA.  UPON  COMPLETION  OF 
CURVE  DIGITIZER,  THE  USER  MAY  UTILIZE  THE  PROGRAM 
EMCADIN  TO  TRANSFORM  THE  DATA  INTO  A  FORM  WHICH  CAN 
BE  USED  BY  THE  PROGRAM  ENCAD  WRITTEN  BY  DR.  N.  A.  MORGAN. 


PLEASE  PRESS  ENTER  TO  CONTINUE. 


FIG.  F-4.  WELCOME  TO 
CRVDIGIN  MENU 
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*  Input  Filename  (Figure  F-5) 

This  screen  prompts  the  user  for  the  name  of  the  file 
containing  the  input  data  for  CRVDIGIN.  The  file  may  have 
any  extension  that  the  user  desires,  however,  the  extension 
must  be  entered  by  the  user  as  part  of  the  filename.  The 
input  data  file  must  be  in  the  same  directory  as  CRVDIGIN. 

*  Output  Filename  (Figure  F-6) 

This  screen  prompts  the  user  for  the  name  of  the  file  that 
will  hold  the  output  data  from  CRVDIGIN.  Because  of  the 
input  requirements  of  CURVE  DIGITIZER,  the  extension  of  this 
file  must  be  .dat,  and  must  be  included  as  part  of  the 
output  filename. 

*  Resolution  (Figure  F-7) 

CRVDIGIN  uses  linear  interpolation  to  provide  the  capability 
to  generate  output  data  with  better  resolution  than  the 
input  data.  This  screen  prompts  the  user  for  the  desired 
output  data  resolution,  in  the  units  of  degrees. 

E .  PROGRAM  OUTPUT 

A  portion  of  an  example  of  CRVDIGIN  output  is  shown  in 
Figure  F-8.  Notice  the  999990  999990  codes  separating  the 
different  layers  of  data. 
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INPUT  FILENAME 


INPUIFN  IS  THE  INPUT  FILENAME  OF  THE  FILE  CONTAINING 
THE  DATA  OF  INTEREST  WHICH  WILL  BE  INTERPOLATED  ON  AND 
COHUERTED  TO  A  FORM  UHICH  CAN  SUBSEQUENTLY  BE  USED  BY 
THE  CURUE  DIGITIZER  PROGRAM. 


PLEASE  ENTER  THE  NAME  OF  THE  INPUT  FILE.  THE  EXTENSION 
OF  THE  FILENANE  MUST  BE  INCLUDED,  I.E.  FILENANE.DAT 


FIG.  F-5.  CRVDIGIN  PROMP 
FGR  INPUT  FILENAME 
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OUTPUT  FILENAME 


OUTPUTFN  IS  THE  OUTPUT  FILENAME  OF  THE  FILE  CONTAINING 
THE  DATA  UHICH  HAS  BEEN  INTERPOLATED  AND  CONVERTED  TO  A 
FORM  UHICH  WILL  BE  USED  BV  THE  CURVE  DIGITIZER  PROGRAM. 


PLEASE  ENTER  THE  NAME  OF  THE  OUTPUT  FILE.  THE  EXTENSION 
OF  THE  FILENAME  MUST  BE  .DAT,  I.E.  FILENAHE.DAT 


FIG.  F-6. 
FOR  OUTPUT 


CRVDIGIN  PROMPT 
FILENAME 
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RESOLUTION 

CELTHE  IS  THE  USER  INPUT  VALUE  OF  THE  DESIRED  THETA 
RESOLUTION  IN  DEGREES. 

PLEASE  ENTER  THE  DESIRED  DEIA  THETA  VALUE  IN  DEGREES. 


FIG.  F-7.  CRVDIGIN  PROMP 
FOR  OUTPUT  DATA  RESOLUTION 
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4.69975948 
4.35180850 
4.00664711 
3.66403151 
3.32372093 
2.98551798 
2.64917397 
2.31447673 
1.98121083 
1 . 64916909 
1 . 31813526 
0.98790431 
0.65828347 
0.32906568 

4 . 32930647E-05 
999990 

0 . 00000000E-01 
0.40459359 
0  75102264 
1.05115640 
1.31384826 
1.54581928 
1.75229287 
1.93734324 
2.10425401 
2.25564241 
2.39365172 


-18.84953690 
-18.84956930 
-18.84958080 
-18.84959030 
-18.84950830 
-18.84953120 
-18.84954830 
-18.84953880 
-18 . 84956740 
-18.84959410 
-18.84956740 
-18.84953120 
-18.84950830 
-18.84952740 
-18.84959980 
999990 

25.13270000 

23.17916870 

21.50649070 

20.05727580 

18.78891940 

17.66880800 

16.67196660 

15.77840710 

14.97255710 

14.24157710 

13.57508280 


FIG.  F-8.  PORT  I CN  OF  SAMPLE 
OUTPUT  DATA  FROM  CRVDIGIN 
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APPENDIX  G 


CSV D I  GIN  SOURCE  CODE 


PROGRAM  CRVDIGIN 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

WELCOME  TO  CRVDIGIN  WRITTEN  BY  LT  E.  M.  CONNOLLY. 

CRVDIGIN  IS  A  PROGRAM  DESIGNED  TO  ACCEPT  PROPERLY  FORMATTED 
INPUT  DATA,  AND  CONVERT  IT  TO  A  FORM  WHICH  CAN  BE  USED  BY 
THE  COMMERCIAL  SHAREWARE  PRODUCT  CURVE  DIGITIZER. 

CURVE  DIGITIZER  IS  A  CAD  PACKAGE  WHICH  WILL  ALLOW  USER 
MANIPULATION  OF  THE  DATA.  UPON  COMPLETION  OF  CURVE  DIGITIZER, 

THE  USER  MAY  UTILIZE  THE  PROGRAM  EMCADIN,  ALSO  WRITTEN  BY 
LT  E.  M.  CONNOLLY  TO  TRANSFORM  THE  DATA  INTO  A  FORM  WHICH  CAN 
BE  USED  BY  THE  PROGRAM  EMCAD  WRITTEN  BY  DR.  M.  A.  MORGAN. 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

VARIABLE  DEFINITIONS 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


INPUTFN 


OUTPUTFN 

I,  J,  K, 
L,  N 

NUMBROWS 

NUMBCOLS 

SEPRATN 

RADARRAY 

MINRAD 

MAXRAD 

RADIUS 

PI 

DELTHE 

DELTHETA 

THETAOLD 

THETA 

DEGTORAD 

X 


THE  INPUT  FILENAME  OF  THE  FILE  CONTAINING 
THE  DATA  OF  INTEREST  WHICH  WILL  BE  INTERPOLATED 
ON  AND  CONVERTED  TO  A  FORM  WHICH  CAN 
SUBSEQUENTLY  BE  USED  BY  THE  CURVE  DIGITIZER 
PROGRAM 

THE  OUTPUT  FILENAME  OF  THE  FILE  CONTAINING 
THE  DATA  WHICH  HAS  BEEN  INTERPOLATED  AND 
CONVERTED  TO  A  FORM  WHICH  WILL  BE  USED  BY 
THE  CURVE  DIGITIZER  PROGRAM 
INTEGER  COUNTERS 

THE  NUMBER  OF  ROWS  CONTAINED  IN  THE  DATA 
MATRIX  IN  THE  FILE  INPUTFN 
THE  NUMBER  OF  COLUMNS  CONTAINED  IN  THE  DATA 
MATRIX  IN  THE  FILE  INPUTFN 

A  FLAG  USED  BY  CURVE  DIGITIZER  TO  DISTINGUISH 
BETWEEN  DATA  POINTS  OR  SETS  OF  DATA 
A  WORKING  ARRAY  USED  TO  HOLD  ALL  THE  VALUES 
OF  THE  ARRAY  SPECIFIED  ABOVE  BY  NUMBROWS  AND 
NUMBCOLS 

THE  MINIMUM  RADIUS  VALUE,  NOT  USED  BY  THIS 
PROGRAM  BUT  REQUIRED  FOR  EMCAD 
THE  MAXIMUM  RADIUS  VALUE,  NOT  USED  BY  THIS 
PROGRAM  BUT  REQUIRED  FOR  EMCAD 
A  WORKING  VARIABLE  USED  TO  HOLD  THE  CURRENT 
VALUE  OF  THE  RADIUS  OBTAINED  FROM  RADARRAY 
A  CONSTANT 

THE  USER  INPUT  VALUE  OF  THE  DESIRED  THETA 
RESOLUTION  IN  DEGREES 

THE  USER  INPUT  VALUE  OF  THE  DESIRED  THETA 
RESOLUTION  IN  RADIANS 

A  WORKING  VARIABLE  USED  TO  HOLD  THE  CURRENT 
VALUE  OF  THETA  OBTAINED  FROM  RADARRAY 
AN  ARRAY  HOLDING  THE  DESIRED  VALUES  OF  THETA 
AS  DETERMINED  BY  DELTHE 

A  CONVERSION  FACTOR  FROM  DEGREES  TO  RADIANS 
THE  VALUE  OF  X  OBTAINED  BY  CONVERTING  THE  POLAR 
COORDINATES  RADIUS, THETA  TO  X,Y  COORDINATES 
THE  X , Y  COORDINATE  SYSTEM  IS  REQUIRED  FOR  THE 
DATA  TO  BE  COMPATIBLE  WITH  THE  CURVE  DIGITIZER 
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oooo  oooo  oooo  oooo  oooo  ooooooooooo 


PROGRAM 

Y  THE  VALUE  OF  Y  OBTAINED  BY  CONVERTING  THE  POLAR 

COORDINATES  RADIUS, THETA  TO  X,Y  COORDINATES 
THE  X , Y  COORDINATE  CVSTEM  IS  REQUIRED  FOR  THE 
DATA  TO  BE  COMPATIBLE  WITH  THE  CURVE  DIGITIZER 
PROGRAM 

INTER1,  INTER2  INTERMEDIATE  CALCULATION  VARIABLES 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
VARIABLE  DECLARATION 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
CHARACTERX12  INPUTFN ,  OUTPUTFN,  DUMMY 
INTEGER  I,  J,  K,  L,  N 

INTEGER  NUMBROWS,  NUMBCOLS,  NUMPTS,  SEPRATN 
REAL  RADARRAYC1000,10) 

REAL  MINRAD,  MAXRAD,  RADIUS,  PI,  DELTHE,  DELTHETA,  THETAOLD 
REAL  THETAC1000),  DEGTORAD,  X,  Y 
REAL  INTER1 ,  INTER2 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
INITIAL  VALUES 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
PI  =  3.14159 
DEGTORAD  =  PI  /  180. 

SEPRATN  =  999990 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
INPUT  DATA  OF  INTEREST 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRI TEC  x , 1 1 0  ) 

READC  x, 100 )  DUMMY 
WRITEC  x, 120  ) 

READC  x , 1 00 )  INPUTFN 
WRITECX, 130) 

READC  x, 100 )  OUTPUTFN 
WRITECX, 140) 

READC  x , X )  DELTHE 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
READ  DATA  FILE 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
0PENC7 , FILE  =  INPUTFN) 

OPEN  C  3 , FI L  E  =  OUTPUTFN) 

READC  7 , X  )  NUMBROWS,  NUMBCOLS 
READC  7 , x )  MINRAD,  MAXRAD 

READC  7 , X )  C  C  RADARRAY  C I , J ) ,  J  =  1 , NUMBCOLS) ,  I  =  1, NUMBROWS) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
INTERMEDIATE  CALCULATIONS 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

NUMPTS  =  INTC180.  /  DELTHE)  +  1 

DELTAOLD  =  C180.  /  CNUMBROWS  -  1))  x  DEGTORAD 

DELTHETA  =  DELTHE  x  DEGTORAD 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
INTERPOLATION  AND  TRANSFORMATION 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
J  =  1 

DO  10  K  =  2,  NUMBCOLS 
DO  20  M  »  1,  NUMPTS 
THETACN)  =  CN-1)  X  DELTHETA 
DO  30  I  =  1,  NUMBROWS 

THETAOLD  =  C  RADARRAY  C I , J  )  -  1)  x  DELTAOLD 
I FC  ABSC  THETAC  N )  -  THETAOLD)  . L  T  .  .00001  )  THEN 
RADIUS  =  RADARRAY  C I , K ) 

X  =  RADIUS  x  SINC  THETACN)  ) 
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Y  -  RADIUS  *  C0S(  THETA(N)  ) 
l-.'RI TEC  8 , X )  X,Y 

GOTO  20 

ELSE  IF  (THETA(N)  .GE.  (THETAOLD  -  DELTAOLD)  .AND. 

1  THETACN)  . L  E .  THETAOLD)  THEN 

L  -  I  -  1 

INTER1  =  RADARRAY  (  L  ,  K  )  -  RADARRAY ( I ,  K  ) 

INTER2  -  (THETA(N)  -  THETAOLD)  /  (-DELTAOLD) 

RADIUS  =  INTER1  *  INTER2  +  RADARRAY ( I , K ) 

REMOVES  NORMALIZING  FACTOR  OF  2*PI  GENERATED  BY  SPHERE,  CONE, 
ETC..,  FACTOR  IS  RESTORED  IN  EMCADIH 
X  =  RADIUS  *  SI N (  THETA ( N )  )  *  1  /  (  2  *  PI  ) 

Y  =  RADIUS  *  COS (  THETA(N)  )  X  1  /  (  2  *  PI  ) 

WRITEC 8 , x )  X , Y 

GOTO  20 
ELSE 

DUMMY  BRANCH,  TRY  NEXT  VALUE  OF  THETAOLD 

RADIUS  =  0 

ENDIF 

30  CONTINUE 
20  CONTINUE 

WRITEC S, X)  SEPRATN, SEPRATN 
10  CONTINUE 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
100  FORMAT ( A ) 

110  FORMAT {///////////////// f / /SI,  'WELCOME  TO  CRVDIGIN ' ,//5X, 

1 ' THIS  PROGRAM  ACCEPTS  PROPERLY  FORMATTED  INPUT  DATA',/5X, 

2 'AND  CONVERTS  IT  TO  A  FORM  WHICH  CAN  BE  USED  BY  THE  ',/5X, 

3 'COMMERCIAL  SHAREWARE  PRODUCT  CURVE  DIGITIZER .' ,/5X, 

4'CURVE  DIGITIZER  IS  A  CAD  PACKAGE  WHICH  WILL  ALLOW  USER',/5X, 

5 'MANIPULATION  OF  THE  DATA.  UPON  COMPLETION  0F',/5X, 

6 'CURVE  DIGITIZER,  THE  USER  MAY  UTILIZE  THE  PROGRAM  ',/5X, 

7 ' EMCADIN  TO  TRANSFORM  THE  DATA  INTO  A  FORM  WHICH  CAN',/5X, 

8 ' BE  USED  BY  THE  PROGRAM  EMCAD  WRITTEN  BY  DR.  M.  A.  MORGAN.', 
9////5X, 'PLEASE  PRESS  ENTER  TO  CONTINUE. ' ,//5X) 

120  F0RMATC////////////////////5X, 'INPUT  FILENAME' ,//////5X, 

1 ' INPUTFN  IS  THE  INPUT  FILENAME  OF  THE  FILE  CONTAINING' ,/5X, 

2 ' THE  DATA  OF  INTEREST  WHICH  WILL  BE  INTERPOLATED  ON  AND',/5X, 

3 'CONVERTED  TO  A  FORM  WHICH  CAN  SUBSEQUENTLY  BE  USED  BY',/5X, 

4 ' THE  CURVE  DIGITIZER  PROGRAM .', ///5X , 

5 ' PL  EASE  ENTER  THE  NAME  QF  THE  INPUT  FILE.  THE  EXTENSION ',/ 5X , 

6 ' 0 F  THE  FILENAME  MUST  BE  INCLUDED,  I.E.  FI L ENAME . DAT ’ , //5X ) 

130  FORMAT (////////////////////5X,  'OUTPUT  FILFNAME' ,//////5X, 

1 ' OUTPUTFN  IS  THE  OUTPUT  FILENAME  OF  THE  FILE  CONTAINING ',/ 5X , 

2 ' THE  DATA  WHICH  HAS  BEEN  INTERPOLATED  AND  CONVERTED  TO  A',/5X. 

3 ' FORM  WHICH  WILL  BE  USED  BY  THE  CURVE  DIGITIZER  PROGRAM .',/// 5X , 

4 ' PL  EASE  ENTER  THE  NAME  OF  THE  OUTPUT  FILE.  THE  EXTENSION' ,/5X, 

6 ' 0 F  THE  FILENAME  MUST  BE  .DAT,  I.E.  FI L ENAME . DAT ’,// 5X ) 

140  FORMAT (//////////////////// 5X, 'RESOLUTION' ,//////5X, 

l'DELTHE  IS  THE  USER  INPUT  VALUE  OF  THE  DESIRED  THETA', /5X, 

2 ' RESOLUTION  IN  DEGREES .',/// 5X , 

3 ' PL  EASE  ENTER  THE  DESIRED  DELTA  THETA  VALUE  IN  DEGREES .' ,//5X) 
CL0SEC7) 

CLOSEC8) 

END 
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APPENDIX  H 


MCADGRAF  USER'S  MANUAL 

Welcome  to  MCADGRAF.  This  program  was  written  to  be 
used  in  conjunction  with  the  FORTRAN  program  entitled  EMCAD. 
MCADGRAF  adds  graphics  features  to  support  the  analysis  of 
data  generated  and  formatted  by  EMCAD. 

A.  BACKGROUND 

MCADGRAF  was  developed  and  written  by  LT.  E.  M. 
Connolly  using  the  programming  language  of  Ryan-McFarland 
FORTRAN  (Appendix  M) .  Microcompatibles  GRAFMATIC  (Appendix 
M) ,  and  later  Jewell  Technologies  GRAFPLUS  (Appendix  M)  were 
chosen  as  the  commercially  available  screen  graphics  and 
hard  copy  support  for  EMCAD  output  graphics.  Any  questions 
concerning  colors,  formats,  changes  to  the  program,  etc. 
should  be  referred  to  the  respective  user's  manuals,  Dr.  M. 
A.  Morgan  or  LT.  E.  M.  Connolly  (Appendix  M) . 

B.  COMPUTER  SETUP 

MCADGRAF  was  written  for  the  following  computer 
configuration : 

*  IBM  AT  or  compatible,  including  80386  based 
compatibles 

*  hard  disk 

*  EGA  graphics  card 
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EGA  monitor 


* 

*  math  coprocessor  -  not  required 

*  memory  -  512K  RAM  minimum 

MCADGRAF  can  be  placed  in  any  directory,  at  the  user's 
discretion,  however,  it  is  highly  recommended  that  MCADGRAF 
and  EMCAD  be  placed  in  the  same  directory.  MCADGRAF  input 
is  taken  directly  from  EMCAD  output.  This  output  is 
required  to  be  in  the  same  directory/ sub-directory  as 
MCADGRAF,  or  on  the  same  floppy  disk  if  using  a  floppy 
system  instead  of  a  hard  disk.  The  only  other  requirement 
for  initial  setup  deals  with  printer  initialization.  In 
order  to  obtain  a  printer  graphical  output,  GRAFPLUS  must  be 
installed  on  the  computer,  configured  for  the  specific 
printer  of  use,  prior  to  running  MCADGRAF.  An  EPSON  LX- 
800  series  printer  was  used  when  writing  this  code.  This 
printer,  as  with  any  other  printer  to  be  used,  was 
configured  by  inserting  two  lines  in  the  AUTOEXEC.BAT  file. 
The  first  of  these  lines  was  a  path  command  to  tell  the 
computer  where  to  find  the  appropriate  files  for  GRAFPLUS. 
Assuming  all  the  files  for  GRAFPLUS  are  contained  in  a 
directory  on  the  C  disk  entitled  Grafplus,  an  appropriate 
path  entry  would  be 

path  c: \;c: \Grafplus 
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This  setup  will  also  work  with  other  IBM  graphics  compatible 
printers  such  as  the  Okidata  ML  -  192/193.  The  second  of 
these  lines  inserted  actually  established  the  initialization 
values  for  the  printer  to  be  used.  The  second  line  must  be 
inserted  anywhere  in  the  AUTOEXEC.BAT  after  the  path 
command.  The  printer  initialization  used  for  the  LX-800  was 

GRAFPLUS  =  #115 

The  equals  sign  (=)  directs  the  printer  to  print  the 
graphics  horizontally.  This  means  that  the  top  of  the  image 
is  at  the  top  left  corner  of  the  page.  For  vertical 
printing,  the  object  appears  to  be  rotated  with  the  top  of 
the  picture  at  the  bottom  left  corner  of  the  page.  The 
pound  sign  (#)  directs  the  printer  to  print  the  text  screens 
as  graphics  screens.  The  I  conveys  the  information  that  the 
screen  mode  is  EGA  -  o40x350  resolution,  with  16  colors 
available.  The  number  15  is  the  number  assigned  by  GRAFPLUS 
to  correspond  with  the  EPSON  LX-800.  Other  options 
available  are  the  capabilities  to  print  in  color  and  to 
reverse  printing.  All  of  these  options  are  discussed  in 
much  greater  detail  on  pages  3  through  8  of  the  GRAFPLUS 
user's  manual.  As  a  note  of  caution,  some  GRAFPLUS  printer 
definitions  are  defined  such  that  "graying"  is  used  on  a 
black  and  white  print-out  to  represent  different  colors.  On 
the  LX-800,  the  "graying"  effect  reduced  printer  clarity  and 
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resolution  dramatically.  This  problem  was  rectified  by 
running  the  EDITGRAF.COM  program  contained  in  GRAFPLUS , 
answering  all  the  questions  without  changes  until  the 
question  "Do  you  wish  to  print  in  Black  &  White  only? 
(disable  gray  scale)  "  appeared  on  the  screen.  Changing  the 
answer  from  no  to  yes  removed  the  "graying"  and  improved  the 
resolution.  A  much  more  detailed  explanation  of  EDITGRAF 
can  be  found  on  pages  9  through  21  of  the  GRAFPLUS  user's 
manual.  Once  these  changes  to  the  printer  definitions  and 
AUTOEXEC.BAT  have  been  made,  it  is  necessary  to  re-boot  the 
computer  to  initialize  these  values.  This  re-boot  can  be 
accomplished  by  either  turning  the  machine  off,  waiting  15 
to  20  seconds,  and  turning  it  back  on,  or  by  hitting  the 
Ctrl,  Alt  and  Del  keys  all  at  the  same  time. 

C.  PROGRAM  INPUT 

As  previously  mentioned,  MCADGRAF  obtains  its  input  from 
the  output  of  EMCAD.  This  EMCAD  output  is  four  files  in  the 
form  of  " FILENAME ".TET,  "FILENAME" . TEP,  "FILENAME" . TMT  and 
"FILENAME" .TMP.  The  quotation  marks  indicate  that  the 
FILENAME  is  a  value  that  was  previously  chosen  by  the  user 
while  running  EMCAD.  The  four  extensions  of  FILENAME  have 
significant  meaning.  The  TET  refers  to  TE  incidence, 
F-Theta;  the  TEP  refers  to  TE  incidence,  F-Phi;  the  TMT 
refers  to  TM  incidence,  F-Theta;  the  TMP  refers  to  TM 
incidence,  F-Phi.  The  values  for  these  four  files  were 
taken  directly  from  a  much  larger,  more  comprehensive  EMCAD 
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output  data  file.  The  scattering  calculation  output  and 
graphics  relationship  is  depicted  in  Figure  H-l. 

D.  PROGRAM  EXECUTION 

MCADGRAF  is  a  menu  driven  output  graphics  program.  Once 
the  correct  disk  or  directory  has  been  selected,  MCADGRAF 
can  be  invoked  by  typing 

MCADGRAF 

at  the  DOS  prompt.  The  program  will  load  itself,  and  the 
user  will  see  the  following  screens: 

*  Welcome  To  MCADGRAF  (Figure  H-2) 

This  screen  presents  some  brief  background  information. 

*  EMCAD  Output  Graphics  Presentation,  Part  1 

(Figure  H-3) 

This  menu  affords  the  user  the  opportunity  to  exit  the 
program  in  the  event  that  the  printer  initialization  has  not 
been  accomplished. 

*  EMCAD  Output  Graphics  Presentation,  Part  2 

(Figure  H-4) 

This  screen  explains  the  means  of  obtaining  an  output  graph 
on  the  printer.  The  user  will  be  presented  with  four  graphs 
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SCATTERING  CALCULATION 
OUTPUT  AND  GRAPHICS 


ttWHHHHHHHt  WELCOflE  10  NCADGRAF  WHHHHHHHW 

THIS  PROGRAM  WAS  WRITTEN  TO  BE  USED  IN  CONJUNCTION  WITH 
THE  FORTRAN  PROGRAM  ENTITLED  ENCAD.  NCADGRAF  ADDS 
GRAPHING  FEATURES  TO  SUPPORT  THE  ANALYSIS  OF  DATA 
GENERATED  AND  FORMATTED  BY  ENCAD.  ENCAD  WAS  WRITTEN  BY 
DR.  N.  A.  MORGAN,  WITH  MODIFICATIONS  BY  HIMSELF  AND 
LT  E.  N.  CONNOLLY.  NCADGRAF  WAS  DEUELOPED  AND  WRITTEN  BY 
LT  E.  N.  CONNOLLY  TO  HELP  FULFILL  THESIS  REQUIREMENTS  WHILE 
UNDER  THE  GUIDANCE  AND  THE  SUPPORT  OF  DR.  NORGAN.  NCADGRAF 
WAS  WRITTEN  AND  COMPILED  USING  RYAN  NCFARLAND  FORTRAN, 
GRAFNATIC,  AND  GRAFPLUS.  ANY  QUESTIONS  CONCERNING  COLORS, 
FORMAT,  CHANGES  10  THE  PROGRAM,  ETC...  SHOULD  BE  REFERRED 
TO  THE  RESPECTIVE  USERS  NANUALS,  DR.  N.  A.  NORGAN,  OR 
LT  E.  H.  CONNOLLY. 

PLEASE  PRESS  ENTER  TO  CONTINUE. 


FIG. 


H-2. 


WELCOME  TO  MCADGRAF 
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EHCAD  OUTPUT  GRAPHICS  PRESENTATION 


NOTES  CONCERNING  GRAPHICAL  HARDCOPY  OUTPUT,  PART  1 

This  progran  allows  the  user  to  obtain  printer  hardcopies 
of  the  graphs  being  viewed.  Hardcopy  graphic  outputs  are 
obtained  through  the  use  of  the  progran  entitled  GRAFPLUS. 
Proper  utilization  of  HCADGRAF  requires  prior  setup  of 
GRAFPLUS.  Directions  for  proper  setup  are  contained 
in  the  GRAFPLUS  users  nanual. 

1.  CONTINUE  WITH  THE  HEX!  SCREEN 

2.  EXIT  THE  PROGRAM 

Indicate  your  selection  by  the  choice  of  1  or  2. 


FIG.  H-3.  EMCAD  OUTPUT 
GRAPHICS  PRESENTATION  PART  1 
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ENCAD  OUTPUT  GRAPHICS  PRESENTATION 


NOTES  CONCERNING  GRAPHICAL  HARDCOPY  OUTPUT,  PART  2 

The  user  nakes  the  decision  as  to  whether  a  printer 
graphical  output  is  desired  after  viewing  the  graph.  Once 
a  graph  is  drawn  on  the  screen,  the  user  nay  request  a 
printer  hardcopy  by  typing  either  a  "P"  or  "p"  and  pressing 
the  "ENTER'’  key.  Selection  of  any  other  key  will  simply 
advance  to  the  next  graph.  Spacing  on  8  1/2  by  11  paper 
allows  up  to  2  graphs  per  page. 


PLEASE  PRESS  ENTER  TO  CONTINUE. 


FIG.  H-4.  EMCAD  OUTPUT 
GRAPHICS  PRESENTATION  PART  2 
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for  each  data  run.  As  each  graph  is  presented  on  the 
screen,  the  user  is  given  the  opportunity  to  obtain  a 
printer  copy  of  the  screen.  The  user  decides  to  obtain  a 
printer  copy  by  depressing  a  P  or  a  p,  followed  by  ENTER. 
The  reminder  line  telling  the  user  what  to  do  to  produce  a 
hardcopy  is  erased,  and  the  graph  is  sent  to  the  printer. 
Any  other  character  input,  followed  by  ENTER  will  simply 
advance  to  the  next  screen.  Setting  up  the  paper  in  the 
printer  is  the  sole  responsibility  of  the  user  for  each  and 
every  graph.  If  one  graph  is  to  placed  on  every  page,  it  is 
recommended  that  the  paper  be  advanced  a  few  inches  so  that 
the  graph  is  printed  closer  to  the  center  of  the  page.  For 
the  LX— 800 ,  it  is  possible  to  place  2  graphs  on  every  page, 
if  the  paper  is  positioned  properly  in  the  printer.  The 
best  advice  is  to  experiment. 

*  Data  File  Name  Input  Menu  (Figure  H-5) 

If  the  data  file  to  be  used  was  generated  by  EMCAD,  the  user 

need  only  choose  option  1.  If  the  file  to  be  used  is  a 

self-generated  file,  it  must  be  in  the  exact  format  of  an 

EMCAD  generated  file,  and  it  must  have  an  extension  of  .TET, 

.TEP,  . TMT  or  . TMP  .  If  the  file  is  not  in  the  correct 

> 

format,  or  if  the  file  does  not  contain  the  proper 
extension,  the  user  must  exit  the  program  and  correct  the 
problem  prior  to  continuing. 
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ENCAD  OUTPUT  GRAPHICS  PRESENTATION 
DATA  FILE  NAME  INPUT  MENU 

If  utilizing  a  file  op  series  of  files  created  by  ENCAD, 
please  select  option  1.  If  utilizing  any  file  other  than 
one  created  by  ENCAD,  please  ensure  that  the  file  is  in 
proper  format,  and  that  the  filename  extension  is  one  of 
the  four  extensions  utilized  by  this  program.  These 
extensions  are:  .tet  .tep  .tnt  and  .tup 
If  you  are  unsure  as  to  whether  your  file  is  in  the  proper 
foroat  with  the  proper  extension,  please  exit  the  progran 
and  verify  before  continuing. 

1.  INPUT  THE  FILE  OF  USER  CHOICE 

2.  EXIT  THE  PROGRAM 

Indicate  your  selection  by  the  choice  of  1  or  2. 


FIG.  H-5.  DATA  FILE  NAME 
INPUT  MENU 
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*  Data  File  Name  Input  (Figure  H-6) 

This  menu  prompts  the  user  for  the  filename  of  the  input 
data  file.  Since  the  extension  is  limited  to  four  cases, 
and  has  previously  been  selected  in  the  Data  File  Type 
Selection  Menu,  no  extension  should  be  added  to  the  filename 
to  be  typed  in. 

*  Data  File  Type  Selection  Menu  (Figure  H-7) 

As  previously  discussed,  MCADGRAF  requires  that  the  four 
EMCAD  output  data  files  be  located  in  the  same  directory. 
This  menu  gives  the  user  the  chance  to  choose  which  of  the 
four  files  the  data  will  come  from. 

This  menu  also  affords  the  user  the  opportunity  to  change 
the  default  setting  of  the  aspect  ratio.  The  menu  for  this 
change  is  shown  in  Figure  H-8 .  Aspect  ratio  is  a  term  that 
relates  the  number  of  available  pixel  columns  on  the  screen 
to  the  number  of  available  pixel  rows  on  a  screen.  Since 
most  computer  screens  are  not  square,  this  representation 
becomes  very  important.  In  simple  terms,  it  deals  with 
making  a  square  look  like  a  square,  and  not  like  a 
rectangle,  or  making  a  circle  look  like  a  circle,  and  not  an 
oval . 

GRAFPLUS  comes  with  a  program  to  help  determine  the 
aspect  ratio  for  any  monitor.  In  the  case  of  the 
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ENCAD  OUTPUT  GRAPHICS  PRESENTATION 
DATA  FILE  NANE  INPUT 

Please  input  the  nane  of  the  file  containing  data  that 
you  wish  to  work  with .  Please  input  only  the  file  name 
i.e.  do  not  include  the  file  extension.  The  file 
extension  will  be  selected  in  the  DAIA  FILE  IYPE 
SELECTION  MENU. 


FIG.  H-6.  DATA  FILE  NAME  INPUT 
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ENCAD  OUTPUT  GRAPHICS  PRESENTATION 
DATA  FILE  TVPE  SELECTION  HENU 


Please  select  the  type  of  data  file  containing  the  information  of 
interest,  ENCAD  automatically  created  these  four  data  file  type  s. 


1. 

TE  INCIDENCE, 

F-IHETA 

2. 

TE  INCIDENCE, 

F-PHI 

3. 

IN  INCIDENCE, 

F-THETA 

4, 

TH  INCIDENCE, 

F-PHI 

5. 

CHANGE  DATA  FILE 

6. 

CHANGE  ASPECT  RATIO 

7. 

EXIT  THE  PROGRAM 

Indicate  your  selection  by  the  choice  of  1,  2,  3,  4,  5,  b  or  7 


FIG.  H-7.  DATA  FILE  TYPE 
SELECTION  MENU 
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**********  aspect  ratio  ********** 

IN  ORDER  TO  AUOID  DISTORTION  ON  THE  SCREEN,  IT  IS  NECESSARY 
TO  DEFINE  THE  NUMBER  OF  ROWS  PER  INCH  AND  THE  NUMBER  OF 
COLUMNS  PER  INCH.  THESE  DEFINITIONS  CHANGE  AS  THE 
RESOLUTION  OF  THE  SCREEN,  AS  WELL  AS  THE  TYPE  OF  SCREEN 
CHANGE.  THIS  PROGRAM  WAS  WRITTEN  IN  NODE  16  (EGA)  USING 
A  GB-1  UIDEO  BOARD  AND  A  NEC  MULTISYNC  MONITOR.  THE 
OPTIMUM  ASPECT  RATIO  FOR  THIS  CONFIGURATION  WAS  .65 


PLEASE  INPUT  ASPECT  RATIO. 


FIG.  H-8.  ASPECT  RATIO 
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development  of  this  program,  the  aspect  ratio  given  by 
GRAFPLUS  delivered  erroneous  results.  As  such,  this  screen 
of  MCADGRAF  allows  the  user  to  input  and  change  the  aspect 
ratio  of  choice.  For  users  unsure  of  their  aspect  ratio, 
start  with  an  initial  value  of  1,  and  work  up  or  down  with 
increments  of  .1  or  .2  until  close  to  the  desired  results, 
then  fine  tune  as  much  as  needed. 

*  Selection  of  Incidence  Angle  (Figure  H-9) 

One  of  the  input  parameters  for  EMCAD  was  the  number  of 
incident  fields,  and  their  respective  angles.  MCADGRAF 
reads  the  input  data  file,  and  presents  in  this  menu  the 
number  of  incidence  angles,  and  their  respective  angles  in 
degrees.  It  is  up  to  the  user  to  select  the  incidence  angle 
of  interest. 

*  Selection  of  Phi  Angle  (Figure  H-10) 

As  with  incidence  angles,  phi  angles  were  input  parameters 
for  EMCAD.  MCADGRAF  presents  in  this  menu  the  number  of  phi 
angles,  and  their  respective  angles  in  degrees.  It  is  up  to 
the  user  to  select  the  phi  angle  of  interest. 
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ENCAD  OUTPUT  GRAPHICS  PRESENTATION 

SELECTION  OF  INCIDENCE  ANGLE 
The  data  file  you  have  selected  contains  data  for  5 
incidence  angles.  Please  select  your  choice  from  the 
following  nenu.  Note,  you  nay  only  work  uith  one 
incidence  angle  at  a  tine. 

1  INCIDENCE  ANGLE  =  8.00 

2  INCIDENCE  ANGLE  =  45,00 

3  INCIDENCE  ANGLE  =  90.00 

4  INCIDENCE  ANGLE  =  135.00 

5  INCIDENCE  ANGLE  =  188.00 

6  RETURN  TO  THE  DATA  FILE  TVPE  SELECTION  HENU 

7  EXIT  THE  PROGRAN 

INDICATE  YOUR  SELECTION  BY  THE  CHOICE  OF  A  NUHBER  I.E.  1  2  3  ...) 


FIG.  H-9.  SELECTION  OF 
INCIDENCE  ANGLE 
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ENCAD  OUIPUI  GRAPHICS  PRESENTATION 
SELECTION  OF  PHI  ANGLE 

The  data  file  you  have  selected  contains  data  for  3 
phi  angles.  Please  select  your  choice  frow  the  following 
nenu.  Mote,  you  nay  only  work  with  one  phi  angle  at  a  tine. 

1  PHI  ANGLE  =  0.00 

2  PHI  ANGLE  =  90.08 

3  PHI  ANGLE  =  180.00 

4  RETURN  TO  THE  DATA  FILE  IVPE  SELECTION  NENU 

5  EXIT  THE  PROGRAM 

INDICATE  YOUR  SELECTION  BY  THE  CHOICE  OF  A  HUMBER  I.E.  1  2  3 


FIG.  H-1D.  SELECTION  OF 


PHI  ANGLE 


E .  PROGRAM  OUTPUT 


For  each  separate  data  set,  MCADGRAF  produces  four 
different  graphs. 

*  Magnitude  Versus  Theta  (Figure  H-ll) 

As  can  be  seen  in  the  first  line  of  the  title  for  this 
graph,  as  well  as  the  other  three  graphs,  a  label  specific 
to  this  data  set  appears.  This  label  was  used  as  an  input 
parameter  to  EMCAD  to  assist  the  user  in  differentiating 
this  set  of  data  from  every  other  set.  In  the  third  line, 
the  values  of  the  incidence  angle  and  phi  angle,  as  well  as 
the  type  of  incidence  are  listed.  Either  F-Phi  or  F-Theta 
appears  in  the  fourth  line  followed  by  a  short  description 
of  the  type  of  graph.  The  line  at  the  bottom  referring  to 
hardcopy  is  simply  a  reminder  to  the  user  of  what  is  needed 
to  request  a  hardcopy  output.  This  line  is  erased  before 
the  printer  graph  is  produced.  This  first  graph  depicts  the 
scattering  magnitude  plotted  versus  the  180  degrees  of 
Theta . 

*  Phase  Versus  Theta  (Figure  H-12) 

The  only  change  in  this  graph  from  the  first  graph  is  that 
phase  instead  of  magnitude  is  being  graphed  versus  Theta. 
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Composite  Body 

Scattered  Field  Complex  Amplitudes 
111  Incidence  =  135.  Degrees  Phi  =  180.  Degrees 
Phase  F-Theta  Phase  Versus  Iheta 


PIC.  H-12.  PHASE  VERSUS 


THETA 


*  Magnitude  Versus  Theta,  Polar  Plot 

(Figure  H-13) 

This  graph  represents  the  same  data  as  the  first  graph,  but 
from  a  different  viewpoint.  The  first  graph  was  drawn  using 
rectangular  coordinates.  This  graph  utilizes  polar 
coordinates.  This  graph  type  is  the  one  predominately  used 
in  the  industry  to  illustrate  scattering. 

*  Magnitude  Versus  Theta,  Polar  Plot 
Axis  Rotated  Minus  Ninety  Degrees 

(Figure  H-14) 

The  only  change  in  this  graph  from  the  third  graph  is  the 
rotation  of  the  axis.  The  limitations  of  GRAFMATIC  made  it 
difficult  to  generate  the  data,  maintain  spacing 
restrictions  and  still  display  all  the  data  necessary  to 
adequately  describe  the  situation.  Rotating  the  axis  ninety 
degrees  displays  the  data  slightly  differently  from  the 
third  graph,  giving  the  user  a  bit  more  useful  information. 


Offset  Sphere  Validation  Ran 
Scattered  Field  Complex  Amplitudes 
IE  Incidence  =  90.  Degrees  Phi  =  0.  Degrees 

Polar  Plot,  F-Phi  Magnitude  Versus  Theta 


FIG.  H-13.  MAGNITUDE  VERSUS 
THETA.  POLAR  PLOT 


168 


Composite  Body 

Scattered  Field  Coup  lex  Amplitudes 
If!  Incidence  =  135.  Degrees  Phi  =  1B8.  Degrees 
Polar  Plot,  F- Theta  Magnitude  Versus  Theta 
Axis  Rotated  Minus  Ninety  Degrees 


S  Theta=90 


Theta  =  0  Theta=180 

Hardcopy  — >  P  or  p 


FIG.  H-14.  MAGNITUDE  VERSUS 
THETA,  POLAR  PLOT,  AXIS 
ROTATED  MINUS  NINETY  DEGREES 
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PPENDIX  I 


MCaDGRAF  SOURCE  CODE 


PROGRAM  MCADGRAF. FOR 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX** 

THIS  PROGRAM  WAS  WRITTEN  TO  BE  USED  IN  CONJUNCTION  WITH  THE 
FORTRAN  PROGRAM  ENTITLED  EMCAD.  MCADGRAF  ADDS  GRAPHING  FEATURES 
TO  SUPPORT  THE  ANALYSIS  OF  DATA  GENERATED  AND  FORMATTED  BY  EMCAD 
EMCAD  WAS  WRITTEN  BY  DR.  M.  A.  MORGAN,  WITH  SUBSEQUENT  UPDATE  BY 
HIMSELF  AND  LT  E.  M.  CONNOLLY  WHILE  AT  THE  NAVAL  POSTGRADUATE 
SCHOOL  (MODS  MAR  87  -  MAR  88).  MCADGRAF  WAS  DEVELOPED  AND  WRITTEN 
BY  LT  E.  M.  CONNOLLY  TO  HELP  FULFILL  THESIS  REQUIREMENTS  WHILE 
UNDER  THE  GUIDANCE  AND  THE  SUPPORT  OF  DR.  MORGAN.  MCADGRAF  WAS 
WRITTEN  AND  COMPILED  USING  RYAN  MCFARLAND  FORTRAN,  GRAFMATIC,  AND 
GRAFPLUS.  ANY  QUESTIONS  CONCERNING  COLORS,  FORMAT,  CHANGES  TO  THE 
PROGRAM,  ETC...  SHOULD  BE  REFERRED  TO  THE  RESPECTIVE  USERS 
MANUALS,  DR.  M.  A.  MORGAN,  OR  LT  E.  M.  CONNOLLY. 

A  NOTE  OF  CAUTION,  TO  OBTAIN  CORRECT  PRINTOUTS  FROM  MCADGRAF, 
GRAFPLUS  MUST  3E  PROPERLY  CONFIGURED  TO  SUPPORT  TOUR  HARDWARE. 
QUESTIONS  CONCERNING  THE  PROPER  HARDWARE  SETUP  SHOULD  BE  REFERRED 
TO  THE  GRAFPLUS  USER'S  MANUAL. 


xxx**xxx**xxxxxxxxxxxxx***xxx*xxxxx**xxx*xxxxxxxxx***x**xxx**x**x* 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
PROGRAM  STATEMENT  DEFINITIONS 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


QSMODE  (mode) 


Clears 

the 

screen  and 

sets 

the 

mode 

of  the  1 

display 

May  be  used 

in  either 

the 

text  or 

graphic 

s  mode.  The  desire 

d  val 

ue  of 

mode  i: 

chosen 

from 

the  following 

table 

MODE 

RESOLUTION  (CxR) 

COLOR  T/G 

0 

40 

X  25 

B/W 

T 

1 

40 

X  25 

C 

T 

2 

80 

X  25 

B/W 

T 

3 

80 

X  25 

C 

T 

4 

320 

X  200 

C 

G 

5 

320 

X  200 

B/W 

G 

6 

640 

X  200 

B/W 

G 

7 

80 

X  25 

3/W 

T  ( M 

13* 

320 

X  200 

C 

G 

14* 

640 

X  200 

C 

G 

15* 

640 

X  350 

50HZ 

MONO 

16* 

640 

X  350 

4/16 

COLORS 

X - > 

Requires  EGA 

carij 

QPLOT  (  jcoll,  Jcol2,  jrowl,  jrow2,  xmin,  xmax,  ymin,  ymax, 
xorg,  yorg,  iopt,  yoverx,  aspect) 
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j  col  1 
j  col2 
0  rowl 
j  r  ow2 
xmin 

xmax 

ymi  n 

ymax 

xorg 

yorg 

i  opt 
yoverx 

aspect 

QXAXIS 

xst 

xf  in 
xmaj or 

minor 


label 

ndec 

QYAXIS 


Plot  initialization  routine:  sets  plot  window, 
corresponding  x  and  y  ranges,  and  scaling 
factors . 

left  hand  Pixel  column  defining  plot  region 
right  hand  pixel  column  defining  plot  region 
lower  pixel  row  defining  plot  region 
upper  pixel  row  defining  plot  region 
value  of  x  variable  assigned  to  correspond  to 
column  jcoll 

value  of  x  variable  assigned  to  correspond  to 
column  jco!2 

value  of  y  variable  assigned  to  correspond  to 
row  jrowl 

If  iopt  =  0,  then  ymax  is  the  value  of  the  y 
variable  assigned  to  row  =  jrow2.  If  iopt  =  1, 
then  on  input  ymax  is  the  assumed  maximum  value 
of  the  y  variable  required  by  the  plot, 
vertical  (y)  axis  intersects  the  horizontal  (x) 
axis  at  x  =  xorg 

horizontal  (x)  axis  intersects  the  wertical  (y) 
axis  at  y  =  yorg 

sets  scaling  method  option,  see  ymax 
number  of  y  unit  changes  for  a  unit  change  in 
value  of  x,  ie. .  the  change  in  the  y  coordinate 
per  inch  divided  by  the  change  in  the  x 
coordinate  per  inch 

aspect  ratio  for  the  user's  CRT,  for  NEC 
Multisync,  aspect  =0.6  for  mode  =  16 


xst,  xfin,  xmajor,  minor,  label,  ndec) 

Draws  the  x  axis  along  the  horizontal  direction, 
including  tic  marks  and  tic  mark  labels 

value  of  x  at  the  start  of  the  x  axis; 

If  a  logarithmic  scale  is  being  used,  then 
xst  and  xfin  must  be  one  of  the  following 
with  xfin  >  xst  : 

0.001,  0.01,  1.0,  10.0,  100.0,  1000.0,  10000.0 

value  of  x  at  the  end  of  the  x  axis 

distance  between  the  major  (large)  tic  marks; 

If  xmajor  is  positive/negative,  then  the  tic 
marks  will  be  placed  above/below  the  x  axis 
respectively.  If  xmajor  is  zero,  then  no  tic 
marks  are  inserted 

integer  flag  for  use  of  minor  tic  marks  or 
logarithmic  scale;  If  minor  =  0,  then  no  tic 
marks  are  inserted.  If  minor  =  +1,  then  minor 
tic  marks  are  inserted  halfway  between  the 
major  tic  marks,  but  are  not  labeled.  If 
minor  =  -1,  then  a  logarithmic  set  of  tic  marks 
are  inserted. 

integer  flag  for  auto  labeling  of  tic  marks; 

If  label  =  0,  then  there  is  no  auto  labeling. 

If  label  =  1,  then  the  ma  or  tic  marks  are 
automatically  supplied  with  a  label.  If  label 
=  -1,  then  the  origin  is  labeled  as  well, 
number  of  decimal  places  for  auto  labeling; 

The  format  of  the  label  is  FlO.n  with  the 
decimal  point  lined  up  with  the  major  tic  marks. 


yst,  yfin,  ymajor,  minor,  label,  ndec  ) 


171 


nononnoononnoooooonnonoonnooooonnnnnoonoonononnoonooononnooDonnnoo 


yst 
yf  i  n 
yn a :  or 

minor 

label 

ndec 

QPREG  ( 

i  reg 
i  va  1 

QPTXT  ( 

ncha  r 
char 

kol  or 
i  col 
i  row 

QPTXTB 


QPTXTC 


QCMOV  ( 


QSETUP 


ndots 


lko 1  or 
isymbl 


Creates  the  y  axis  for  a  plot. 

value  of  y  at  the  start  of  the  y  axis 
value  of  y  at  the  end  of  the  y  axis 
distance  between  major  (large)  tic  marks  on 
the  y  axis 

integer  flag  for  use  of  minor  tic  marks  or 
logarithmic  scale,  see  QXAXIS 

integer  flag  for  auto  labeling  of  tic  marks, 
see  QXAXIS 

number  of  decimal  places  for  auto  labeling, 
see  QXAXIS,  0  <=  ndec  <=  3 

ireg,  ival) 

Changes  the  color  table  by  modifying  the 
Pallette  Registers,  EGA  card  only 

pallette  register  index  (0-15) 
value  to  be  placed  in  register  (0-63) 


nchar,  char,  kolor,  icol,  irow  ) 

Writes  plot  text.  The  character  variable  char 
of  length  nchar  is  written  at  (  column,  row  )  - 
(  icol,  irow  )  in  specified  color. 

length  of  character  variable 

character  variable  corresponding  to  desired 

text 

specifies  color  of  the  text 
starting  text  column  of  text 
text  row  containing  text 

(  nchar,  char,  kolor  ) 

Writes  plot  text  at  position  B.  A  character 
variable  cha-  of  length  nchar  is  written  at 
position  ("B")  just  beyond  the  positive  X  axis 

(  nchar,  char,  kolor  ) 

Writes  plot  text  at  position  C.  A  character 
variable  char  of  length  nchar  is  written  at 
position  ("C")  just  above  and  centered  on  the 
Y  axis 

icol,  irow  ) 

Moves  the  cursor  to  a  specified  column  and 
row 

(  ndots,  lkolor,  isymbl,  klrsym  ) 

Sets  up  line  type  (solid  or  dashed)  and  symbols 
to  be  used  by  QTABIE  when  plotting  data 

Number  of  pixels  (dots)  specifying  length  of  the 
dash.  If  ndots  -  Q,  then  a  solid  line  is  drawn. 
The  input  ndots  -  1  is  not  allowed.  By  varying 
ndots  >  -  2,  you  may  create  a  line  composed  of 
variable  length  dashes 
Specifies  color  of  the  line 

Specifies  symbol  to  be  inserted  at  data  points 
QTABL.  If  isymbl  =  -1,  then  a  cross  is  drawn  at 
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k 1 rsyn 

QTA5L  (  itype. 


1  type 


npt 

xarray 
yar ray 
selectl 


sea  1 exy 
xymax 

f namel 

fnane2 

xtensi on 


gra f pri nt 

dummy 
ti t 1  el 
ti t 1 e2 
t :.  1 1  e  3 
+l tl e4 
xarray 


'/array 


xamax 

yamax 

ganax 

Pi 

f  name 
numsets 


nunpts 
a  ,  b  ,  c 


l ncangi 


the  data  points.  If  isymbl  -  -2,  then  no  symbol 
is  drawn.  Otherwise,  the  character  with  the 
ascii  number  "l synb 1 "  is  drawn 
Specifies  color  of  the  symbol 

npt,  xarray,  '/array  ) 

Plots  a  tabulated  function 

Specifies  point  or  line  plot.  When  itype  =  0, 
then  the  plot  is  a  series  of  symbols  at  the 
tabulated  points.  When  itype  =  1,  tne  symbols 
are  joined  by  straight  lines 
The  number  of  points  in  the  function  table 
X  values  in  the  function  table 
Y  values  in  the  function  table 
A  character  string  used  to  hold  the  irput 
values  of  incidence  angle  and  phi  angle 
selection 

A  scale  factor  used  in  three  of  the  four  graphs 
The  maximum  scaled  value  out  of  ail  the  scaled 
X  values  and  all  of  the  scaled  Y  values 
A  character  string  used  to  hold  the  value  of 
the  input  data  file  name  without  extension 
A  character  string  used  to  hold  the  value  of 
the  input  data  file  name  with  extension 
A  character  string  used  to  hold  the  value  of 
the  input  data  filename  extension,  i.e. 

.TET,  . TEP .  . TMT,  or  . TMP 

A  character  string  used  to  allow  the  user  to 
print  out  a  graph  from  the  screen 
An  input  variable  used  to  advance  the  screen 
The  first  line  of  the  title  for  each  graph 
The  second  line  of  the  title  for  each  graph 
The  third  line  of  the  title  for  each  graph 
The  fourth  line  of  the  title  for  each  grapn 
Data  files  contain  numbers  in  THET,  MAG,  PHASE 
format.  This  polar  coordinate  format  can  be 
transformed  to  rectangular  coordinate  format. 
XARRAY  is  the  array  established  to  hold  all  of 
the  transformed  X  values. 

Data  files  contain  numbers  in  THET,  MAG,  PHASE 
format  This  polar  coordinate  format  can  be 
transformed  to  rectangular  coordinate  format. 
YARRAY  is  the  array  established  to  hold  all  of 
the  transformed  Y  values. 

The  maximum  value  of  the  X  array  (XARRAY) 

The  maximum  vulue  of  the  Y  array  (YARRAY) 

The  maximum  value  of  the  X  array  (XARRAY) 

and  the  Y  array  (YARRAY) 

A  constant  a pproxi ma tel y  equal  to  22/7 
A  character  string  which  equals  either  F-Phi  or 
F-Theta 

A  number  used  for  indexing  purposes,  eaual  to 
number  of  incidence  angles  times  the  number  of 
phi  angles 

The  number  of  theta  angles  per  given  incidence 
angle  and  phi  angle 

Dummy  variables  used  for  incrementation, 

a  Index  for  incidence  angles 

b  Index  for  phi  angles 

c  Index  for  theta  angles 

A  dummy  variable  used  to  read  and  hold  the 
value  of  the  incidence  angle  presently  being 
utilized 
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incangl ( a ) 

A  matrix  containing  the  values  of  all  the 
incidence  angles 

i 

A  dummy  integer  used  as  a  counter 

phi  i 

A  dummy  variable  used  to  read  and  hold  the 
value  of  the  phi  angle  presently  peing  ut  ■.  eed 

incangl d 

A  dummy  variable  containing  an  incidence  -  gle 
value  used  for  comparison  purposes 

thetai 

A  dummy  variable  used  to  read  and  hold  the 
value  of  the  theta  angle  presently  being 
utilized 

magi 

The  magnitude  of  the  output  for  the  existing 
incidence  angle,  phi  angle,  and  theta  angle 

phasei 

A  dummy  variable  used  to  read  and  hold  the 
value  of  the  phase  angle  presently  being 
uti li zed 

ti tl e6 

A  dummy  character  string  used  to  hold  an  input 
title  to  be  used  in  an  output  graph 

ti tlo7 

A  dummy  character  string  used  to  hold  an  input 
title  to  be  used  in  an  output  graph 

tl  1 1  e8 

A  dummy  character  string  used  to  hold  an  input 
title  to  be  used  in  an  output  graph 

tarrayC  c ) 

A  matrix  containing  the  values  of  all  the  theta 
angles  in  degrees  for  a  given  incidence  angle 
and  phi  angle 

rmarray( c ) 

A  matrix  containing  tne  values  of  all  the 
magnitudes  for  a  given  incidence  angle 
and  phi  angle 

parray( c ) 

A  matrix  containing  the  values  of  all  the 
phase  angles  for  a  given  incidence  angle 
and  phi  angle 

ca leva  r 

A  variable  used  in  intermediate  ca 1 cula ti oi ns 
to  determine  a  graphical  scale  factor  (scalefac) 

sea  1 efac 

A  graphical  scale  factor 

dummyva  r 

A  dummy  variable  used  as  a  limit  in  scale  factor 
calculations  for  visual  graphic  asthetics 

theta (a , b,  c ) 

A  matrix  containing  the  values  of  all  the  theta 
angles  for  all  the  incidence  angles  and  all  the 
phi  angles 

mag(a,b,c) 

A  matrix  containing  the  values  of  all  the 
magnitudes  for  all  the  incidence  angles  and 
all  the  phi  angles 

phase(a,b,c) 

A  matrix  containing  the  values  of  all  the  phase 
angles  for  all  the  incidence  angles  and  all  the 
phi  angles 

rad( c  ) 

A  matrix  containing  the  values  of  all  the  theta 
angles  in  radians  for  a  given  incidence  angle 
and  phi  angle 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
VARIABLE  DECLARATIONS 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX***** 

INTEGER  A,  B.  C,  I,  IREG,  IVAL,  SCALEFAC 

INTEGER*2  ICOL,  IOPT,  IRON,  ISYMBl,  JC0L1,  JC0L2,  JR0W1 

INTEGERX2  JR0W2,  KI.RSYM,  K.OLOR,  LABEL,  LK.OLOR,  MINOR 

I  NT EGER *2  MODE,  NCHAR,  NDEC ,  SELECT1 

CHARACTERX2  SELECT,  XTENSION 

CHARACTERX3  GRAFPRNT 

CHARACTERX3  DUMMY,  FNAME1 ,  TITLE7 

CHARACTERx 1 2  FNAME2 

CHARACTERX20  CHAR 

CHARACTER#64  TITLE1,  T I TL  E2 ,  TITLE3,  TITLE4,  TITLE5 
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CHARACTER*64  TITLES .  TITL E6 ( 5 , 1 9 , 1 ) 

REAL  ASPECT ,  CALCVAR,  DUMMYVAR,  GAMAX,  GMAX,  INCANGI 
REAL  INCANGLD,  PHASEI,  PHII,  PI,  MAGI,  MOOTS,  MUMPTS 
REAL  NUM3ETS ■  RAD,  THFTAI 

REAl  XAMAX,  X  FI M ,  XMAJOR,  XMAX,  XMIN,  XORG,  XST 

REAL  YAMAX,  YFIM,  YMAJOR,  YMAX ,  YMIM,  YORG,  YOVERX,  YST 

REAL  AARRAY C  37  )  ,  BARRAY C  57  )  ,  INCANGL ( 5 ) ,  PHI ( 19  ) 

REAL  MAGC  5,19,37),  PHASEC 5 , 19 , 37  )  ,  THETA ( 5 , 1 9 , 37 ) 

COMMON/ONE/  RMARRAY (371,  SCAIEFAC,  TARRAY C  37 ) 

COMMON/TWO/  PARRAY ( 37  ) 

COMMON/THREE/  XARRAY ( 37 ) ,  YARRAY C  37  ) ,  SCALEXY,  XYMAX 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXS (XXXXXXXXXXXXXXXXXXXXXXXXXXX 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX* 

READ  DATA  FILES,  LOAD  ARRAYS,  CALCULATE  MAXIMUM  VALUES 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
ESTABLISH  INITIAL  VALUES 
pi  -  3.1416 
aspect  =  .65 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
ESTABLISH  INITIAL  SCREEN  MODE  AND  COLORS 
mode  -  3 

CALL  qsmodeCmode) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WELCOME  TO  MCADGRAF 

writeC*,103) 

read(x,100)  dummy 

CALL  qsmodeCmode) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
NOTES  CONCERNING  GRAPHICAL  OUTPUT  PART  1 
60  writeCX, 101  ) 

readC*,100)  select 
CALL  qsmodeCmode) 

ifCselect  ,ne.  '1'  .and.  select  .ne.  '2')  go  to  60 
ifCselect  .eq.  *2')  so  to  1 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

NOTES  CONCERNING  GRAPHICAL  OUTPUT  PART  2 

wri teC  *,102) 

read! *,100)  dummy 

CALL  qsmodeCmode) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
DATA  FILE  NAME  INPUT  MENU 
4  wri te ( * , 9 1  ) 

read(x,100)  select 

ifCselect  .ne.  ’I’. and.  select  .ne.  ’2')  go  to  4 
ifCselect  .eq.  '2')  go  to  1 
CAl L  asnodeCmode) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
INPUT  DATA  FILE  NAME 
45  wri teC  X , 92 ) 

readCx,100)  fnamel 
CALL  qsmodeCmode) 

xx*xxxxxx**xx**xxxxxxxxx*xx*xxxxxxx*xxx*xxxx*xxx***xx*xxxx**xxxx*x 
DATA  FILE  TYPE  SELECTION  MENU 
44  wri teC  * , 90  ) 

readC*,10Q)  select 
xtension  -  select 
CALL  qsmodeCmode) 

ifCselect  .ne.  '1'  .and.  select  .  ne .  '2*  .and.  select  .ne.  '3' 

1  and  select  .ne.  '4'  .and.  select  . ne .  '5'  .and  select  .  ne .  '6' 

2. and.  select  .ne.  '7')  go  to  44 
ifCselect  . eq .  '7')  go  to  1 

ifCselect  .eq.  '6')  go  to  52 
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iftselect  .eq.  ’5’)  go  to  45 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx** 
OPEN  DATA  FILE 
if (xtension  . eq .  *1')  then 

fname2  =  fnanel/'''  '  . tet ' 
else  if(xtension  eq.  *2')  then 
fnane2  =  fnamel// '  . tep  1 
else  lfCxtension  . eq .  *3')  then 

frame2  =  f nanel// 1 . tmt ' 
else 

fname2  =  f namel// ' . tmp 1 
endi  f 

CALL  qsmode(mode) 
open(4,  file=fname2) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

READ  DATA  FILE 

read! 4,100)  titleS 

r ead( 4 , 1 00 )  title7 

read(4,86)  numsets 

read(4,36)  numpts 

a  =  1 

incangld  =  -1000  . 
b  =  1 

do  16  i  =  1 ,  int( numsets ) 
c-  1 

read(4,86)  incangi 
lfCincangi  .ne.  incangld)  b  =  1 
iflincangi  .ne.  incangld)  a  =  a  +  1 
if(i  .eq.  1)  a  =  1 
lncangl(a)  =  incangi 
incangld  =  incangi 
read(4,86)  phii 
phi(b)  =  phii 
read(4,100)  ti 1 1 e6 ( a , b , c  ) 
do  17  c  -  1 , int ( numpts ) 
read(4,88)  thetai,  magi,  phasei 
theta(a,b,c)  -  thetai 
mag(a,b,c)  -  magi 
phase(a,b,c)  3  phasei 
17  continue 
b  =  b  +  1 
16  continue 
closet  4  ) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SELECT  DATA  FILE  SUBSET 

CI.E.  SELECT  INCIDENCE  ANGLE  AND  PHI  ANGLE) 

5  wri te( X , 93 )  a 
do  25  i =1 ,  a 

write(*,94)  i,incangl(i) 

25  continue 
write(x,95)  a  +  l>  a  +  2 
wntet  X ,  96  ) 
read'X,X)  selectl 
CALL  qsmode(mode) 

iftselectl  .ne.  1  .and.  selectl  .ne.  2  .and.  selectl  .ne.  3 
1 . and .  selectl  .ne.  4  .and.  selectl  .ne.  5  .and.  selectl  .ne.  6 
2  .and.  selectl  .ne.  7)  go  to  5 
iftselectl  .eq.  a+1)  go  to  44 
ifCselectl  .eq.  a+2)  go  to  1 
a  z  selectl 

6  wri te( * , 97 )  b-l 
do  26  i  =  l , b-l 

write(X,98)  i,phi(i) 

26  continue 

wri te( x , 99 )  b,  b+1 
wri te( x , 96  ) 
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rea 

d ( x , x )  selectl 

CAL 

L  qsmode 

( mode ) 

CAL 

L  qpregt 

0,32; 

) 

if( 

selectl 

.  eq  . 

b)  go  to 

44 

if( 

selectl 

.  eq  . 

b+1)  go 

to  1 

b  -  selectl 

C  *XXXXXX*XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX*XXXX*XXXXXX#X 

C  LOAD  ARRAYS 

c  =  1 

xamax  -  0.0 
yamax  =  0.0 
gamax  =  0.0 

title2  -  ti tl e6 ( a , b, c ) 
do  18  c  =  1 , int( numpts ) 
tarray(c)  =  theta < a , b , c  ) 
rmarr ay(c)  =  mag(a,b,c) 
parrayCc)  =  phase(a,b,c) 
rad  =  tarray(c)  *  pi  /  180 
xarray(c)  =  rmarray(c)  *  sin(rad) 
yarray(c)  =  rmarray(c)  x  cos(rad) 

C  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C  CALCULATE  MAXIMUM  VALUES 

if  (abs( xarray( c ) )  . gt .  xamax)  xamax  =  absCxarrayC c ) ) 
if  (abstyarray(c) )  .gt.  yamax)  yamax  -  abs( yarrayC c ) ) 

18  continue 
75  continue 

C  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
C  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
C  CALCULATE  SCALE  FACTORS 
gamax  =  abs(xamax) 

if  (yamax  . gt .  xamax)  gamax  =  abs(yamax) 
if  (abs(gamax  -  0)  .It.  .000001)  gamax  =  .000001 
calcvar  =  a  1 ogl 0 ( gamax ) 
scalexy  =  int(calcvar) 

if  (gamax  .It.  1.0)  scalexy  =  int(calcvar)  -  1 
scalefac  -  int(scalexy) 
dummyvar  =  gamax/ ( 1 Qxxscal exy ) 
xymax  =  int( dummyvar  +  1) 

C  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C  CLEAR  SCREEN 

mode  =  16 
CALL  gsmo de(mode) 

C  XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXjCXXXXXXX 

CALL  rectmvt(title2,title7,title8. numpts, aspect ) 

C  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C  CLEAR  SCREEN 

node  =  16 
CALL  qsmode(mode) 

C  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

CALL  rectpvt(title2,title7,title8, numpts, aspect) 

C  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C  CLEAR  SCREEN 

mode  -  16 
CALL  qsmode(mode) 

C  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

CALL  pola rmvt ( title2,title7,title8, numpts , aspect ) 

C  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C  CLEAR  SCREEN 

mode  =  16 
CALL  qsmode(mode) 

C  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

CALL  polarrot(title2,title7,title8, numpts, aspect ) 

C  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C  CLEAR  SCREEN 

mode  =  16 
CALL  qsmode(mode) 
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C  XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX#X 

goto  44 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
ASPECT  RATIO 
52  wn  t  e  (  *  ,  1  0  4  ) 

read( x  ,  x  )  aspect 
mode  =  3 

CALL  qsmode(mode) 
goto  44 

C  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

1  continue 
C  CLEAR  SCREEN 
mode  -  3 

CALL  qsmode(mode) 

C  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
80  format(6X,F5 . 0, 3X, Ell . 3,6X,F7 . 2) 

86  forma t ( F5 . 0  ) 

88  formate  F5.0,E11.3,F7.2) 

90  format(/5X, ' EMCAD  OUTPUT  GRAPHICS  PRESENTATION' ,/5X, 

1 ’ DATA  FILE  TYPE  SELECTION  MENU',//5X, 

2'Please  select  the  type  of  data  file  containing  the  information  of 

3' ,/5X, 

4'interest.  EMCAD  automatically  created  these  four  data  file  type 
5s. ' ,//8X, 

6 ' 1 .  TE  INCIDENCE,  F-THETA ' , //8X, 

7 ' 2 .  TE  INCIDENCE,  F-PHI  ',//8X, 

8 ' 3 .  TM  INCIDENCE,  F-THETA //8X , 

9 1 4 .  TM  INCIDENCE,  F-PHI  ',//8X, 

1*5.  CHANGE  DATA  FILE',//8X, 

2 ' 6 .  CHANGE  ASPECT  RATIO', //8X, 

3 ' 7 .  EXIT  THE  PROGRAM ', //5X , 

4'Indicate  your  selection  by  the  choice  of  1,  2,  3,  4,  5,  6  or  7 ' , 
5//5X) 

91  format(//5X, 'EMCAD  OUTPUT  GRAPHICS  PRESENTAT I  ON ' , //5X , 

1 ' DATA  FILE  NAME  IN PUT  MENU ' , //5X, 

2'If  utilizing  a  file  or  series  of  files  created  by  EMCAD, ',/5X, 
3'please  select  option  1.  If  utilizing  any  file  other  than',/5X, 
4'one  created  by  EMCAD,  please  ensure  that  the  file  is  in',/5X, 
5'proper  format,  and  that  the  filename  extension  is  one  of',/5X, 
6'the  four  extensions  utilized  by  this  program.  These', /5X, 
/'extensions  are;  . tet  , tep  . tmt  and  .tmp',/5X, 

8'If  you  are  un*=ure  as  to  whether  your  file  is  in  the  proper',/5X, 

9 'format  with  the  proper  extension,  please  exit  the  program ',/ 5X , 
l'and  verify  before  cont inui ng . ' , // 8X , 

2 ' 1 .  INPUT  THE  FILE  OF  USER  CHOICE ' ,//8X, 

3 ' 2 .  EXIT  THE  PROGRAM ',// 5X , 

4'Indicate  your  selection  by  the  choice  of  1  or  2.'//5X) 

92  format(//5X, 'EMCAD  OUTPUT  GRAPHICS  PRESENTATION ' ,//5X, 
l'DA”A  FILE  NAME  INPUT  ',//5X, 

2'Please  input  the  name  of  the  file  containing  data  that',//5X, 
3'you  wish  to  work  with.  Please  input  only  the  file  name',//5X, 
4'i.e.  do  not  include  the  file  extension.  The  file',//5X, 
5'extension  will  be  selected  in  the  DATA  FILE  TYPE',//5X, 

6 ' SEL  ECTI ON  MENU. ',//5X) 

93  format(//5X, 'EMCAD  OUTPUT  GRAPHICS  PRESENTATION* ,//5X, 

1 ' SEL  ECTION  OF  INCIDENCE  ANGLE', /5X, 

2'The  data  fiie  you  have  selected  contains  data  for  ',I2,/5X, 
3'incidence  angles.  Please  select  your  choice  from  the',/5X, 
4'following  menu.  Note,  you  may  only  work' with  one',/5X, 
5'incidence  angle  at  a  time.',/) 

94  forma t( 5X, II , 5X, ' INCIDENCE  ANGLE  =  ',F6.2) 

95  format(5X, II , 5X, 

1 ' RETURN  TO  THE  DATA  FILE  TYPE  SELECTION  MENU' ,/5X, II , 5X, 

2 ' EXIT  THE  PROGRAM' ) 

96  format(/5X, ' INDICATE  YOUR  SELECTION  BY  THE  CHOICE  OF  A  NUMBER  I.E 

1.1  2  3  . . . ) ' ,//5X) 
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97  format(/5X, 'EMCAD  OUTPUT  GRAPHICS  PRESENTATION ', //5X, 

1 ' SELECTION  OF  PHI  ANGLE', /5X, 

2'The  data  file  you  have  selected  contains  data  for  ’>12, /5X, 

3'phi  angles.  Please  select  your  choice  from  the  f o 1 1  owl ng 1 , / 5X , 
4'menu.  Note>  you  may  only  work  with  one  phi  angle  at  a  time. 1 ,/ ) 
93  format(5X, 12, 5X, 'PHI  ANGLE  =  ' , F6 . 2 ) 

99  f o rna t ( 5X , 1 2 , 5X , 

1 ' RETURN  TO  THE  DATA  FILE  TYPE  SELECTION  MENU ' , / 5X , 1 2 , 5X , 

2 ’ EXIT  THE  PROGRAM' ) 

100  formats  A) 

101  format(/5X, 'EMCAD  OUTPUT  GRAPHICS  PRESENTATION' ,//5X, 

1 'NOTES  CONCERNING  GRAPHICAL  HARDCOPY  OUTPUT,  PART  l',//5X, 

2'This  program  allows  the  user  to  obtain  printer  hardcopies' ,/5X, 
3’of  the  graphs  being  viewed.  Hardcopy  graphic  outputs  are',/5X, 
4'obtained  through  the  use  of  the  program  entitled  GRAFPLUS . ' , /5X, 
5'Proper  utilization  of  MCADGRAF  requires  prior  setup  of',/5X, 
6'GRAFPLUS.  Directions  for  proper  setup  are  contained' ,/5X, 

7'in  the  GRAFPLUS  users  manual .' ,//5X, 

8 ' 1 .  CONTINUE  WITH  THE  NEXT  SCREEN  * ,//5X, 

9 ' 2 .  EXIT  THE  PROGRAM' , //5X, 

l'Indicate  your  selection  by  the  choice  of  1  or  2.V/5X) 

102  format(/5X,  ’ EMCAD  OUTPUT  GRAPHICS  PRESENTATION' ,//5X, 

1 'NOTES  CONCERNING  GRAPHICAL  HARDCOPY  OUTPUT,  PART  2',//5X, 

2’The  user  makes  the  decision  as  to  whether  a  printer  ',/5X, 

3’ graphical  output  is  desired  after  viewing  the  graph.  Qnce',/5X, 

4 'a  graph  is  drawn  on  the  screen,  the  user  may  request  a  ',/5X, 
S'printer  hardcopy  by  typing  either  a  "P"  or  "p"  and  pressing' ,/5X, 
6'the  "ENTER"  key.  Selection  of  any  other  key  will  simply', /5X, 
7'advance  to  the  next  graph.  Spacing  on  3  1/2  by  11  paper  ',/5X, 
8'allows  up  to  2  graphs  per  page . ' , ///5X, 

9 ' PLEASE  PRESS  ENTER  TO  CONTINUE .', 5X) 

103  format(/5X, 'xxxxxxxxxx  WELCOME  TO  MCADGRAF  XXXXXXXXXX ' , // 5X , 

1 ' THIS  PROGRAM  WAS  WRITTEN  TO  BE  USED  IN  CONJUNCTION  WITH  ',/5X, 
2'THE  FORTRAN  PROGRAM  ENTITLED  EMCAD.  MCADGRAF  ADDS  ',/5X, 

3 ' GRAPHING  FEATURES  TO  SUPPORT  THE  ANALYSIS  OF  DATA  ’,/5X, 

4 'GENERATED  AND  FORMATTED  BY  EMCAD.  EMCAD  WAS  WRITTEN  BY',/5X, 
S'DR.  M.  A.  MORGAN,  WITH  MODI FICATIONS  BY  HIMSELF  AND  ',/5X, 

6'LT  E.  M.  CONNOLLY.  MCADGRAF  WAS  DEVELOPED  AND  WRITTEN  BY',/5X, 
7’LT  E.  M.  CONNOLLY  TO  HELP  FULFILL  THESIS  REQUIREMENTS  WHILE', /5X, 
8 ' UNDER  THE  GUIDANCE  AND  THE  SUPPORT  OF  DR.  MORGAN.  MCADGRAF’ ,/5X, 
9'WAS  WRITTEN  AND  COMPILED  USING  RYAN  MCFARLAND  FORTRAN ,',/ 5X , 

1 ' GRAFMAT I C ,  AND  GRAFPLUS.  ANY  QUESTIONS  CONCERNING  COL ORS , ' , / 5X , 

2 ' FORMAT ,  CHANGES  TO  THE  PROGRAM,  ETC...  SHOULD  BE  REFERRED  ',/5X, 

3 ' TO  THE  RESPECTIVE  USERS  MANUALS,  DR.  M.  A.  MORGAN,  0R',/5X, 

4 ' L T  E.  M.  CONNOLLY. ' ,//5X, 

5 'PLEASE  PRESS  ENTER  TO  CONTINUE .', 5X) 

104  f ormat(/5X, ' XXXXXXXXXX  ASPECT  RATIO  XXXXXXXXXX *  , //5X, 

1 ' IN  ORDER  TO  AVOID  DISTORTION  ON  THE  SCREEN,  IT  IS  NECESSARY ', /5X , 
2 ' TO  DEFINE  THE  NUMBER  OF  ROWS  PER  INCH  AND  THE  NUMBER  0F',/5X, 

3 ' COLUMNS  PER  INCH.  THESE  DEFINITIONS  CHANGE  AS  THE‘,/5X, 

4 ' RESOLUTION  OF  THE  SCREEN,  AS  WELL  AS  THE  TYPE  OF  SCREEN', /5X, 

5 ' CHANGE ,  THIS  PROGRAM  WAS  WRITTEN  IN  MODE  16  (EGA)  USING', /5X, 

6 ' A  GB - 1  VIDEO  BOARD  AND  A  NEC  MULTISYNC  MONITOR.  THE',/5X, 

7 'OPTIMUM  ASPECT  RATIO  FOR  THIS  CONFIGURATION  WAS  .65', ///5X, 

8 'PLEASE  INPUT  ASPECT  RATIO . ' , //5X ) 

105  forma t( F5 . 3 ) 
end 

SUBROUTINE  RECTMVT (title2, title! , titleS , numpts , aspect ) 

“XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
MAGNITUDE  VERSUS  THETA 
REC'ANGULAR  COORDINATE  GRAPH 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx* 
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XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

VARIABLE  DECLARATIONS 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

INTEGER  A ,  B,  C,  I,  IREG.  IVAL  ,  SCALEFAC 

INTEGERS  ICOL,  IOPT,  IRON,  ISYMBL,  JC0L1,  JC0L2 ,  JR0W1 

INTEGERX2  JR0N2,  KLRSYM,  KOLOR,  LABEL,  LKOLOR ,  MINOR 

INTEGERS  MODE,  NCHAR ,  NDEC,  SELECT1 

CHARACT ER*2  SELECT,  XTENSION 

CHARACTERS  GRAFPRNT 

CHARACTERS  DUMMY,  FNAME1 ,  TITLE? 

CHARACT ER*12  FNAME2 
CHARACTER*20  CHAR 

CHARACTER*64  TITLE1,  TITLE2,  TITLE3,  TITLE4,  TITLES 
CHARACTER*64  TITLE8,  TI TL  E6 1 5 , 19,1) 

REAL  ASPECT,  CALCVAR ,  D'JMMYVAR ,  GAMAX,  GMAX,  INCANGI 
REAL  INCANGLD,  PHASEI,  PHIi,  PI,  MAGI,  NDOTS ,  NUMPTS 
REAL  NUMSETS,  RAD,  THETAI 

REAL  XAMAX ,  XFIN,  XMAJOR,  XMAX,  XMIN,  XORG,  XST 

REAL  YAMAX,  YFIN,  YMAJOR,  YMAX ,  YMIN,  YORG,  YOVERX,  YST 

REAL  AARRAYC  37 ) ,  BARRAY I  37 ) ,  INCANGL 15),  PHI ( 19 ) 

REAL  MAGI  5,19,37),  PHASEI 5 , 1 9 , 37  ) ,  THETAI 5 , 1 9 , 37  ) 

COMMON/ONE/  RMARRAY I  37 ) ,  SCALEFAC,  TARRAYI 37 ) 

COMMON/TWO/  PARRAYI 37  ) 

COMMON/THREE/  XARRAY (  37  )  ,  YARRAY I  37  ) ,  SCALEXY,  XYMAX 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SET  DISPLAY  MODE 

mode  =  16 

CALL  qsmode(mode) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
SET  AXIS  COLOR 
CALL  qpreg!7,6) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

DEFINE  GRAPH  REGION 

jcoll  =  080 

jcol2  =  520 

jrowl  =  40 

jrow2  =  320 

xmin  =  -0.5 

xmax  -  180 

ymin  -  -0.5 

ymax  -  xymax 

xorg  =  0.0 

yorg  =  0.0 

iopt  =  1 

yoverx  -  . 7 X I ymax-ymi n )/ 1 xmax-xmi n ) 

CALL  qplotIjcoll,jcol2,j  rowl ,  j  row2,  xmi  n ,  xmax,  ymin,  ymax,  xorg, 
lyorg, iopt , yoverx, aspect ) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SET  UP  LINE  TYPE 

ndots  =  0 

Ikolor  =  3 

isymbl  =  -2 

klrsym  =  4 

CALL  qsetupl ndots, Ikolor, isymbl , kl rsym ) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

DRAN  X  AXIS 

xst  =0.0 

xfin  =  180 

xmajor  =  45 

minor  =  0 

label  =  -1 

ndec  =  1 

CALL  qxaxislxst , xf in , xmaj or .minor , label, ndec) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
DRAW  Y  AXIS 
yst  =0.0 
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yfin  =  xymax 
yna:or  =  2 

lfCxynax  . eq.  1  .or.  xymax  eq.  2  .or.  xymax  .eq.  3  .or.  xymax 
l.eq.  4  .or.  xymax  . eq .  5)  ymajor  =  1 
minor  -  1 
label  =  -1 

CALL  qyaxi s  Cyst,yfin,y major, minor, label, ndec ) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WRITE  TEXT  ON  SCREEN 

nchar  -  20 

char  =  'Theta' 

kolor  =  14 

CALL  qptxtbC nchar ,  char,  kolor} 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WRITE  TEXT  ON  SCREEN 

nchar  =  20 

char  =  'Magnitude' 

kolor  =  14 

CALL  qptxtcC nchar ,  char,  kolor) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WRITE  TEXT  ON  SCREEN 

nchar  =  64 

kolor  =  10 

icol  =  15 

irow  "  24 

CALL  qptxtCnc ha r,title8, kolor, icol, irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITE  TEXT  ON  SCREEN 
nchar  -  64 

titlel  =  'Scattered  Field  Complex  Amplitudes' 
kolor  z  10 
icol  -  15 
irow  =  23 

CALL  qptxtCnc ha r, titlel, kolor,  icol,  irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WRITE  TEXT  ON  SCREEN 

nchar  -  64 

kolor  =  10 

icol  =  15 

irow  =  22 

CALL  qptxtCnc ha r , ti t le2 , kolor, icol, irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WRITE  TEXT  ON  SCREEN 

nchar  =  7 

kolor  =  10 

icol  =  15 

irow  -  21 

CALL  qptxtCnc ha r , ti t 1 e7 , kolor, icol, irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITE  TEXT  ON  SCREEN 
nchar  -  40 

title3  ~  '  Magnitude  Versus  Theta' 
kolor  =  10 
icol  =  30 
irow  ~  21 

CALL  qptxt(nchar,title3,kolor,icol .irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
SCALE  "Y"  VALUES 
do  30  i  =  1 , nunpts 

aarray(i)  =  rmarrayCi )/( 10**scalexy) 

30  continue 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
GRAPH  DATA  POINTS 
itype  =  1 
npt  -  nunpts 

CALL  qtablCi type ,npt,tarray,aarray) 


13 
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XXXXXXXXXXXXXXXXXXX  XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX* 

WRITE  SCALE  FACTOR  ON  SCREEN 

nchar  -  1 

ti  1 1  e4  =  'S' 

kolor  *  14 

icol  ~  1 

irow  =  18 

CALL  qptxtCnc ha r,title4, kolor, icol, irow) 
ti  t 1  e4  =  '  c  ' 
irow  =  17 

CALL  qptxt(nchar,title4,kolor,icol,irow) 
ti  tl e4  =  '  a  ' 
irow  =  16 

CALL  qptxt(nchar,title4,kolor,icol,irow) 
ti tle4  -  *  1 ' 
irow  -  15 

CALL  qptxt ( nchar , title4, kolor, icol, irow) 
title4  =  'e' 
irow  =  14 

CALL  qptxt(nchar,title4, kolor, icol, irow) 
ti tle4  -  'X' 
irow  =  12 

CALL  qptxt (  nchar, title4, kolor, icol, irow) 
nchar  =  2 
title4  *  '10' 
irow  =  10 

CALL  qptxtCnc ha r,title4, kolor, icol, irow) 
nchar  -  1 
title4  =  1  * 1 
irow  -  9 

CALL  qptxt(nc ha r,title4, kolor, icol, irow) 
title4  =  '  X  1 

irow  -  8 

CALL  qptxt ( nchar,title4,kolor,icol ,irow) 

CALL  qcnov ( 2,8) 

write(x,ll)  scalefac 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
HARD  COPY  OF  GRAPH  ? 
nchar  =  20 

title3  -  'Hardcopy  - >  P  or  p' 

kolor  =  14 
icol  =  1 
irow  =  1 

CALL  qp txt(nc ha r,title3, kolor, icol, irow) 

CAL  L  qcmov( 25,1) 
read(X,100)  grafprnt 
nchar  =  40 

title3  =  '  ' 

kolor  =  10 
icol  =  1 
irow  =  1 

CALL  qptxt ( nchar , title3, kolor, icol, irow) 
ifCgrafprnt  .eq.  'P'  .or.  grafprnt  .eq.  'p')  CALL  prtsc 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx* 
11  format(I3) 

100  format(A) 
r eturn 
end 

SUBROUTINE  RECTPVTC  title2,title7,title8, numpts .aspect ) 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx:.xxxxxxxxxxxxxxxxxxxxxxx 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
PHASE  ANGLE  VERSUS  THETA 
RECTANGULAR  COORDINATE  GRAPH 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
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c 

C  ****x#xxxx**xxx***x*****xxxx***x***x*xx*xxx*xxx*xxxxx**xx*x******* 

C  VARIABLE  DECLARATIONS 

C  ***x******x**xxxxxxx*x***xx*******x*****xx*x*x***xxxxxx**xx**x**** 

INTEGER  A,  B,  C,  I.  IREG.  IVAL,  3CALEFAC 
INTEGERX2  ICOL,  IOPT.  IRON,  I3YMBL,  JCOL1,  JC0L2,  JROW1 
I  NT  EGER'*  2  JR0W2,  KLRSYM,  KOLCR,  LABEL,  LKOLOR,  MINOR 
INTEGERX2  MODE,  NCHAR ,  NDEC,  SELECT1 
CHARACTERX2  SELECT,  XTENSICN 
CHARACTER*3  GRAFPRNT 
CHARACTERX8  DUMMY,  FNAME1 ,  TITLE7 
CHARACTERX12  FNAME2 
CHARACTERX20  CHAR 

CHARACTERX64  TITLE1 ,  TITLE2,  TITLE3,  TITLE4,  TITLE5 
CHARACTERX64  TITLE8,  TITL E6 C 5 , 19 , I) 

REAL  ASPECT,  CALCVAR,  DUMMYVAR,  GAMAX,  GMAX,  INCANGI 
REAl  INCANGLD,  PHASEI,  PHir,  PI,  MAGI,  NDOTS,  NUMPTS 
REAL  NUMSETS,  RAD,  THETA  I 

REAL  XAMAX ,  XFIN,  XMAJOR,  XMAX,  XMIN,  XORG,  XST 

REAL  YAMAX,  YFIN,  vMAJOR,  YMAX ,  YMIN ,  YORG,  YOVERX,  YST 

REAL  AARRAYC  37 ) ,  BARRAV f  37 ) ,  INCANGL ( 5 ) ,  PHI ( 1 9 ) 

REAL  MAG( 5,19,37),  PHASE< 5 , 19 , 37  ) ,  THETAC 5 , 1 9 , 37 ) 

COMMON/ONE/  RMARRAY ( 37  ) ,  SCALEFAC,  TARRAY ( 37 ) 

COMMON/TWO/  PARRAY ( 37  ) 

COMMON/THREE/  XARRAY ( 37  ) ,  YARRAY (  37  )  ,  SCALEXY,  XYMAX 
C  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx* 

C  SET  DISPLAY  MODE 

mode  =  16 
CALL  qsmode(mode) 

C  **x*xxx*xxxxxxxx*xxxxxxxxxxxxxxxxxxxx**x*xxx**xx*xxxxx*xxx*xxxx*xx 

C  SET  AXIS  COLOR 

CALL  qpregt 7,6) 

C  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx** 

C  DEFINE  GRAPH  REGION 

jcoll  =  030 
jcol2  =  530 
growl  -  30 
jrow2  =  320 
xmin  =  -0.5 
xmax  =  120 
ymin  =  -130 
ynax  -  ISO 
xorg  -  0.0 
yorg  =  0.0 
iopt  -  1 
yoverx  -  1.3 

CALL  qplot( jcoll, jcol2,j  rowl , j  row2 , xmin , xmax, ymin, ymax, xorg, 
lyorg, icpt, yoverx, aspect ) 

C  *xxxxx*xxxxxxxxxxxx*xxxx**xxxxxxxxxx*xxxxxxxxxxx**xxxx*xxx****xxxx 

C  SET  UP  LINE  TYPE 

ndots  =  0 
Ikolor  =  3 
isymbl  =  -2 
klrsym  =  4 

CALL  qsetup(ndots, Ikolor, i symbl , kl rsym ) 

C  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C  DRAM  X  AXIS 

xst  -  0.0 
xfin  =  180 

xmajor  =  45  t 

minor  -  0 
label  -  1 
ndec  =  1 

CALL  qxaxisC  xst , xfin , xma  jor, minor, label, ndec ) 

C  xxx*x*x?:****xxxxxx**xxxx*xxx**x*xxx**xxx*xxxxx*xx**xxx*x*xxxx*xxxx 

C  DRAW  Y  AXIS 
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yst  =  -180 
yfin  =  150 
yrajor  =  90 
minor  =  1 
label  -  -1 

CALL  qy a xis( yst, yfin, y major, minor, label,  ndec) 

X*XXX*XXXXXXXXXXXX*XXXXXXXXXXXXXXXXXX*XXXXXXXXXX*XXXXXXXXXXXXXXX** 

WRITE  TEXT  ON  SCREEN 
nchar  -  20 
char  =  'Theta ' 
kolor  =  14 

CALL  qptxtbC nchar ,  char,  kolor) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx».xxxxxxxxxxxxxxxxxx 

WRITE  TEXT  ON  SCREEN 

nchar  =  20 

char  =  'Phase' 

kolor  -  14 

CALL  qptxtcC nchar ,  char,  kolor) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WRITE  TEXT  ON  SCREEN 

nchar  =  64 

kolor  =  10 

icol  =  15 

irow  =  24 

CALL  qptxtCnc ha r,title8, kolor, icol, irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITE  TEXT  ON  SCREEN 
nchar  -  64 

title)  -  'Scattered  Field  Complex  Amplitudes' 
kolor  -  10 
icol  =  15 
irow  =  23 

CALL  qptxtCnchar.titlel, kolor, icol, irow) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXxXXXXXXXXXXXXXXXXXXXXXXXX 

WRITE  TEXT  ON  SCREEN 

nchar  =  64 

kolor  1  10 

icol  1  15 

irow  =  22 

CALL  qp txtCnc ha r,title2, kolor, icol, irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WRITE  TEXT  ON  SCREEN 

nchar  -  7 

kolor  -  10 

icol  =  15 

irow  =  21 

CALL  qptxt(nchar, title?, kolor, icol, irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITE  TEXT  ON  SCREEN 
nchar  =  40 

title3  -  '  Phase  Versus  Theta' 
kolor  =  10 
icol  =  30 
irow  =  21 

CALL  qptxtCnc ha r, title!, kolor, icol , irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx* 
GRAPH  DATA  POINTS 
itype  -  1 
npt  =  numpts 

CALL  qtabl ( itype ,npt,tarray,parray) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX,';XXXXXXXXXXXXXXXXXXXXXXX 

HARD  COPY  OF  GRAPH  ? 
nchar  =  20 

title3  -  'Hardcopy  - >  P  or  p' 

kolor  =  14 
icol  =  40 
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irow  =  1 

CALL  qptxt(nchar,title3,kolor,icol,irow) 

CALL  qcmov (65,1) 
read(x,100)  grafprnt 
nchaf  -  40 

ti tle3  =  '  ' 

kol o r  -  10 
i c  o 1  =  4  0 
l row  -  1 

CALL  qptxt(nchar,title3,kolor,icol,irow) 
if(grafprnt  . eq .  'P'  .or.  grafprnt  . aq.  'p')  CALL  prise 

C  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

100  formatCA.' 
return 
end 

SUBROUTINE  POL ARMV T( ti t 1 e2 , ti t 1 e7 , ti tl e8 , numpts , aspect ) 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
MAGNITUDE  VERSUS  THETA 
POLAR  COORDINATE  GRAPH 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
VARIABLE  DECLARATIONS 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXxXXXXXXXXXXXXXXXWXXXXXXXXXXXXXXXXXX 

INTEGER  A,  B,  C,  I,  IREG,  IVAL,  SCAI.EFAC 
INTEGERX2  ICOL,  I0PT,  IROW,  ISYMBL,  JC0L1,  JC0L2,  JR0W1 
INTEGERX2  JR0W2,  KLRSYM,  K0L0R,  LABEL,  LKOLOR,  MINOR 
INTEGERX2  MODE,  NCHAR ,  NDEC ,  SELECT1 
CHARACTERX2  SELECT,  XTENSION 
CHARACTERX3  GRAFPRNT 
CHARACTERX8  DUMMY,  FNAME1 ,  TITLE? 

CHARACTERX12  FNAME2 
CHARACTERX20  CHAR 

CHARACTERX64  TITLE1 ,  TITLE2,  TITLE3,  TITLE4,  TITLE5 
CHARACTERX64  TITLE8 ,  TITLES (5, 19 , 1 ) 

REAL  ASPECT,  CALCVAR ,  DUMMYVAR ,  GAMAX,  GMAX,  INCANGI 
REAL  INCANGLD,  PHASEI ,  PHII,  PI,  MAGI,  NDOTS,  NUMPTS 
REAL  NUMSETS,  RAD,  THETAI 

REAL  XAMAX,  XFIN,  XMAJOR,  XMAX,  XMIN,  XORG,  XST 
REAL  YAMAX,  YFIN ,  YMAJOR ,  YMAX ,  YMIN,  YORG,  YOVERX,  YST 
REAL  AARRAY ( 37  )  ,  BARRAY ( 37  )  ,  INCANGL ( 5 ) ,  PHIC19) 

REAL  MAGC5, 19, 37 ) ,  PHASEI 5,19,37),  THET A C 5 , 1 9 , 37 ) 

COMMON/ONE/  RMARRAY ( 37  )  .  SCALEFAC,  TARRAY ( 37  ) 

COMMON/TWO/  PARRAY ( 37 ) 

COMMON/THREE/  XARRAY ' 37 ) ,  YARRAY ( 37 ) ,  SCALEXY,  XYMAX 
xxx*xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
SET  DISPLAY  MODE 
mode  =  16 
CALL  qsmode(mode) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
SET  AXIS  COLOR 
CALL  qpregt 7,6) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
SCALE  "X"  AND  "Y"  VALUES 
do  10  i  =  1, numpts 

aarray(i)  -  xa rray < i )/ ( 1 GXXsca 1 exy ) 
barray(i)  =  ya  may  ( l )/ ( 1 0  xxsca  1  exy ) 

10  continue 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


DEFINE 

GRAPH  REGION 

i  c  o  1 1 

=  080 

j  c  o  1 2 

-  400 

j  r  owl 

=  10 
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jrow2  =  280 
xmin  =  -xymax/2 
xnax  -  xymax 
ymn  =  -xynax 
ynax  =  xynax 
xorg  =  0.0 
yorg  -  0.0 
iopt  =  1 

yoverx  -  0.55 

CALL  qplot(jcoll,jcol2,jrowl,j row2 , xmi n , xnax, ymin.y max, xorg, 

1 yorg, iopt, yoverx, aspect) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX* 

SET  UP  LINE  TYPE 
ndots  =  0 
lkolor  =  3 
isymbl  =  -2 
klrsym  =  4 

CALL  qsetup(ndots, lkolor, i symbl , k 1 rsym ) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

DRAW  X  AXIS 

xst  =0.0 

xfin  =  xymax 

xmajor  =  2 

ifCxymax  le.  5)  xmajor  =  1 
minor  =  1 
label  =  1 
ndec  =  1 

CALL  qxaxi s( xst , xf i n , xmajor, mi  nor, label, ndec) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

DRAW  Y  AXIS 

yst  =  -xymax 

yfin  =  xymax 

ynajor  =  xmajor 

label  =  Cl 

CALL  qyaxisCyst, yfin, yma j or, minor, label ,  ndec ) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx** 

WRITE  TtXT  ON  SCREEN 

nchar  =  10 

char  =  'Theta=90' 

kolor  =  14 

CALL  qptxtbC  nchar,  char,  kolor) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WRITE  TEXT  ON  SCREEN 

nchar  =  10 

char  =  ' Theta=Q ' 

kolor  =  14 

i col  =  10 

irow  =  19 

CALL  qptxt (  nchar,  char,  kolor,  icol,  irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WRITE  TEXT  ON  SCREEN 

nchar  =  10 

char  =  'Theta=180' 

kolor  =  14 

i col  =  10 

irow  =  1 

CALL  qptxtC  nchar,  char,  kolor,  icol,  irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WRITE  TEXT  ON  SCREEN 

nchar  =  64 

kolor  =  10 

icol  =  15 

i  roi  •  =  24 

CALL  qptxtC  nchar,  title8,  kolor,  icol,  irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WRITE  TEXT  ON  SCREEN 
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Scattered  Field  Complex  Amplitudes 


nchar  1  64 
titlel  =  ' 
kolor  =  10 
icol  -  15 
irow  -  23 

CALL  qptxtC  nchar,  titlel,  kolor,  icol,  irow) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

HR  ITS  TEXT  ON  SCREEN 

nchar  =  64 

kolor  -  10 

icol  =  15 

irow  =  22 

CALL  qptxt(nchar,title2, kolor, icol, irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITE  TEXT  ON  SCREEN 
nchar  -  11 

title3  =  'Polar  Plot,  * 
kolor  =  10 
icol  =  15 
irow  -  21 

CALL  qptxt(nchar,title3, kolor, icol, irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WRITE  TEXT  ON  SCREEN 

nchar  -  7 

kolor  =  10 

icol  =  27 

irow  =  21 

CALL  qp t xt (nc ha r,title7, kolor, icol, irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITE  TEXT  ON  SCREEN 
nchar  -  40 

titleS  -  '  Magnitude  Versus  Theta' 
kolor  =  10 
icol  -  35 
irow  =  21 

CALL  qptxtC nchar, title3, kolor, icol, irow) 

xxxx*xxxx**xx*xxxxx*x*xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WRITE  SCALE  FACTOR  ON  SCREEN 

nchar  =  1 

title4  =  'S' 

kolor  =  14 

icol  s  10 

irow  =  17 

CALL  qp t xt (nc ha r,title4, kolor, icol, irow) 
ti tl e4  =  'c' 
irow  =  16 

CALL  qptxtC  nchar , titled, kolor, icol, i row ) 
title4  -  'a' 
irow  =  15 

CALL  qptxt(nc ha r,title4, kolor, icol, irow) 
titie4  =  '  1  ' 
irow  =  14 

CALL  qptxt( nchar , titl e4 , kolor , : ^  >l,irow) 
title4  =  'e' 
irow  =  13 

CALL  qptxt(nchar,title4, kolor, icol, irow) 
title4  =  'X* 
irow  =  11 

CALL  qptxtCnchar, title4, kolor, icol, i row) 
nchar  =  2 
title4  =  '10' 
irow  -  9 

CALL  qptxtC  nchar , tit  1 e4, kolor , icol , irow) 
nchar  =  1 
ti t 1 e4  = 
irow  =  8 
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CALI  qptxt?  nchar , titled, kolor,icol,i row) 
title4  =  1  *  1 
irow  -  7 

CALL  qptxt?nc ha r,title4, kolor, icol, irow) 

CALL  qcmov ( 5,7) 
write?*, 12)  scalefac 

See  discussion  of  QCMOV  for  problem  requiring  this  fix 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

GRAPH  DATA  POINTS 
itype  =  1 
npt  =  numpts 

CALL  qtabl ? i type, npt,aarray,barray) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
HARD  COPY  OF  GRAPH  ? 
nchar  =  20 

title3  =  'Hardcopy  - >  P  or  p' 

kolor  =  14 
icol  -  50 
irow  »  1 

CALL  qptxt ( nchar,title3,kolor,icol,irow) 

CALL  qcmov (75,1) 
read(*,10Q)  grafprnt 
nchar  =  26 

title3  =  '  ' 

kolor  =  10 
icol  -  50 
irow  =  1 

CALL  qptxt?nchar,title3, kolor, icol, irow) 
ifCgrafprnt  .eq.  'P'  .or.  grafprnt  . eq .  'p')  CALL  prtsc 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
12  forma t ( 9x , 13) 

100  formats  A ) 
return 
end 

SUBROUTINE  POLARROT? title2 , ti tle7 , title8 , numpts, aspect ) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

MAGNITUDE  VERSUS  THETA 

POLAR  COORDINATE  GRAPH 

AXIS  ROTATED  MINUS  NINETY  DEGREES 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

xxxxxXxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
VARIABLE  DECLARATIONS 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

INTEGER  A,  B,  C,  I,  IREG,  IVAL,  SCALEFAC 

INTEGERS  ICOL,  IOPT,  IROW,  ISYMBL,  JC0L1,  JC0L2,  JR0W1 

INTEGERX2  JR0W2,  KLR3YM,  KOLOR,  LABEL,  LKOLOR,  MINOR 

INTEGERS  MODE,  NCHAR,  NDEC ,  SELECT1 

CHARACTERS  SELECT,  XTENSI0N 

CHARACTERX3  GRAFPRNT 

CHARACTER*8  DUMMY,  FNAME1 ,  TITLE7 

CHARACTERX12  FNAME2 

CHARACTERX20  CHAR 

CHARACTERX64  TITLE1,  TITLE2,  TITLE3,  TITLE4,  TITLE5 
CHARACTERX64  TITLE8,  TITL E6 ( 5 , 19 , 1  ) 

REAL  ASPECT,  CALCVAR ,  DUMMYVAR ,  GAMAX,  GMAX,  INCANGI 
REAL  INCANGL D,  PHASEI,  PHII,  PI,  MAGI,  NDOTS,  NUMPTS 
REAL  NUMSETS,  RAD,  THETAI 

REAL  XAMAX,  XFIN,  XMAJ0R,  XMAX,  XMIN,  X0RG,  XST 

REAL  YAMAX,  YFIN,  YMAJOR,  YMAX,  YMIN,  YORG,  YOVERX,  YST 

REAL  AARRAY ?  37  ) ,  BARRAY ( 37  ) ,  INCANGL ( 5  ) ,  PHI ( 1 9  ) 

REAL  MAG?  5,19,37),  PHASE?  5,19,37),  THETA? 5 , 19 , 37 ) 

COMMON/ONE/  RMARRAY?  37 ) ,  SCALEFAC,  T  ARRAY  ?  37 ) 
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COMMON/TWO/  PARRAY ( 37 ) 

COMMON/THREE/  XARRAY C  3 7  )  ,  YARRAYC  37  ) ,  SCALEXY,  XYMAX 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

SET  DISPLAY  MODE 

node  =  16 

CALI  qsmode(mode) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx******** 
SET  AXIS  COLOR 
CALL  qpreg( 7.6) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

DEFINE  GRAPH  REGION 

jcoll  =  080 

jcol2  =  530 

jrowl  -  10 

jrow2  =  320 

xmin  =  -xymax 

xmax  =  xymax 

ymin  =  -xymax/4 

ymax  =  xymax 

xorg  =  0.0 

yorg  =  0.0 

iopt  =  1 

yoverx  =  0.5 

CALL  qplot(jcoll>jcol2<j rowl , jrow2, xmi n , xmax, ymin, ymax, xorg, 
lyorg, iopt, yoverx, aspect) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SET  UP  LINE  TYPE 

ndots  -  0 

Ikolor  =  3 

isymbl  =  -2 

klrsym  *  4 

CALL  qsetup( ndots, Ikolor, isymbl , kl rsym) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

DRAW  X  AXIS 

xst  =  -xymax 

xfin  ~  xymax 

xmajor  -  2 

i f (xymax  .le.  5)  xmajor  -  1 
minor  =  1 
label  =  0 
ndec  =  1 

CALL  qxaxisCxst , xf in, xma j or , minor, label ,  ndec ) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

DRAW  Y  AXIS 

yst  =  0.0 

yfin  -  xymax 

ymajor  =  xmajor 

label  =  1 

CALL  qyaxis (yst, yfin.ymaj or .minor , label .ndec) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

WRITE  TEXT  ON  SCREEN 

nchar  =  10 

char  -  'Theta^SO' 

kolor  =  14 

icol  =  65 

irow  =  2 

CALL  qptxt( nchar , char. kolor, icol .irow) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

WRITE  TEXT  ON  SCREEN 
nchar  -  10 
char  =  'Theta=90' 
kolor  =  14 

CALL  qptxtc(  nchar,  char,  kolor) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITE  TEXT  ON  SCREEN 
nchar  -  10 


189 


on  on  oo  on  on  oo  oo  oo 


char  -  'Theta  =  O' 
kolor  -  14 
icol  3  6 
irow  =  2 

CALL  qp txt(nc ha r, char, kolor, icol, irow) 

C  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C  WRITE  TEXT  ON  SCREEN 

nchar  =  64 
kolor  3  10 
icol  *  15 
irow  -  24 

CALL  qptxt(  nchar,  title8,  kolor,  icol,  irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITE  TEXT  ON  SCREEN 
nchar  3  64 

titlel  3  'Scattered  Field  Complex  Amplitudes' 
kolor  -  10 
icol  *  15 
irow  3  23 

CALL  qptxt(  nchar,  titlel,  kolor,  icol,  irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITE  TEXT  ON  SCREEN 
nchar  3  64 
kolor  3  10 
icol  3  15 
irow  3  22 

CALL  qptxt(nchar,title2, kolor, icol, irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITE  TEXT  ON  SCREEN 
nchar  3  11 

title3  3  'Polar  Plot,' 
kolor  3  10 
icol  3  15 
irow  3  21 

CALL  qptxt(nchar,title3, kolor, icol, irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITE  TEXT  ON  SCREEN 
nchar  3  7 
kolor  3  10 
icol  3  27 
irow  3  21 

CALL  qptxt(nchar,title7, kolor, icol, irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITE  TEXT  ON  SCREEN 
nchar  3  40 

title3  3  '  Magnitude  Versus  Theta' 
kolor  3  10 
icol  3  35 
irow  3  21 

CALL  qptxt(nchar,title3, kolor, icol, irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITE  TEXT  ON  SCREEN 
nchar  3  64 

title5  3  'Axis  Rotated  Minus  Ninety  Degrees' 
kolor  3  10 
icol  3  15 
irow  3  20 

CALL  qptxtC  nchar , title5, kolor, icol, irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
SCALE  "X"  AND  "Y"  VALUES 
ao  20  i  ‘  1  numpts 

aarray(i)  3  -yarray( i )/ ( lGxXscalexy) 
barray(i)  3  xarrayi i )/ ( 1 0**scal exy ) 

20  continue 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
GRAPH  DATA  POINTS 
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itype  -  1 
npt  =  numpts 

CALL  qtabl (  itype,  npt ,  aarray ,  ba may  ) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  XXXXXXXXXXXXXXXXXXXXXXXX 

WRITE  SCALE  FACTOR  ON  SCREEN 
nchar  =  1 
title4  =  'S’ 
kolor  =  14 
icol  =  5 
irow  -  16 

CALL  qptxt(nchar, titl e4, kolor, icol, irow) 
title4  =  'c' 
irow  =  15 

CALL  qptxt( nchar, title4, kolor, icol, irow) 
title4  =  'a' 
irow  -  14 

CALL  qptxt( nchar , title4,kolor, icol, irow) 
title4  =  '1' 
irow  r  13 

CALL  qptxt ( nchar, title4,kolor, icol, irow) 
title4  =  'e' 
irow  =  12 

CALL  qptxt (nchar, title4,kolor, icol, irow) 
ti tle4  =  'X' 
irow  =  10 

CALL  qptxt ( nchar , title4,kolor, icol, irow) 
nchar  -  2 
titl e4  =  «10' 
irow  =  8 

CALL  qptxt<  nchar , title4,kolor,icol,i row) 
nchar  -  1 
ti tl e4  =  '*< 
irow  -  7 

CALL  qptxt ( nchar, titl e4, kolor, icol, irow) 
title4  =  'x* 
irow  ~  6 

CALL  qptxt (nchar, titl e4,kol or, icol , irow) 

Call  QCMOV (4,6) 
wri te( X , 1 3 )  scalefac 

See  discussion  of  QCMOV  for  problem  requiring  this  fix 

X**SXX*XXXXXXXXXXX*XXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXX*XX* 

HARD  COPY  OF  GRAPH  ? 
nchar  =  20 

title3  =  'Hardcopy  - >  P  or  p' 

kolor  -  14 
icol  -  30 
irow  =  1 

CALL  qptxt ( nchar , title3, kolor, icol, irow) 

CALL  qcmov(55,l) 
read(x,lQ0)  grafprnt 
nchar  =  40 

title3  =  '  ' 

kolor  =  10 
icol  =  30 
irow  =  1 

CALL  qptxt ( nc ha r,title3, kolor, icol, irow) 
if(grafprnt  ,eq.  'P'  .or.  grafprnt  .eq.  'p')  CALL  prtsc 
C  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
13  format(4x, 13) 

100  format(A) 
return 
end 
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APPENDIX  J 


SPHERE  SOURCE  CODE 


PROGRAM  SPHERE 

GENERATING  THE  RADII  8  Er,  Ur  AT  ROUS  I  8  1+1 
FOR  A  LAYERED  DIELECTRIC  SPHERE  HAVING 

Z-AXIS  OFFSET/LAMBDA  =  D.  Written  By  M.  Morgan  2/14/88 

Modified  By  E.  Connolly  4/3/88 

REAL  RAD( 5 ) , RO! 5  ) 

COMPLEX  ER , UR 

CHARACTERX64  FNAME, HDR1 , HDR2 
PI=3. 1415927 
WRITE!*, 110) 

READ!*, 100)  FNAME 

OPEN! 1 , FI L E=FNAME, STATUS=' UNKNOWN ' ,ACCESS= 'SEQUENTIAL ' , 

1  F0RM=' FORMATTED' ) 

WRITE!*, Ill) 

READ!*, 100)  FNAME 

OPEN! 2 , FI  LE  =  FNAME, STATUS3 ' UNKNOWN' , ACCESS* ’ SEQUENTIAL  ' , 

1  FORM3 ' FORMATTED' ) 

WRITE!*, 112) 

READ!*, 100)  HDR1 
WRITE!*, 113) 

READ!*, 100)  HDR2 
WRITE!1,105)  HDR1 
WRITE! 1,105)  HDR2 
WRITE!*, 114) 

READ!*,*)  NL 
WRITE! 1,101)  NL 
DO  11  L  =  1 , NL 

WRITE!*,*)  'Layer  No:  *,L 
WRITE!*,*) 

WRITE!*,*)  'Enter  Radius/Lambda :  ' 

READ!*,*)  RO ! L  ) 

WRITE!*,*)  'Enter  Real  Er  of  This  Layer:  ' 

READ!*,*)  ERR 

WRITE!*,*)  'Enter  Imag  Er  of  This  Layer:  ' 

READ!*,*)  ERI 
ER=CMPLX! ERR, ERI ) 

UR=GMPLX!1.0,-1 . E-6 ) 

WRI i E! 1 , 1 02 )  ER, UR 
11  CONTINUE 

WRITE!*,*)  'Enter  Z-Axis  Offset/Lambda:  * 

READ!*,*)  D 

WRITE!*,*)  'Enter  IBIG:  ' 

READ!*,*)  IBIG 
WRITE! 2, 104)  IBIG,NL+1 
RMIN=2.*PI*!R0! 1)  -D) 

RMAX=2 . *PI*!R0!NL )+D) 

WRITE! 2 , 102)  RMIN, RMAX 
DT=PI/! IBIG-1 . 0) 

DO  33  1=1, IBIG 
T  =  ! I — 1 . 0 ) *DT 
DS2=!D*SIN!T))**2 
C=C0S!T) 

DO  22  L  =  1 ,  NL 

22  RAD!L)=2.0*PI*! D*C+SQRT! RO ! L ) **2-DS2 ) ) 

WRITE! 2,103)  I , ! RAD! L ) , l =1 , NL  ) 

100  FORMAT ! A ) 
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101  FORMAT  (15) 

102  FORMAT (4(E14.6)) 

100  FORMAT (I5,5(E14.6)) 

104  FORMAT (2(15)  ) 

ig:  formatc '  ’,a) 

110  FORMAT(//7X, 'TWO  FILES  ARE  REQUIRED  TO  HOLD  PROGRAM  OUPUT ' ,//7X, 
i ' DATA .  PLEASE  ENTER  THE  NAME  OF  THE  FILE  TO  BE  USED  TO  ',//7X, 
2'HOLD  THE  PROGRAM  MATERIAL  PARAMETER  DATA.  PLEASE  ',//7X, 
3'INCLUDE  THE  EXTENSION.  I.  E.  FIL ENAME . DAT ' ,//7X) 

111  F0RMATC//7X, 'PLEASE  ENTER  THE  NAME  OF  THE  FILE  TO  BE  USED  ',//7X, 
1 ' TO  HOLD  THE  PROGRAM  OUTPUT  DATA.  PtFASE  INCLUDE  THE',//7X, 

2 ' EXTENSION .  I.  E.  FILENAME. DAT' ,//7X) 

112  FQRMATC//7X, 'THE  HEADERS  ALLOW  THE  USER  TO  IDENTIFY  THIS',//7X, 

1 ' SET  OF  DATA  FROM  ALL  OTHER  SETS.', 

2//7X, 'PLEASE  ENTER  HEADER  #1  (64  CHARACTERS  MAX)',//7X) 

113  FORMAT(//7X, 'PLEASE  ENTER  HEADER  #2  (64  CHARACTERS  MAX)',//7X) 

114  FORMAT(//7X, 'ENTER  THE  NUMBER  OF  LAYERS  (.LE.  5)',//7X) 

03  CONTINUE 

CLOSEd) 

CL  OSE( 2 ) 

STOP 

END 
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APPENDIX  K 


CONE  SOURCE  CODE 


PROGRAM  CONE 

XX*XX*XXXXXXXXX*XXX*X**XXXXXXXX*X**X*XXXXXXX*XXXXXXXXX*XXXXXXX*XXX 

GENERATING  THE  RADII  S  Er,  Ur  AT  ROWS  I  8  1+1 
FOR  A  LAYERED  DIELECTRIC  CONE 
Written  3y  E.  Connolly  4/19/88 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
VARIABLE  DECLARATIONS 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  RAD(1000,5),  PI,  DELTHETA,  DEGTORAD 

REAL  THETA,  HC  5 ) ,  X!5),  YC  5  5 

REAL  ALPHAK5),  ALPHA!  5) ,  B,  M,  RL, GAMMA 

INTEGER  I  BIG ,  I 

COMPLEX  ER, UR 

CHARACTERS  FNAME,  HDR1 ,  HDR2 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
INITIAL  VALUES 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

RMIN  =  10000 

RMAX  =  -10000 

PI =3 . 1415927 

DEGTORAD  =  PI/180. 

x***xxx*x*xxxxx**xx*xxx**xxx*xxxx***xx**xxxxxxx**x*x**x**x**xx**x* 
INPUT  NECESSARY  DATA 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITE! X, 110) 

READ!*, 100)  FNAME 

0PEN(1, FIL E= FNAME, STATUS* ’ UNKNOWN*, ACCESS=* SEQUENTIAL ' , 

1  FORM^ ' FORMATTED ' ) 

WRITE(x.lll) 

READ!*, 100)  FNAME 

OPEN! 2 , FI LE= FNAME, STATUS = ' UNKNOWN ' , ACCESS = ’SEQUENTIAL ’ , 

1  F0RM= ' FORMATTED ' ) 

WRITE!*, 112) 

READ( * , 100 )  HDR1 
WRITE!*, 113) 

READ! X , 1 00  )  HDR2 
WRITE(1,105)  HDR1 
WRITE! 1,105)  HDR2 
WRITE!*, 114) 

READ!*,*)  NL 
WRITE! 1,101)  NL 
DO  11  L  =  1 ,  NL 

WRITE!*,*)  '  Layer  No:  ',  L 
WRITE!*,*) 

WRITE!*, 10) 

READ!*,*)  ALPHA1U) 

ALPHA! L )  =  ! ALPHA1 ! L )/2 . )  *  DEGTORAD 
WRITE!*,*) 

WRI , E!*, 15) 

READ!*,*)  H! L  ) 

WRITE!*,*)  '  Enter  Real  Er  of  This  Layer:  ' 

READ!*,*)  ERR 

WRITE!*,*)  '  Enter  Imag  Er  of  This  Layer:  ' 
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READ(x,X)  ERI 
ER -CMPLX( ERR , ERI ) 

UR=CMPLX(1 .0,-1 . E-6 ) 

WRITEC1.102)  ER ,  UR 
11  CONTINUE 

WRITE(x.X)  'Enter  IBIG:  ' 

READ! * , x  )  IBIG 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

CALCULATIONS 

XXX*XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX*XXXXXXXX*XXXXXXX**XXXXXXX**XX 

DELTHETA  =  PI/(IBIG-1) 

DO  22  1=1, IBIG 

THETA  =  (1-1)  X  DELTHETA 

DO  33  L  =  1 , NL 

Y(L)  =  H(L)/2 

X  (  L  )  =  (Y(L))XTAN(ALPHA(L)) 

CHECK  FOR  HORIZONTAL  LINE 
CAMMA  =  ATAN(X(L)/Y(L)  ) 

IF! (PI-GAMMA-THETA)  .LT.  .0001)  THEN 
YFOUND  =  -YCL) 

XFOUND  =  YFOUND  x  TAN(  THETA  ) 

RAD( I , L '  =  2XPIXSQRTCXF0UNDXX2  +  YF0UND*x2 ) 

ELSE 

SLOPE  OF  LINE  CONTAINING  CORNERS  OF  CONE 
M  =  (  ~H(L)/XCL)  ) 

Y-INTERCEPT  OF  LINE  CONTAINING  CORNERS  OF  CONE 
B  =  v( L  ) 

RECALL  EQUATION  OF  ANY  LINE  IS  Y  =  MX  +  B 

THE  LINE  CONTAINING  THE  I-TH  AND  K-TH  DATA  POINTS  INTERSECTS 
THE  LINE  DRAWN  FROM  THE  ORIGIN  THROUGH  THE  DESIRED  DATA  POINT 
AT  THE  DESIRED  DATA  POINT  XFOUND  AND  YFOUND.  THE  EQUATION  FOR 
THE  LINE  CONTAINING  THE  I-TH  AND  K-TH  DATA  POINTS  MUST  THEN 
ALSO  WORK  AT  THE  DESIRED  DATA  POINT.  WRITING  XFOUND  AND 
YFOUND  IN  TERMS  OF  THEIR  ANGLE  THETACN),  AND  SUBSTITUTING 
YFOUND  =  B  /  (  1  -  (  M  x  TAN(  THETA  )  )  ) 

XFOUND  =  YFOUND  X  TAN(  THETA  ) 

RAD( I , L )  =  2XPIXSQRTCXF0UNDXX2  +  YFOUNDXX2 ) 

ENDIF 

IF(RAD(I,L)  . GT .  RMAX )  RMAX  =  RADII, L) 

I F( RAD( I , L )  .LT.  RMIN )  RMIN  =  RADCI-L) 

33  CONTINUE 
22  CONTINUE 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

OUTPUT 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITE! 2 , 104 )  IBIG.NL+l 
WR I  TEC  2 , 102)  RMIN,  RMAX 
DO  44  I  =  1 , IBIG 

WRITE! 2, 103)  I , ( RAD! I , L ) , L = 1 , NL ) 

44  CONTINUE 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
10  F0RMATC/8X, ' Please  enter  the  apex  angle,  ALPHA1,  in  degrees,  of  ’ 
l/3X,'the  offset  cone.  Note,  apex  angle  is  the  angle  from  one  side' 
2/8X,'of  the  cone  to  the  other  side  of  the  cone.  ') 

15  F0RMAK/8X,  '  Please  enter  the  height,  H,  of  the  '/ 

13X , '  cone .  '  ) 

100  FORMAT ( A ) 

101  FORMAT ( 15 ) 

102  FORMAT ( 4! E14 . 6  ) ) 

103  F0RMAT(I5,5(E14.6)) 

104  FORMAT (2(15)) 

105  FORMAT ( '  ' ,A) 

110  F0RMATC//7X, 'TWO  FILES  ARE  REQUIRED  TO  HOLD  PROGRAM  OUPUT',//7X, 
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1 ' DATA ,  PLEASE  ENTER  THE  NAME  OF  THE  FILE  TO  BE  USED  TO  ' ,//7X, 
2'HOLD  THE  PROGRAM  MATERIAL  PARAMETER  DATA.  PLEASE  './/IX, 
3'INCLUDE  THE  EXTENSION.  I.  E.  FILENAME. DAT  1  ,/v7X) 

111  FORMATC//7X, 'PLEASE  ENTER  THE  NAME  OF  THE  FILE  TO  BE  USED  './/IX, 
1 1  TO  HOLD  THE  PROGRAM  OUTPUT  DATA.  PLEASE  INCLUDE  THE‘,//7X. 

2 ’ EXTENSION .  I.  E.  FILENAME . DAT* ,//7X) 

112  FORMAT { / /IX ,  ' THE  HEADERS  ALLOW  THE  USER  TO  IDENTIFY  THIS',//7X, 

1 1  SET  OF  DATA  FROM  ALL  OTHER  SETS.’, 

2//7X, ’PLEASE  ENTER  HEADER  #1  (64  CHARACTERS  MAX)',//7X) 

113  FORMAU//7X, 'PLEASE  ENTER  HEADER  *2  (64  CHARACTERS  MAX)’,//?X) 

114  FORMATC//7X, 'ENTER  THE  NUMBER  OF  LAYERS  (.L£.  5)’,//7X) 

CLOSE(l) 

CL  OSE( 2 ) 

STOP 

END 
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APPENDIX  L 

CYLINDER  SOURCE  CODE 


PROGRAM  CYLINDER 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX* 

GENERATING  THE  RADII  &  Er,  Ur  AT  ROWS  I  &  1+1 
FOR  A  LAYERED  DIELECTRIC  CYLINDER 
Written  By  E.  Connolly  4/19/88 

xxxxxxxxxxxxxxxxxx^xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
VARIABLE  DECLARATIONS 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
REAL  RADC 1000,5) .  PI,  DELTHETA,  DEGTORAD 
REAL  THETA,  HC5  ,  RC5),  YC5) 

REAL  RL,  GAMMA 
INTEGER  I  BIG,  1 
COMPLEX  ER, UR 

CHARACTERX64  FNAME, HDR1 , HDR2 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
INITIAL  VALUES 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

RMIN  =  10000 

RMAX  =  10000 

PI=3. 1415927 

DEGTORAD  =  PI/180. 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
INPUT  NECESSARY  DATA 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITEC  X , 1 1 0  ) 

READC*. 100)  FNAME 

OPENC  1 ,  FI L E- FNAME,  STATUS = 'UNKNOWN • , ACCESS= • SEQUENTIAL • , 

1  FORM= ' FORMATTED ' ) 

WRITECX, 111 ) 

READC*. 100)  FNAME 

OPENC2, FI LE  =  FNAME,STATUS= 'UNKNOWN '  ,ACCESS=' SEQUENTIAL  • , 

1  F0RM= ' FORMATTED ' ) 

WRITEC*, 112) 

READC  X, 100  )  HDR1 
WRITEC*, 113) 

READC  X, 1 00 )  HDR2 
WRITEC 1 , 1 05 )  HDR1 
WRITEC 1,105)  HDR2 
WRITECX, 114) 

READC  X , X )  NL 
WRITECl.lOl)  NL 
DO  11  L  =  1 , NL 

WRITECX, X)  'Layer  No :  ' ,  L 
WRITEC  X, X ) 

WRITECX, 15) 

READC  x , x )  HCL) 

WRITECX, 16) 

READC  * , x  )  RCL) 

WRITECX, x)  'Enter  Real  Er  of  This  Layer:  ' 

READC*,*)  ERR 

WRITEC*,*)  'Enter  Imag  Er  of  This  Layer:  ' 

READC*,*)  ERI 
ER=CMPLXC  ERR, ERI ) 
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UR=CMPLX( 1 .0,-1 . E-6 ) 

WR1TEC 1,102)  ER,UR 
11  CONTINUE 

WRITEC*/*)  'Enter  IBIG:  ' 

READC  * , *  )  IBIG 

X********''******************************************************** 

CALCULATIONS 

*XXX*X*X-.*X**XXX**X***X*XXXXXX*XXXXX*X*X**X**XXXXX*XXXX*XXX****XXX 
DELTHETA  =  PI/CIBIG-1) 

DO  22  1=1, IBIG 

THETA  =  (1-1)  *  DELTHETA 

DO  33  L  =  1 , NL 

YCL)  =  H(L)/2 

CHECK  FOR  HORIZONTAL  LINE 

GAMMA  =  ATANC  R( L )/Y( L )  ) 

I F( ( PI -GAMMA -THETA)  .LT.  .0001)  THEN 
YFOUND  =  -YCL) 

XFOUND  =  YFOUND  *  TAN(  THETA  ) 

RAD( I , L )  =  2*PI*SQRT(XF0UND**2  +  YF0UND**2) 

CHECK  FOR  HORIZONTAL  LINE 

ELSE  IFC (GAMMA-THETA)  . GT .  0)  THEN 

YFOUND  =  YCL) 

XFOUND  =  YFOUND  *  TAN(  THETA  ) 

RADCI/L)  =  2XPIXSQRT (XFOUND* *2  +  YF0UND**2) 

ELSE 

CHECK  FOR  VERTICAL  LINE 
XFOUND  =  RCl ) 

YFOUND  =  XFOUND  /  TAM(  THETA  ) 

RAD(I,L)  =  2*PI*SQRT(XF0UND**2  +  YF0UND**2) 

ENDIF 

IFC RAD(I,L)  .GT.  RMAX )  RMAX  =  RADC I , L ) 

IFC RAD(I,L)  .LT.  RMIN)  RMIN  =  RADCI.L) 

33  CONTINUE 
22  CONTINUE 

****************************************************************** 

OUTPUT 

****************************************************************** 
WRITEC2, 104)  IBIG,NL+1 
WRITE(2,102)  RMIN,  RMAX 
DO  44  I  =  1 , IBIG 

WRITEC2,103)  I, (RAD(I,L), L=1,NL) 

44  CONTINUE 

************************ ********************* ********************* 

15  FJRMATC/8X, ’Please  enter  the  height,  H,  of  the  V 

18X,'  cylinder  (from  base  to  base).  ') 

16  FORMATC/8X, 'Please  enter  the  radius  R,  of  the  V 

18X, '  cylinder .  '  ) 

100  FORMAT ( A ) 

101  FORMAT ( 15 ) 

102  FORMAT ( 4( E14 . 6 ) ) 

103  F0RMAT(I5,5(E14.6)) 

104  FORMAT (2(15)) 

105  FORMAT ( '  • ,A) 

110  F0RMATC//7X, 'TWO  FILES  ARE  REQUIRED  TO  HOLD  PROGRAM  OUPUT',//7X, 

1 ' DAT  A .  PLEASE  ENTER  THE  NAME  OF  THE  FILE  TO  BE  UStD  TO  *,//7X, 
2'HOLD  THE  PROGRAM  MATERIAL  PARAMETER  DATA.  PLEASE  ',//7X, 

3 ' INCLUDE  THE  EXTENSION.  I.  E.  FI L ENAME . DAT ' ,//7X ) 

111  FORMATC//7X, 'PLEASE  ENTER  THE  NAME  OF  THE  FILE  TO  BE  USED  ',//7X, 

1 ' TO  HOLD  THE  PROGRAM  OUTPUT  DATA.  PLEASE  INCLUDE  THE',//7X, 

2' EXTENSION .  I.  E.  FI L ENAME . DAT //7X) 

112  FORMATC//7X, 'THE  HEADERS  ALLOW  THE  USER  TO  IDENTIFY  THIS',//7X, 

1 ' SET  OF  DATA  FROM  ALL  OTHER  SETS.', 

2//7X, 'PLEASE  ENTER  HEADER  #1  (64  CHARACTERS  MAX)',//7X) 
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113  F0RMAT(//7X, 'PLEASE  ENTER  HEADER  12  (64  CHARACTERS  MAX)’,//7X) 

114  FORMAT (//7X, 'ENTER  THE  NUMBER  OF  LAYERS  C . L E .  5)',//7X) 

CLOSE(l) 

CL0SEC2) 

STOP 

END 
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APPENDIX  M 


SOFTWARE  SOURCES 


Ryan-McFarland  FORTRAN 

609  Deep  Valley  Drive 
Rolling  Hills  Estates,  CA  90274 
(213)  541  -  4828 

Microcompatibles,  Inc.  GRAFMATIC 
301  Prelude  Drive 
Silver  Spring,  MD  20901 
(301)  593  -  5151 

Jewell  Technologies,  Inc.  GRAFPLUS 
Seattle,  WA 
(206)  937  -  1081 

West  Coast  Consultants  CURVE  DIGITIZER 
4202  Genesee  Avenue,  Suite  309 
San  Diego,  CA  92117 
(619)  565  -  1266 

Microsoft  FORTRAN 

16011  NE  36th  Way 
Box  97017 

Redmond,  WA  98073  -  9717 

Prof.  M.  A.  Morgan 
Code  62Mw 

Naval  Postgraduate  School 
Monterey,  CA  93943 

LT.  E.  M.  Connolly 

c/o  J.  J.  Connolly 
208  Kent  Drive 
Greenville,  NC  27858 
(919)  752  -  2085 

OR 

c/o  C.  W.  Snell 
228  Windsor  Road 
Greenville,  NC  27858 
(919)  756  -  4918 
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